Nuclear Spins as Quantum Testbeds: 
Singlet States, Quantum Correlations, and 
Delayed-choice Experiments 

A thesis 

Submitted in partial fulfillment of the requirements 
Of the degree of 
Doctor of Philosophy 

By 

Soumya Singha Roy 
20083009 



IISERPUNE 

INDIAN INSTITUTE OF SCIENCE EDUCATION AND RESEARCH PUNE 



August, 2012 



Certificate 



Certified that the work incorporated in the thesis entitled ''Nuclear Spins as Quantum 
Testbeds: Singlet States, Quantum Correlations, and Delayed-choice Experiments" , 
submitted by Soumya Singha Roy was carried out by the candidate, under my super- 
vision. The work presented here or any part of it has not been included in any other 
thesis submitted previously for the award of any degree or diploma from any other Uni- 
versity or institution. 



Date Dr. T. S. Mahesh 



i 



ii 



Declaration 



I declare that this written submission represents my ideas in my own words and where 
others' ideas have been included, I have adequately cited and referenced the original 
sources. I also declare that I have adhered to all principles of academic honesty and 
integrity and have not misrepresented or fabricated or falsified any idea/data/fact/source 
in my submission. I understand that violation of the above will be cause for disciplinary 
action by the Institute and can also evoke penal action from the sources which have 
thus not been properly cited or from whom proper permission has not been taken when 
needed. 



Date Soumya Singha Roy 

Roll No.- 20083009 



iii 



iv 



Acknowledgement 



This thesis would have not been possible without the help and support by many 
individuals whom I met during my PhD life at USER Pune. I have been very privileged 
to have so many wonderful friends and collaborators. 

First of all, I am grateful to my advisor Dr. T. S. Mahesh for teaching me every- 
thing about NMR and Quantum Information Processing starting from the scratch. His 
cheerful guidance and deep understanding in NMR-QIP are the driving forces behind 
this thesis. Being his first PhD student, I consider myself to be very fortunate to get all 
of his support, care and aff"ection. 

I would like to thank Dr. Vikram Athalye for his valuable set of lectures which 
ultimately led to a fruitful collaborative work. I thank Prof. G. S. Agarwal for some 
constructive discussions with him and visiting us in spite of his busy schedule. I am 
thankful to Prof. Apoorva Patel for useful discussions and collaboration. Very special 
thanks to Prof. Anil Kumar for his great support and insightful conversations on various 
scientific problems and my future career path. 

I thank all the members of NMR Research Center -past and present- with whom I 
have worked. It's been a great pleasure working with Abhishek (whom we fondly call 
'Shuklaji') for all his 'complicated' queries and thoughtful discussions that I had with 
him. His famous 'I am not saying this, but I am saying that' is really unforgettable. 
It was highly exciting to work with Hemant who became an integral part of the NMR 
center ever since he joined the lab. Discussions on 'Quantum weirdness' were always 
intriguing with Manvendra. Working with Swathi was interesting because of her thor- 
ough theoretical understandings. Short stay of Philipp left many sweet memories of the 
discussions that I had with him. Discussions with Sheetal and Abhijeet were interesting 
and that actually made me learn many aspects of NMR and Quantum Computing. I also 
thank Pooja for her organized way of managing the spectrometers for years. Sachin was 
always there to look after the spectrometers and he made sure that it's up all the time. 

V 



I thank my RAC members- Dr. A. Bhattacharyay, Dr. T. G. Ajitkumar, and Dr. K. 
Gopalakrishnan for all their support. I also thank Dr. R. G. Bhat, Dr. V. G. Anand, Dr. 
H. N. Gopi, for all those help and affection. I am thankful to Prof. K. N. Ganesh for 
providing all the necessary experimental facilities in the lab. I thank Dr. V. S. Rao for 
all his help whenever I needed, be it academic or non-academic. I would like to thank 
CSIR and USER for the graduate scholarships that I received during my PhD. 

Life outside lab was also highly enjoyable for many many friends and I thank them 
all. I thank all my friends in Physics department, especially my 2008 batch mates for all 
their support and good times we had together. Arthur's casual approach and Murthy's 
'mass' approach created a funny contrast that I had thoroughly enjoyed all these years. 
Mayur always surprises me with his witty comments and jokes. Arun, Kanika, Padma, 
Ramya, and Resmi provided jovial company always. I thank all my friends in Chemistry 
department for helping me with many chemical compounds whenever I needed. 

I thank Anurag, Harsha, JP, and Amar for all those adventurous weekend treks and 
cheerful company. Dinner table was always full of noise, argument, and fun because 
of Biplab, Abhigyan, and Dada and I thank them all for making it so fascinating. Dada 
also has been my very close friend and room-mate for all these years. 

I thank all my long-time and long-distance friends- Sudipta, Diganta, Kalyan, Nan- 
dan, Swarup and Souravda for their support and encouragement. Conversations were 
always cheerful and motivating with them. 

My research career would have not been possible without the active support of my 
family. I thank my mother, father and sister for their love, support and encouragement 
throughout all the time. 



vi 



Contents 



Certificate 

Declaration 

Acknowledgement 

List of Figures 

List of Publications 

Abstract 

1 Introduction 

1 . 1 Nuclear Magnetic Resonance 

1.1.1 A nuclear spin under a static magnetic field 

1.1.2 Radiofrequency field 

1.1.3 Nuclear spin interactions 

1.1.4 Systems of spin- 1/2 nuclei 

1.1.5 NMR Relaxation 

1.2 Quantum Information Processing 

1.2.1 Computational science 

1.2.2 Quantum Information 

1.2.3 Quantum Bits 

1.2.4 Quantum Gates 

1.2.5 Quantum Algorithms 

1.2.6 Experimental implementations of QIP . . . 

1.3 NMR QIP 

vii 



Contents 



1.3.1 NMR- A suitable candidate for QIP 

1.3.2 NMRQubits 

1 .3.3 Initialization of NMR Qubits 

1.3.4 NMR Quantum Gates 

1.3.5 Numerically optimized quantum gates 

1.3.6 Measurement 

1.3.7 Coherence order 

1 .3.8 Limitations of NMR-QIP 

1 .4 NMR - An ideal platform for studying quantum mechanical phenomena 

2 Density Matrix Tomography of Long Lived Singlet States 

2.1 Introduction 

2.2 Long-lived singlet states 

2.2.1 Why singlet state is long lived ? 

2.2.2 Singlet Preparation in NMR 

2.3 Density Matrix Tomography 

2.4 Singlet State Characterization 

2.4.1 Observing through antiphase magnetization 

2.4.2 Tomography under varying spin-lock duration 

2.4.3 Offset dependence 

2.5 Long lived singlet states in multi-spin systems 

2.5. 1 Long lived singlet states in a 3-spin system 

2.5.2 Long lived singlet states in a 4-spin system 

2.6 Conclusions 

3 Preparation of Pseudopure States Using Long Lived Singlet States 

3.1 Introduction 

3.1.1 A pure state and a mixed state 

3.1.2 Necessity of Pure states in QIP 

viii 



Contents 



3.1.3 Pure states in NMR 

3.1.4 Pseudopure states 

3.2 Methods for preparing pseudopure state 

3.2.1 Temporal averaging 

3.2.2 Logical labeling 

3.2.3 Spatial averaging 

3.3 Preparation of pseudopure states using Long-Lived Singlet States .... 

3.3.1 preparation of singlet states 

3.3.2 Initializing NMR Registers 

3.4 Experiments 

3.4.1 2-qubit register 

3.4.2 3-qubit register 

3.4.3 4-qubit register 

3.5 Conclusions 

4 Storing Entanglement Via Dynamical Decoupling 

4.1 Introduction 

4.1.1 Decoherence 

4.1.2 Dynamical decoupling 

4.2 Uhrig dynamical decoupling 

4.2.1 Efficiency of UDD over CPMG 

4.3 Preparation of Entanglement HOC 

4.3.1 Preparation of singlet states 100 

4.3.2 Preparation of other Bell states from singlet states 102 

4.4 Storage of entanglement by UDD IM 

4.4.1 Different orders of UDD 

4.4.2 Performance of UDD over CPMG sequence 106 

4.4.3 Decay of magnetization during various dynamical decouplings . Il08 

ix 



Contents 



4.4.4 Efficiency of UDD over CPMG for a non-entangled state and 

various Bell states 

4.5 Conclusions 

5 Violation of Leggett-Garg Inequality 

5.1 Introduction 

5.2 Leggett-Garg inequality 

5.2.1 Spin- 1/2 precession 

5.3 Evaluating TTCCs using network proposed by Moussa et al 

5.4 Experiment 

5.4.1 Confirmation of dichotomic nature of x-component of nuclear 
spin observable 

5.4.2 Violation of LGI for 3 measurement case 

5.4.3 Violation of LGI for 4 measurement case 

5.5 Conclusion 

6 Quantum Delayed- Choice Experiment 

6.1 Introduction 

6.2 Studying wave-particle duality by interferometers 

6.2.1 Mach-Zhender Interferometer 

6.2.2 Wheeler's delayed-choice experiments 

6.2.3 Quantum delayed-choice experiments 

6.3 Theory 

6.4 Experiment 

6.4.1 Open and closed interferometers 

6.4.2 Quantum delayed-choice experiment 

6.5 Conclusions 

A Density Matrix tomography for a pair of spin- 1/2 homonuclear system 



X 



Contents 

B Density Matrix tomography for a three spin- 1/2 homonuclear system 153 
Bibilography Il57 



xi 



List of Figures 



1 . 1 Bloch sphere representation 12 

1.2 Quantum gates 22 

1 .3 Switching ON and OFF of interactions for the first qubit |31 

2.1 spectrum of 5-bromothiophene-2-carbaldehyde 52 

2.2 The pulse sequences for the preparation of singlet states and detection via |53 

2.3 Data characterizing the singlet state under CW spin-lock at an RF am- 
plitude |55 

2.4 Data characterizing the singlet state under WALTZ- 16 (CPD) spin- lock 

at an RF amplitude [56 

2.5 Bar plots showing traceless part p., of the theoretical singlet state density 
matrix |57 

2.6 Correlations calculated using the density matrix tomography of singlet 
states prepared with different spin-lock conditions: [59 

2.7 Pulse sequence for the creation of long lived singlet states in a 3 spin 
system ^ 

2.8 Experimental results of 3-spin LLS [61 

2.9 Pulse sequence for the creation of long lived singlet states in a 4 spin 
system ^ 

2.10 Experimental results of 4-spin LLS [63 

3.1 Bloch sphere representation of pure states, mixed states |67 

xiii 



List of Figures 



3.2 Simulations of Grover's search algorithm 69 

3.3 Population distribution of a two-spin system at room temperature . ... 72 

3.4 Preparation of |00) pseudopure state by using temporal averaging method 73 

3.5 Preparation of pseudopure state by using logical labeling technique . . . [75 

3.6 Pulse sequence for preparing |00) pseudopure state by using spatial av- 
eraging method 76 

3.7 Pulse sequence for the preparation and detection of singlet states .... 79 

3.8 Pulse sequence for the creation of |01) pseudopure state 80 

3.9 Circuit diagram for the preparation of a 3 qubit pseudopure state .... 81 

3.10 Circuit diagram for the preparation of a 4 qubit pseudopure state .... 83 

3.11 The circuit diagrams for initializing n-qubit register 84 

3.12 Experimental results for a 2 qubit register 86 

3.13 Experimental results for a 3 qubit register 88 

3.14 Experimental results for a 4 qubit register 89 

3.15 Pictorial description of the evolution of PPS over time |90 

4.1 NMR spectrum and the molecular structure of 5-chlorothiophene-2- 
carbonitrile 

4.2 NMR pulse sequence to study dynamical decoupling on Bell states . . . 

4.3 Density matrix tomography of Bell states 

4.4 Pulse sequences for various orders of Uhrig Dynamical Decoupling . . 

4.5 Experimental correlations (circles) of singlet state as a 

4.6 correlation exceeding 0.9 for 

4.7 The decay of the singlet spin-order measured by 

4.8 Experimental correlations of the product state and various Bell states . . 



5.1 The protocols for evaluating ^3 = Cn + C23 - C13 . . 

5.2 Quantum network for the evaluation of TTCCs . . . 

5.3 The energy level diagram NMR spectra of and '^C 

5.4 Intensity of decoupled ^^C spectrum as a function 

xiv 



101 



101 



103 



105 



107 



108 



109 



111 



117 



119 



122 



123 



List of Figures 



5.5 Correlations versus 

5.6 Decay of ^3 w. r. t. time 

5.7 The individual correlations C12, C23, C34, and 

6.1 Different types of Mach-Zehnder interferometer setups 

6.2 Molecular structure of chloroform (a) and pulse- sequences 

6.3 The experimental spectra obtained after the open 

6.4 The experimental intensities 5 p_o (particle) and 5 w,o (wave) 

6.5 The experimental spectra obtained after the quantum delayed choice 
experiment 

6.6 The intensities 5 wp,o(Q^> <P) versus phase 



XV 



List of Publications 



1. S. S. Roy and T. S. Mahesh, 

Density Matrix Tomography of Singlet States, 
J. Magn. Reson. 206, 127 (2010). 

2. S. S. Roy and T. S. Mahesh, 

Initialization ofNMR Quantum Registers using Long-Lived Singlet States, 
Phys. Rev. A 82, 052302 (2010). 

3. S. S. Roy, T. S. Mahesh, and G. S. Agarwal, 

Storing Entanglement of Nuclear Spins via Uhrig Dynamical Decoupling, 
Phys. Rev. A 83, 062326 (201 1). 

4. V. Athalye, S. S. Roy, and T. S. Mahesh, 

Investigation of Leggett-Garg Inequality for Precessing Nuclear Spins, 
Phys. Rev. Lett. 107, 130402 (2011). 

5. S. S. Roy, A. Shukla, and T. S. Mahesh, 

NMR Implementation of Quantum Delayed-Choice Experiment, 
Phys. Rev. A 85, 022109 (2012) 

6. S. S. Roy and T. S. Mahesh, 

Study ofElectromagnetically Induced Transparency using Long-Lived Singlet States, 
arXiv : quant-ph/1 103.3386 

7. H. Katiyar, S. S. Roy, T. S. Mahesh, and A. Patel, 

Evolution of Quantum Discord and its Stability in Two-qubit NMR Systems, 
Phys. Rev. A 86, 012309 (2012) 

8. S. S. Roy, M. Sharnia, V. Athalye, and T. S. Mahesh, 
Experimental Test of Quantum Contextuality in Nuclear Spin Ensembles, 
(in preparation). 



xvii 



Abstract 



Nuclear Magnetic Resonance (NMR) forms a natural test-bed to perform quantum 
information processing (QIP) and has so far proven to be one of the most successful 
quantum information processors. The nuclear spins in a molecule are treated as quantum 
bits or qubits which are the basic building blocks of a quantum computer. 

The long lived singlet state (LLS) has found wide range of applications ever since it 
was discovered by Carravetta, Johannessen, and Levitt in 2004. Under suitable condi- 
tions, singlet states can live up to minutes or about many times of longitudinal relaxation 
time constant (T 1 ) . For the first time, we have exploited the long lifetime of singlet states 
in NMR to execute several potentially important QIP problems. We were able to pre- 
pare high fidelity pseudopure states (PPS) in multi-qubit systems starting from LLS. We 
developed an efficient scheme of density matrix tomography to study all these quantum 
states. The tomographic study on LLS shows some interesting results. We performed 
experiments, where we created all the four Bell states from LLS and then studied the 
effect of various dynamical decoupling sequences on preserving these states. We found 
that Uhrig dynamical decoupling sequence is better than CPMG sequence in preserving 
Bell states for longer duration under suitable conditions. 

Nuclear spin systems form convenient platforms for studying various quantum phe- 
nomena. We used violation of Leggett-Garg Inequality (LGI) in a two-qubit system to 
study the transition from quantum to macrorealistic behavior. We observed perfect vio- 
lation of LGI for time scales which are much small compared to the spin-spin relaxation 
time scales. However, with the increasing time scales, we notice a gradual transition of 
spin-states from quantum to classical behavior. This steady arrival of classicality can be 
attributed to the decoherence process. In a separate experiment we performed quantum 
delayed choice experiment in nuclear spin ensembles to study the wave-particle duality 
of quantum states. These set of experiments clearly demonstrate a continuous morphing 
of the target qubit between particle-like and wave-like behaviors, thus supporting the 
theoreticians' demand to reinterpret Bohr's complementary principles. 
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Chapter 1 



Introduction 

"If we cannot possibly reach our desired destination, what is the point 
in setting out? This question might seem reasonable, but its premise is too 
restrictive: sometimes one walks not to reach a destination, but to observe 
the scenery along the way, and the pursuit ofNMR quantum computation 
has thrown up some surprising sights." 

- Jonathan A. Jones, 2010 

1.1 Nuclear Magnetic Resonance 

There are four main physical properties in an atomic nucleus: mass, electric charge, 
magnetism, and spin [[H |2l. Most of the macroscopic physical or chemical properties 
of matter depend on the mass and charge characteristics of nucleus. Though it is less 
evident, most of the nuclei are magnetic and behave like a tiny bar magnet [HI. However, 
this nuclear magnetism is very weak and may have little consequence on the matter's 
property. The dynamics of a nuclear spin can not be understood fully under classical 
physics and one has to invoke quantum mechanics. The spin and the associated nuclear 
magnetism provide us the tool to look not only inside the atom but also its microscopic 
world [HI El. 
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The first direct evidence of nuclear magnetism was given by Stem and Gerlach in 
1922 [[3]|- The Stern-Gerlach experiment involves sending a beam of particles through an 
inhomogeneous magnetic field and observing their deflection [3|. Much to the astonish- 
ment of classical physics, the beam splits into only two parts depending on the parallel 
and anti parallel alignment of their respective magnetic moment in the magnetic field. 
The exact measure of proton's magnetic moment was done by a series of experiments 
performed by Frisch, Estermenn, and Stern during 1933-1937 [SI Almost around 
the same time Isidor Rabi was working on the nuclear magnetism using the extended 
version of the Stern-Garlach apparatus. Rabi and co-workers showed the first indication 
of 'nuclear magnetic resonance' in molecular beams Soon after, this resonance ef- 
fect achieved its spectroscopic importance after Bloch |[8l|9l[lQl and Purcell [[m[T2l[T3]| 
independently observed nuclear magnetization in a bulk matter in 1946. 

Since then, NMR has been studied extensively and has found wide field of applica- 
tions in physical, chemical, biological, medical and material sciences. In this section, 
we give a brief overview of the basic principles of NMR. Later sections will describe 
the field of quantum information processing and its physical realization through NMR. 

1.1.1 A nuclear spin under a static magnetic field 

Let us consider the simplest situation where we have a single nucleus (an isolated spin) 
placed in an external magnetic field Bq. The magnetic nuclei will have a characteris- 
tic 'Larmor' frequency of (jJq = -yBo. Here y represents gyromagnetic ratio of the 
particular nuclear isotope. The Zeeman Hamiltonian can be written as 

-H, = -fiB 

= -MzyBo = najolz, (1-1) 

where fi is the nuclear magnetic moment operator and the external magnetic field Bq is 
taken along z direction. denoting the z component of the nuclear spin operator and the 
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relation between the spin operator and magnetic moment can be written as /x = yfilz- 
Below is a table showing comparative study of different properties relevant in NMR for 
various nuclei O. 



Nucleus Spin Natural Gyromagnetic ratio NMR frequency at 1 1 .7 T 
abundance(%) y/ 10^ rad s-'T"' {ooo/ln) in MHZ 





1/2 


-100 


267.522 


-500.000 




1/2 


1.1 


67.283 


-125.725 


19p 


1/2 


-100 


251.815 


-470.470 


31p 


1/2 


-100 


10.394 


-202.606 



Table 1.1: Table of most commonly used nuclear isotopes in NMR 



The eigenvalues of the Zeeman Hamiltonian (Eq. 11.11) represent the energy levels of 
the nucleus and are given by 

Em = -mhojo. (1.2) 

Here m represents the magnetic quantum number and it can take certain discrete values 

m = -I, -I + 1 , , / - 1,1, where / can be integer or half-integer and is known as the 

spin quantum number. /(/ -I- l)tf is the eigenvalues of total spin operator I^. 

While Iz represents a stationary state under the Zeeman Hamiltonian, (I^) and (ly) 
show out of phase oscillations at Larmor frequency (coq)- Higher (positive) m values 
have lower energy state (Eq. 11.21) and thus the ground state is the state with m = I. 
In a semiclassical picture, it can be seen as the nuclear spin that is aligned along the 
static magnetic field direction. On the other hand, highest excited state corresponds to a 
spin- alignment against the magnetic field. 

For an ensemble of nuclear spins at thermal equilibrium , the population distribution 
can be represented by Boltzmann statistics. For spin-1/2 ensemble, there will be only 
two possible energy levels with m = -1/2 and m = -1-1/2. The population ratio of these 
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two levels is detemiined by Boltzmann distribution 

£_ ^ ^-ncoo/ksT ^ (1.3) 

where is the Boltzmann constant and T is the absolute temperature of the ensemble. 
In the case of nuclei at a 10 T magnetic field strength, Hcoq ~ lO'^eV. Whereas 
at room temperature ksT « 2.5 x 10~^eV, hence the ratio Hcoo/kBT « 10"^. So the 
Boltzmann factor g ^''^o/^^sT' almost close to unity. This can naively be interpreted as, 
there are slightly more spins in the parallel direction (lower state) than in the anti-parallel 
direction (upper state) and this slight imbalance in the populations is responsible for the 
'net' nuclear magnetization along the z-direction. This also reveals the fact that, NMR 
is a very low sensitive technique. 

The nuclear magnetization for an ensemble of spin- 1/2 nuclei at thermal equilibrium 
is given by lfT4l 

nny^Jf'Bo 

Mo = — -, (1.4) 

where no is the number of nuclei per unit volume. From above equation it is clearly 
seen that the magnetization increases linearly with the external field strength, whereas it 
is inversely proportional to the temperature. Hence nuclear magnetism is paramagnetic 
in nature and follows Curie's law [[T4ll . Also, the demand for higher field strength can 
be understood from the above equation. However, the temperature of the ensemble can 
not be reduced as per wish, since it is related to the 'state' of the matter and hence 
on its dynamics. Here it can be noted that, electrons also posses paramagnetism and 
the magnitude of electron paramagnetism is three order of magnitude higher than the 
nuclear magnetism. 



1.1.2 Radiofrequency field 

The application of static magnetic field will create a Zeeman splitting according to the 
Eq. 11.11 Now the transitions between the energy levels can be induced by the application 
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of suitable oscillatory magnetic fields with appropriate frequencies. From the Table 
ll.li it is seen that the Larmor frequencies are of the order of MHz in present days' 
magnet of a few Tesla and resonance can be achieved by the application of RF fields. In 
Comparison, typical electron Larmor frequencies are of the order of GHz range. 

The dynamics of nuclear spin excitation due to the application of oscillatory mag- 
netic field can be well understood by considering a time dependent magnetic field, Bi{t) 
applied perpendicular to the static magnetic field Bq. The RF interaction Hamiltonian 
CHrf), can be written in a similar way as the Zeeman Hamiltonian. 



Here Q and ^ are respectively the frequency and phase of the RF field which is along 
the X direction. The strength of the oscillatory magnetic field (Bi{t)) is much smaller 
than the Zeeman field strength (Bq) an hence it is reasonable to treat the RF Hamilto- 
nian CHrf) as a perturbation to the Zeeman Hamiltonian CKj). The dynamics can be 
described by the standard time dependent perturbation theory [fTSll . The result shows 
that, at resonance condition (Q = ojo), there will be induced transitions between the 
eigenstates of 'H, with a transition rate given by the Fermi golden rule [ 16j 



where nii and m2 are two energy eigenstates of the system. It can be seen from the 
above equation that, the transition probability on either way depends on square root of 
gyromagnetic ratio of the nucleus and the magnitude of RF field. The selection rule for 
the allowed transition should be. Am = + 1 . 

Now we will discuss the logic behind choosing the RF magnetic field similar to 
Eq. II. 6[ We can think of a linearly polarized magnetic field Bi{t) as composed of two 
circularly polarized fields with same frequency and amplitude but precessing in opposite 



•Krf = -fi.Biit) = -yhlAlBx cos(Q? + 0)] 



(1.5) 



where, Bi{t) = 2Bi cos(Q.t + (p)x 



(1.6) 




(1.7) 



5 



Chapter 1. Introduction 



directions about z-axis. 



Bi{t) 



2Bi cos(Q.t + (p)x 



Blit) + B-it) 



(1.8) 



where, Bl(t) 



Bi[cos(Q.t + 



+ 



sm(Clt + (f>)y] 



(1.9) 



B^(t) 



5i[cos(Q? + 0)jc 



sm(Qt + (f))y] 



(1.10) 



The RF field interactions can be better described in rotating frame formalism. At res- 
onance condition, (i. e. £1 = coq) the field B^{t) rotates coherently with the nuclear 
Larmor precession along z-axis. Whereas, the field B'l^(t) rotates exactly in opposite 
sense. In a frame which is rotating along with the Larmor frequency, the field B^(t) 
is stationary along with nuclear spin, whereas the field B^(t) rotates with a frequency 
twice the Larmor frequency. Therefore, at high static fields it can safely be assumed that 
only the field By{t) has efi'ect on the nuclear spins. 

Let us assume the on-resonance condition i.e., Q. = ojq. In a frame that is rotating 
with -Bj^(0 with same frequency and direction, the magnetic moment sees a static field, 
say along direction x', and precesses about it. In the case of ofF-resonance conditions 
(i.e. Q. 4^ a>o), the precession of magnetic moments in the rotating frame is around an 
axis defined by an eff'ective magnetic field given by. 



At on-resonance condition, the precession frequency about B^ff is also known as nu- 
tation frequency a»„„, = -yB^ff in analogy with the Larmor frequency. Application 
of an RF pulse for the time duration tp, makes the magnetization shift from its initial 




(1.11) 



The relation between laboratory frame and rotating frame is given by. 



x' = cos(Q.t + 4>)x - sin(D? -I- (p)y 



(1.12) 
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z-direction by a nutation angle given by, 



Op = yBitp. 



(1.13) 



Hence, a n/2 pulse is defined as a pulse which can take the magnetization from longi- 
tudinal direction to transverse plane. One must remember that in laboratory frame the 
magnetization is always precessing around the z- axis in addition to nutating about the 
RF axis. 

1.1.3 Nuclear spin interactions 

So far we have described the nuclear spins in isolated situation without any kind of 
interactions. In practice, nuclear spins are interacting with each other as well as with 
the environment. The interaction of nuclear spins with each other makes NMR a very 
sophisticated tool with versatile applications. However interaction of nuclear spins with 
environment remains a challenge in the field of NMR-QIP and we will discuss this case 
in detail in a later chapter. Here we describe the main interactions involving in the 
nuclear spins under normal conditions [[T6l[T7]| . 

The total nuclear Hamiltonian is given by 



where "Hint represents the internal interactions of the nuclei. Here we will concentrate 
on the l-fint part of the total Hamiltonian. There are several contributors to the internal 
Hamiltonian part based on its physical and chemical characters. In most of the case the 
material in study under NMR is a diamagnetic insulating substance. For this the internal 
Hamiltonian is given by 



(1.14) 



(1.15) 
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where, 'Hcs is the chemical shift interaction, 'Kd is the direct dipolar interaction, 'Hj is 
the indirect spin-spin interaction, and 'Hq is the quadrupolar interaction. 

Chemical Shift 

Though the external magnetic field applied is same for all the nuclei, it is not even 
exactly same for a same type of nuclei in a molecule. The slight change in the mag- 
netic field is due to the modified chemical environment created by the electron density 
surrounding it. The modified magnetic field is given by 

Bloc = {\ -a)Bo, (1.16) 

where a is known as chemical shielding tensor allied to that particular nuclear site. 
Hence the chemical shift Hamiltonian can be written as, 

"Hcs = -M • (-5=5,,) = yncT,,{e, moh- d-^) 

The approximation is known as secular approximation. Now, 

(T^^{6, (p) = (Tn sin^ 6 cos^ (p + 0-22 sin^ 6 sin^ (p + (T33 cos^ 6, (1-18) 

where, (Xi 1,0-22, and 1X33 are the principle values of the chemical shielding tensor cr. 
Here 9 and (p are the azimuthal and polar angle respectively, describing the magnetic 
field Bq in the principle axis system. In isotropic liquid, due to rapid molecular motions, 
shielding tensor get averaged. Hence, the time averaged shielding constant for isotropic 
liquid can be written as, 

O-isO = ^ {O-XX + 0-yy + Cr^^). (LI 9) 

The consequence of the above calculation is the introduction of a shift in resonance 
frequency, 

o) = coo{\ - o-i,o). (1-20) 
8 
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For a monocrystalline material, the above equation will be modified just by replacing 
(Tiso with cr^^. In the case of polycrystalline material or powder samples, the continuous 
distribution of orientations of the several crystallites causes an anisotropic broadening, 
known as chemical shift anisotropy (CSA)[|2l, 

i^a = a^^-]^{a^x + ayy). (1.21) 

The resonance frequency is conventionally expressed by the relative shift from the ref- 
erence resonance frequency {oJref), 

OJ — OJypf 

6 = —. (1.22) 

O^ref 

Here 6 represents the chemical shift of the resonance lines and normally expressed in 
terms of parts per million (ppm). 



Direct dipolar coupling 

Any two magnetic dipole moments interact directly with each other through the mag- 
netic fields created by each one for the others. It provides rich structural information 
about the materials. The dipolar Hamiltonian is defined as. 



k<i 



- E 



k<i 



4n -3 



' ki 



Ik ■ Ii- 3\ilk . ndili . ru) 
r 



(1.23) 
(1.24) 



where D],\ is the dipole-dipole interaction tensor, r^./ is the radius vector connecting the 
two spins. Under secular approximation, the Hamiltonian can be rewritten as. 



c^trunc 



-E 



k<i 



ykji^ 1 



An rl 2 



- (3 cos^ en - l) [3/fe4 - ik . lA , (1.25) 
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where Oki is the angle between rki and z- In case of heteronuclear spin systems (i.e. 
Jk yd, further simplification is possible, 

= ^^(i cos' eu)hJi.. (1.26) 

kl 

Indirect spin-spin coupling 

Indirect spin-spin coupling (also called J-coupling or scalar coupling) is also an interac- 
tion between the nuclear magnetic dipole moments. This type of coupling is not direct 
and being mediated by the electron cloud involved in the chemical bonds between the 
atoms. The J-coupling Hamiltonian is defined as, 

nHj = 2nnY^hJI,, (1.27) 

k<i 

where J is the J-coupling tensor. J-coupling posses an isotropic part which survives 
under random molecular motion in an isotropic substance (e.g. liquid samples), whereas 
direct dipolar coupling is averaged out under similar situation. In the case of solid 
samples, the J-coupling is generally overwhelmed by the strong direct dipolar couplings. 
Under secular approximation, the simplified J-coupling term is written as, 

= 2nnJuIkJu. (1.28) 

The approximation can be carried out when InJij « - (x)j\. It can be seen that this 
approximation holds for all heteronuclear pairs. 

Quadrupolar coupling 

All the nuclei with spin, I > 1/2 are subjected to electrostatic interaction with the neigh- 
boring electrons, ions due to the non-spherical charge distribution of nuclei tlSj. The 
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Hamiltonian form of quadrupolar interaction is defined as, 



N 



(1.29) 



where, Qk is the quadrupolar coupling tensor and it can be expressed in terms of electric 
field gradient tensor Vk at the k''^ nuclear site. 



Qk 



eQk 



24(27,-1) 



V,. 



(1.30) 



Here Qk is the nuclear quadrupolar moment of the k''' nucleus. In the partial axis coor- 
dinate system, the quadrupolar Hamiltonian for the k''' nucleus can be written as. 



"He 



44(27, - 1) 



(1.31) 



where eqk = and 77, defines the assymetry parameter. 



rik 



kyy 



Vk,. 



(1.32) 



1.1.4 Systems of spin- 1/2 nuclei 



The Hamiltonian for N coupled spin- 1/2 nuclear spins in an isotropic medium is given 

by, 

N 

"H^Yj ^'t^fc + Z '^''Jkilk.Iu (1.33) 

k=\ k<i 

where o), is the chemical shift for the k*'^ nucleus and 7,/ is the J-coupling constant 
between the two spins. Considering weak coupling condition i.e. \lnJki\ « - coj}, 
the Hamiltonian can be written as, 

N 

n< = Yj ^kh, + ^^Jkilkzh- (1-34) 

k=l k<l 
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The coupling part of the above equation actually commutes with the Zeeman part and 
hence both the part will share common eigenbasis. All the 2^ eigenstates of 'H can be 
expressed as tensor products of the single spin eigenstates, namely \aa . . .a), \aa . . ./3), 
. . . ,\/3/3 . . -13). Here la) and 1/3) denote | + 1 /2) and | - 1 /2) single-spin eigenstates, which 
are labeled as |0) and |1) in QIP terminology. The NMR spectrum displays A^^ set of 2^~^ 
spectral lines of equal intensity. 

The Hamiltonian for a pair of spin- 1/2 system in an isotropic liquid environment can 
be written as, 

"Hki = cokhz + ojiliz + InJkihh, (1 .35) 

where we have considered the weak coupling condition. This Hamiltonian will have 
four eigenstates and corresponding four eigenenergy values. The four probable transi- 
tion will reflect as four transition line in an NMR spectra. The eigenstates and eigenen- 
ergy are: 

100) :: £oo= -{-co^- oJi + nJ) 12 

101) :: £oi = -{-0Jk + 0Ji-nJ)l2 

(1.36) 

110) :: £io= -{cou-coi-nJ)l2 

111) :: En= - (cok + coi + nJ) /2 

1.1.5 NMR Relaxation 

In equilibrium, the population distribution of the spins follow Boltzmann statistics with 
off diagonal elements are zero for the density matrix of the system. The NMR mecha- 
nism depends on the perturbation of the system from equilibrium situation. For example, 
application of a single n/2 pulse on equilibrium equalizes the populations and also cre- 
ates the coherences. Now, this is clearly a non-equilibrium situation and the disturbed 
state tends to go back to the original equilibrium state through relaxation mechanism of 
the spins. There are two different processes, occurring simultaneously but in general in- 
dependently that can be identified for this relaxation. These two relaxation mechanisms 
known as transverse relaxation and longitudinal relaxation ffT9l[T6]| . 
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Just after the RF pulse, the magnetization is on the transverse plane perpendicular 
to the static magnetic field Bq. The transverse relaxation mechanism makes the mag- 
netization along transverse plane to disappear. The result of the transverse relaxation 
is the loss of coherences among the spins. This happens due to the spread in nuclear 
precession frequencies of the spin ensemble. As shown earlier, the Larmor frequency 
of each spin depends on the external magnetic field as well as locally created magnetic 
field for various reasons. Hence due to this slight variance in the Larmor frequency, 
after some time these spins are oriented in a complete random direction on transverse 
plane and the vector sum of all this magnetization will be zero. The decay of coherences 
due to the inhomogeneous fields is one part of the transverse relaxation process. The 
other important part occurs due to the fluctuations in the local magnetic field |I21 . 

Under normal conditions, the decaying of the transverse component of the mag- 
netization of the nuclear spins ensemble in the rotating frame can be described by a 
phenomenological differential equation given by Bloch |I91I3- 

dt 12 

where T2 is known as transverse or spin- spin relaxation constant. The solution of the 
above equation is simple and can be written as, 

M,,, = Moe-'/^^ (1.38) 

where Mq represents the initial value of the transverse magnetization. Hence from the 
above equation it is seen that the transverse magnetization decays with time in exponen- 
tial fashion. The exact value of T2 depends on the detail of each particular nuclear spin 
system and its environment. 

The longitudinal part of the nuclear magnetization also goes under relaxation simul- 
taneously with transverse relaxation. The mechanism can be understood as follows. Just 
after the nil pulse, the longitudinal magnetization, Mj. = and the population of a two 
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level system is equalized. Since this condition is non-equilibrium, the system will tend 
to go back to its equilibrium condition that is supported by Boltzmann distribution. The 
preferable way towards the equilibrium is by giving up its excess populations in upper 
to lower energy level till the Boltzmann distribution is reestablished. Since this mecha- 
nism involves energy exchange and that happens with the lattice part of the system, this 
relaxation mechanism is also termed as spin-lattice relaxation process [|2l. 

Similar to the transverse case, the longitudinal relaxation mechanism is also de- 
scribed by a phenomenological differential equation given by Bloch |l9l|2l. 



where Ti representing the longitudinal or spin-lattice relaxation constant. The solution 
of the above equation is given by. 



As it can be seen from the above solution, the longitudinal magnetization is gaining 
with time beginning from zero and reaches to the stable magnetization Mq after certain 
time. The exact values of T\ and T2 time constants depend on various factors such 
as physical state of matter (liquid or solid), temperature, molecular mobility, viscosity, 
concentration, external magnetic field etc |l2l. In most of the cases it is found that 
Ti > T2. In case of liquids, T2 values are comparable with and in many cases both 
are almost equal. However in case of solids, Ti is much larger than T2. 

It is worth noting that the above simplistic approach of relaxation formalism in nu- 
clear spins is not straightforward in many complicated situations. The relaxation phe- 
nomenon can be best understood by the elaborative mechanism of Redfield theory lfT4ll . 



dM, Mo - M- 



(1.39) 



dt Tx 




(1.40) 
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1.2 Quantum Information Processing 

Quantum information processing (QIP) is the study of the information processing tasks 
that can be accomplished using quantum mechanical systems [|20l . The idea of utilizing 
quantum systems for the information processing was first introduced by BeniofF in early 
1980s ||2T1I221 . The exponential time required for simulating the dynamics of quantum 
systems using classical computers inspired Feynman to propose exploiting quantum 
systems for such purpose [23J. He was rather skeptical whether a classical computer 
is capable enough to simulate a quantum system and advocated for building a quantum 
computer for this purpose. In 1985 Deutsch gave a decisively important step towards 
quantum computers by presenting the first example of quantum algorithm which utilizes 
the fact of quantum superposition in speeding up computational process [|24ll . He is also 
the pioneer of quantum computer history for introducing the notion of quantum logic 
gate in 1989 [|25l . Since then there has been a good theoretical progress in the field of 
quantum computation and quantum information. Classically intractable problems were 
reduced to tractable regime by treating it in quantum way. It was 1994, when a major 
breakthrough happened, calling the attention of scientific community for the potential 
practical importance of quantum computation and its direct consequence on our society. 
Peter Shor discovered a quantum algorithm which is capable of factorization of prime 
numbers in polynomial time instead of exponential time |[26l l271l . Prime factorization 
being the heart of computational security, draws tremendous attention from computer 
scientists and cryptographers as well. The successful experimental tests of Shor's algo- 
rithm have been performed using a liquid state NMR systems Il28ll . A few years after 
that, in 1997, another important discovery had been made by Lov Grover by introducing 
a quantum search algorithm for searching an unsorted database f29l|. Grover's algorithm 
makes use of quantum superposition and quantum phase interference to find an item in 
an unsorted database, faster than any other classical algorithms. Various schemes on 
error correction has also being developed to counter the faulty outcomes [|30l . In the 
meantime, other branches of QIP, namely quantum teleportation, quantum key distri- 
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bution and quantum cryptography are also being developed. Many of these techniques 
have actually making commercial success and continue to be better OTI . Considering 
the extreme difficulty in controlling a quantum system, there has been modest develop- 
ment towards a practical quantum computer. Nonetheless, commercialization of quan- 
tum computer has been taken very seriously and till date it has already arrived (arguably) 
in the markets |[32]| . This section intend to give a brief theoretical understanding on QIP 
and later its physical realization by various experimental schemes. 

1.2.1 Computational science 

Today, we can not even think a society without the machine called computer. The im- 
pact of a computer is such that, there hardly any field left where we are not using a 
computer directly or indirectly. There is a long history of development of computers 
and the theoretical notion of computation. As put by David Deutsch [[33l . 'Computers 
are physical objects, and computations are physical processes. What computers can or 
cannot compute is determined by the laws of physics alone, and not by pure mathe- 
matics'. Computation is carried out through a procedure called algorithm and it needs 
three basic resources (space, time, energy) for it COl . Space refers to the the computer 
hardware, i.e. the number of logic gates used. Time refers to the computational time 
required and energy refers to the energy spent for the computational work. The basic 
model of a modem day computer was mainly given by Alonzo Church and Alan Turing 
in early 20th century. Later it became famous as Turing machine 0411 . A Turing ma- 
chine is a hypothetical, idealized theoretical model of an actual computer. There is not a 
single computation work which can be done by an actual computer but not by a Turing 
machine. In that sense, a real computer is a physical realization of a Turing machine. It 
consisted of a program, a finite state of control, a memory tape, and a read-write head 
EOl . The Church-Turing thesis calls a problem 'computable' only if it can be done by 
a Turing machine. Quantum computation also obeys the ideology of the Church-Turing 
thesis and hence the notion of 'computable' has not changed, only efficient algorithms 
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could be possible. The efficiency of an algorithm is studied by its asymptotic behavior 
as the size of the input increases [|20ll . Consider the time taken by an algorithm varies as 
f{N), where N is the number of input bits. Now, if f{N) is polynomial, then the highest 
power in f{N), say ^(A^) is known as order of algorithm denoted by O {g{N)). Depend- 
ing on these requirements, computational problems are classified into various classes 
known as 'complexity classes' as shown in table 11.21 A simple addition, subtraction, or 



Class Time Space 



EXP exponential unlimited 

PSPACE unlimited polynomial 

NP exponential polynomial 

P polynomial polynomial 

L logarithmic polynomial 



Table 1.2: Complexity classes in computational science. These classes are related as : 
L c P c NP c PSPACE c EXP 

multiplication are in class L, while division comes under class P. Prime factorization is 
believed to be a class NP problem, however not proven till date. Many of the complexity 
classes are unclear even today. In fact it is a great source of debate whether N = NP or 
N NP and nobody has come up with a concrete prove so far. 

The relationship of energy with information processing has an important physical 
significance OOll . Erasure of information is a dissipative process, as pointed out by Rolf 
Landauer in 1961 [ 35|. Erasure of each bit increases the amount of entropy by fcln 2 and 
the energy dissipates at least by an amount kT In 2. However, this amount is negligible 
compared with the energy dissipated in a modern computer which is of the order of 
SOOkT In 2. All the irreversible gates involve in loss of information and hence dissipates 
energy. Interestingly in 1973, it was found by Charles Bennett that the dissipation of 
energy can be made vanishingly small by making all the gates reversible Il36ll . 
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1.2.2 Quantum Information 

Information always exists as encoded with an physical system and therefore it should 
obey the physical laws. In other words, 'Information is Physical' Il30l [35]| and physi- 
cal systems obey quantum mechanics. Hence the information encoded in such system 
is 'Quantum Information'. Treating some problems in quantum mechanical way can 
actually make it much more efficient than classical way. For example, prime factoriza- 
tion is a 'NP' class problem classically (require exponential time), whereas solving it in 
quantum mechanical way can make it a 'P' class problem (require polynomial time). 

1.2.3 Quantum Bits 

The unit of information in quantum computation and quantum information is known 
as quantum bit or 'qubit'. A qubit can assume a logical values '0' and '1' along with 
a state that is a linear combination of them. Physically a qubit can be represented by 
any well defined distinct eigenstates. For example, qubits can be the polarization states 
of a photon or nuclear spins inside a static magnetic field. Let us consider a two level 
quantum system, where the eigenstates are represented by |0) and |1). The general form 
of a quantum state under this condition can be written as. 



where Q < 9 <n and < < 2n, neglecting the global phase factor. On measurement in 
|0), |1) basis, the probability of getting the state |0) is cos'^ie/l) and for |1) it is sin^(^/2). 
Also this kind of representation allows one to visualize this complex quantum state 
geometrically. The qubit states are designated as some geometrical point on the surface 
of a 'Bloch sphere' (Fig. 11.11) . Any surface point on the Bloch sphere is a 'pure' state 
while any non-surface point represents a 'mixed' state. A more detailed description 
about pure and mixed states is given in chapter 3. The power of quantum computation 
comes from the quantum mechanical laws such as superposition of states of qubits and 




(1.41) 
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the ability to manipulate the quantum states through unitary transformations as will be 
seen in next subsections. 



1.2.4 Quantum Gates 

A classical gate can transform a string of n bits into a string of m bits, 

/:{0,ir ^{0,ir. (1.42) 

Now for / to be a reversible classical gate, it should be one to one (each inout is mapped 
to a unique output). In general n and m are not equal and hence a classical gate is a 
irreversible gate. Quantum gates on the other hand transform a state of quantum system 
from one point in the Hilbert space to another point. A single qubit can be expressed 
by \if/) = a\0) + b\\), where a and b are coefficients having a relationship jap + = 1. 
Quantum gates on a qubit must preserve this normalization condition and thus it can be 
described by a 2 x 2 unitary matrices ll20l . Since all the quantum operations are unitary 
operators, quantum gates must also be reversible. 

19 



Chapter 1. Introduction 



Some important unitary transformations for one qubit are the Pauli matrices, 
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(1.43) 



These three Pauli matrices along with the 2x2 identity matrix form a 2 x 2 basis matrix 
space. Hence any one qubit operation can be decomposed as a linear combination of 
the four matrices. A NOT gate is nothing but the Pauli-x matrix and it flips the |0) to 1 1 ) 
and vice versa. 
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(1.44) 



(1.45) 



(1.46) 



Another very important one qubit gate is Hadamard gate which has no classical analogue 
and it is used for the creation of superposition states as shown below. One important 
property of Hadamard operator is its self -reversibility, i.e. H^ = 1. 



Um = 
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(1.47) 

(1.48) 
(1.49) 



A phase shift gate P selectively introduces a phase to either of the qubit of a superposi- 
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tion state, 



Up 



( \ 
1 



e" 

(a|0) + Z7|l)) (a|0) + Z7y^|l)). 



(1.50) 
(1.51) 



For a two qubit system the dimension of Hilbert space is 4 x 4 and can be realized by 
tensor products among the one qubit states, 



{|0),|1)}®{|0),|1)} = {|00),|01>,|10),|11)}, 



(1.52) 



where. 



100) = 



1 



101) = 



^0^ 



110) = 







|11) = 



^0^ 



(1.53) 



The matrix representation for operators that act only on one of the qubits of a system 
of two qubit can be constructed by tensor product between one qubit operator and 2x2 
identity operator. 

(9«=(9®1; (9fo = 1(8(9 (1.54) 

Here O denoting the Pauli matrix operators. The above given scheme can be worked 
out for any number of qubits in a similar fashion. The most important two qubit gate 
is definitely the CNOT (or controlled not) gate. It can be proved that all the quantum 
operations necessary for quantum computation can be achieved using only CNOT and 
set of one qubit gates Il20l . In that sense CNOT is the universal quantum gate similar 
to NAND gate in classical counterpart. A CNOT gate has a control qubit and a target 
qubit. Depending on the state of control qubit, the status of target qubit get flipped while 
the control qubit remaining same. The operator form of CNOT gate whose control is 'a' 
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Figure 1.2: Quantum gates, (i) Hadamard gate acting on |0> qubit, (ii) CNOT gate, (iii) 
A general two qubit controlled- gate, where U can be any one qubit operator, and (iv) 
Toffoli gate. In the above circuits, the inputs are assumed to be individual basis states. 
If on the otherhand, inputs are in superposition, output may be entangled. 



and target is 'b' (and vice versa) can be written as, 
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(1.55) 



CNOT can also be represented by binary addition of two qubits, i.e. CNOT„|a, Z?) = 
\a,a © b) and CNOThki,b) = |<3 © b,b). Here the symbol © represents the addition 
modulo 2, for which OffiO = 0, 0©1 = 1,1©0 = 1, and 1 © 1 = 1. The application of 
CNOT gate on two qubit states has the following results, 



CNOTalOO) = 100), CA^Or„|01) = 101), 

CA^orjio) = 111), CA^orjii) = |io). 



(1.56) 
(1.57) 



The circuit diagram of a CNOT gate is shown in Fig. 11.21 For a three qubit system, 
TOFOLLI gate is the universal gate which is nothing but a controUed-CNOT gate. 
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1.2.5 Quantum Algorithms 

Quantum algorithm solves problems by exploiting the properties of quantum mechan- 
ics. An efficient algorithm will require minimum resources. Many quantum algorithms 
are much more efficient than any classical algorithms by exploiting the features like su- 
perposition and entanglement. The first quantum algorithm was given by David Deutsch 
and is known as Deutsch algorithm. This algorithm is capable of finding out whether 
a binary function of one qubit is 'constant' or 'balanced' in one go |[24l . The most 
powerful quantum algorithm till date is the prime factorization given by Peter Shor. 
Shor's algorithm can factor a number by exponentially faster than its classical version. 
Grover's search algorithm can search an unsorted database in polynomially faster than 
classical algorithm. 

Adiabatic quantum algorithm gains much attention due to its universality [37]. In 
most cases a quantum algorithm begins with a uniform superposition and ends with an 
eigenstate which is the desired result. Often it is found that the ground state of the final 
Hamiltonian CH/) is the desired answer, however it is not easy to find the answer. Now 
suppose we have a Hamiltonian 'Hi whose ground state can easily be found. Hence by 
evolving the system adiabatically from "TY, to 'Hf, one can reach the ground state of 'Hf 
and hence the desired result. One has to make sure that there is no crossover of the 
ground state with any other state and the evolution process is slow enough that there 
won't be any possible transition. 

1.2.6 Experimental implementations of QIP 

While there is a good amount of progress in the theoretical understanding of QIP, the 
physical realization of a quantum computer is proving extremely challenging. DiVin- 
cenzo laid out five criteria which must be fulfilled for a successful quantum computer 
architecture [|38]| . 

1 . Well defined qubits 

2. Ability to initialize 
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3. Universal set of quantum gates 

4. Qubit-specific measurement 

5. Long coherence time 

Also there are two more criteria that will be needed for quantum communication. Meet- 
ing all these criteria in a single experimental setup is a highly challenging task. Nonethe- 
less, various techniques have been proposed and is being explored for QIP tasks. All the 
techniques have there own advantages and some disadvantages. The major techniques 
available today are : 

1. Nuclear spins in NMR 

2. Trapped ions / atoms 

3. Photons 

4. Spins in semiconductor : 

— Quantum dots 

— NV centers in diamond 

5. Superconducting circuits 

The first two techniques deal with the mutual interaction of quantum particles (atomic 
nucleus, atoms, ions) and controlled by electromagnetic field. Polarization of photons 
can be treated as qubits and it is controlled by optical means. Quantum dots and NV 
center in diamond techniques utilize the much developed semiconductor field in minia- 
turization scale. The well defined 'phase' and 'flux' parameters can serve as qubits in 
a superconducting circuits. Apart from these techniques, there are few more interesting 
techniques which might get much attention in future due to its hybrid approach. These 
methods are exploiting the best features among the available techniques and intend to 
make out a optimized experimental setup. For example, nuclear spins have much larger 
coherence time, but nuclear magnetization is very faint. By transferring magnetization 
from electrons to nuclei, the above problem can be solved and thus integrating the NMR 
with the ESR technique [[39ll . Another approach is integrating NMR with Atomic Force 
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Microscopy (AFM) which is capable of measuring a single atom compare to bulk sam- 
ple measurement done by NMR HOl . A comparative study of all the key techniques is 
given in Table Oll4ni42ll43ll44l. 



1.3 NMR QIP 

Application of NMR for the physical realization of QIP adds one more feather to the 
much colorful NMR application field. Implementation of QIP by NMR was indepen- 
dently proposed by Cory et al f45] and Gershenfeld et al [|46ll in 1997. Since the criteria 
laid by DiVincenzo was fulfilled more or less by NMR, it became an automatic choice 
whilst all other techniques were slowly coming up. One more thing fueled the NMR- 
QIP initiative was the fact that many of the QIP experimental basics are routinely done 
in conventional NMR experiments [|47ll . For example, the selective inversion of popula- 
tions achieved in 1973 is described as a CNOT gate [|48ll . The inversion of zero quantum 
coherence takes the name as SWAP gate Il49ll50ll . However, NMR-QIP gains much of 
its attention after Cory et al and Gershenfeld et al independently showed the preparation 
of 'pseudopure state' in a liquid state NMR at room temperature. NMR-QIP in liquids 
containing small number of spins (preferably spin- 1/2) have been studied extensively 
and its proven to be an excellent testbed for a small scale quantum information proces- 
sor. Many complicated algorithms have been tested and verified. For example, Shor's 
factorizing algorithm has been tested till date only by liquid state NMR ||28]| . How- 
ever, scalability of liquid state NMR is an issue which hurdling the possibility of being 
an 'useful' quantum information processor in long-run. It is unlikely to get more than 
15-20 qubits unless some technological breakthrough occurs [47]. On the other hand 
Solid state NMR has the potential to become a reliable QIP architecture in future, since 
scalability issue and preparing 'true' ground state seems more realistic. Some aspects 
of NMR-QIP are discussed in the following. 
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NMR 


Trapped ions 


Photons 


Semiconductors 


Superconductors 


System 


Nucleus 


Atom 


Photon 


Atom, Vacancy 


Phase, Flux, 
Charge 


Maximum Qubits 
demonstrated 


12 (entangled) in liq- 
uids, >100 (correlated) 
in solids 


lU-lU Stored, 14 
(entangled) 


10 (entangled) 


1 (l^JJS), J(JNV 

centers) 


izo (labncatea), 
3 (entangled) 


Coherence time 


>ls (liquids), ~ 100ms 
(solids) 


>ls 


~ IOOa/s 


1-10 yUS (QDs), 1- 
10 ms (NV) 


~ lOjUS 


Two qubit gates 
(highest fidelity) 


CNOT (>99%) 


CNOT (>99%) 


CNOT 

(>94%) 


-90% (NV cen- 
ters) 


>90% 


Measurement 


Bulk magnetization 


Fluorescence: 
'quantum jump' 
technique 


Optical 


Electric, optical 


SQUID 


Controls 


RF pulses 


Optical, MW, 
electrical 


Optical 


RF, electrical, op- 
tical pulses 


MW, voltages, 
currents 



Table 1.3: Comparison of main features for different available techniques in QIP 



1.3. NMRQIP 



1.3.1 NMR- A suitable candidate for QIP 

A small scale NMR system is an attractive candidate for QIP for several reasons. 

• The fast reorientation of nuclear spins in liquids makes them fairly well isolated 
from the environment and therefore provide a good source of qubits. Since they 
are weakly coupled with the environment, nuclear spins have long coherence 
times of the order of seconds. 

• Nuclear spins can easily be manipulated and controlled using RF pulses. It is 
fairly simple to construct quantum logic gates using RF pulses and evolution of 
couplings. 

• The development of NMR over half a century itself makes huge difference in 
developing and optimizing so many experiments. NMR is a perfect experimental 
tool for a quantum mechanical theorist. 

• Modern NMR spectrometers are well developed. Though they are extremely so- 
phisticated, they can be easily controlled. 

1.3.2 NMR Qubits 

A spin-1/2 nuclei (^H, '^C, '^F, ''^P) in an external magnetic field will be placed either 
'parallel' or 'antiparallel' to the applied magnetic field. These two orthogonal states can 
be labeled as |0) and |1) state of a qubit. A spin-1/2 nuclei is most preferable since it 
is naturally equivalent to a qubit. A multi qubit system should have the individual ad- 
dressing capability and strong enough coupling constant with the farthest qubit. Since 
NMR is an ensemble system, it can address qubits separately depending on the slight 
variance in Larmor frequencies. Coupling is provided by J-interaction (in liquids) or 
dipole-dipole interaction or both (in liquid crystals, solids). A stronger coupling con- 
stant is always welcome since it reduces the time duration of gates operation. Nuclei 
with spin> 1/2 has both advantages and disadvantages. Each nucleus can be treated 
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as multiple qubits. Quadrupolar nuclei having spin / = (2^ - l)/2 (for N > 2) will 
have 2^ states. Thus, a spin 3/2 nucleus has 4 states and can be treated as a 2-qubit 
system provided individual transitions are selectively addressable. This can be achieved 
by introducing first order quadrupolar coupling. Quadrupole systems partially oriented 
in liquid crystals can form an ideal multiqubit system. 

Another possibility is that, one can take advantage of higher number of 'base' (A'^ 
2) into computation work. For example, a spin-1 system has three orthogonal states, 
which makes a 'qutrit' system [|20l . 

1.3.3 Initialization of NMR Qubits 

NMR is an ensemble system which deals with a large number (10'^) of identical spin- 
systems. At room temperature the nuclear energy levels are overwhelmed by the Boltz- 
mann energy distribution. In order to achieve a state like |000...0), all the spins should 
be brought to the ground state which needs extremely low temperatures or extremely 
high magnetic fields. It is almost a impossible task to prepare a 'pure' initial state at 
room temperature NMR. On the other hand for high resolution NMR, we need to have 
the systems to be in solution state. It is however pointed out by Cory et al [|45]| and 
Chuang et al [46] that the problem of preparing a pure initial state can be alleviated by 
preparing a pseudopure state which is a specially 'mixed' state mimicking a pure state. 
We will have a thorough discussion on preparing pseudopure states by various methods 
in Chapter 3. 

1.3.4 NMR Quantum Gates 

Single qubit gates can be thought of as a rotation in a Bloch sphere and can be imple- 
mented by the simple RF pulses |l5Tl[52l|53l. In NMR, it is convenient to understand 
the effect of RF pulses in rotating frame. 
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The NOT gate as given in expression (11.441) can be realized by a :7r pulse. 



hi 



-inly 



-I 



^0 

1 



NOT- 



(1.58) 



The factor -i can be ignored since it produces an undetectable global phase [|20l . Here 
the n rotation is achieved by using a RF pulse of power ojp for duration r and phase x 
such that ojpT = n. Similarly a Hadamard gate (see expression 1 1.471) can be realized by 
applying two RF pulses. 



1 1 
1 -1 



= U 



H- 



(1.59) 



A phase shift gate (Eq. 11.501) is equivalent to rotation by an angle (p about z-axis and can 
be realized as follows. 



= [(7r/2)_, .0, .(7r/2),] 



( \ 

1 



= f/p(0). 



(1.60) 
(1.61) 



Construction of multiqubit gates (e.g. CNOT) is achieved mainly by the proper exploita- 
tion of evolution of coupling and qubit specific RF pulses. The Hamiltonian for a pair 
of weakly coupled system in an isotropic medium is given by 



(1.62) 



where J is the coupling constant and (x>\ and tt>2 are Larmor frequencies of two spins. 
Now these J-coupling constant and Larmor frequencies are time independent fixed quan- 
tities and can not be turned off. But little tricks with the refocusing scheme can make it 
possible to overcome the problem. Historically most of the NMR experiments rely on 
this same refocusing technique. Consider at time ? = 0, the density matrix of the system 
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is given by + then The density matrix after time ? = r is given by 

+ ^ cos(ct>iT) + l] sin(a;iT) + l]. cos(co»2t) + ll ?,m{<x)2T) (1.63) 



X X 

InJllll 
> 



ll cos(a)iT) cos nJT + 2Iyll sin(a»iT) sinTr/r 

+/y sin(a»iT) cos nJr - 2lll^ sin(a)iT) sin nJr 

+1^ cos(a»2T) COS nJr + 21} I^, sin(a»2T) sin nJr 

+Iy sinCcuiT) cos nJr - llllj sin(a)2T) sin nJr (1-64) 



This generalized calculations can be simplified for many practical situations. For ex- 
ample, a n pulse on spin 2 at the middle of r period refocuses the J-coupling as well as 
chemical shift evolution of the spin 2 (see Fig ll.3h ). 

ll + ' ll cos(coit) + ll sm(t02T) + I]. (1 .65) 

Similarly, a n pulse applied on both the spins in the middle of r period refocuses the 
chemical shift while retaining the J-coupling evolution (see Fig. 11.3b ) 

^l + ' ^l cos{7tJt) + ll]!] sm{nJT) + 1] cosinJr) + 2l]ll sin(;r7T). (1 .66) 

Both J-coupling and chemical shift can be refocused over a time r by the pulse program 
shown in Fig. 11.3b . 

Il+Il 'Il+Il (1-67) 

A general method for refocusing has been described by Linden et. al. Il54l[55]| . A CNOT 
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2 



2 




(a) 



(b) 



(c) 



(d) 



Figure 1.3: Switching ON and OFF of interactions for the first qubit (Ii) over a time 
T. (a) Chemical shift ON, J-coupling OFF, (b) chemical shift OFF, J-coupling ON, (c) 
both chemical shift and J-coupling OFF (d) CNOT gate (t = 1/2 J). Filled rectangles 
indicates n/2 pulses and rest are representing a k pulses with phase mentioning on each 
pulse. 



gate can be achieved by the pulse sequence shown in Fig. 11.3d . 
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(1.68) 



where t = jj. Various gates and corresponding unitary operators has been studied in 



this context lEIl 152 EH. 



1.3.5 Numerically optimized quantum gates 

Apart from the above described methods of preparing quantum gates using ordinary 
'hard' RF pulses, there are more robust ways to synthesize any desired gate with high 
fidelity. These kind of pulses are often called 'strongly modulated pulses' (SMPs). 
SMPs are made up with suitable sequences of RF pulses whose amplitude, frequency, 
and phase are made time-dependent in such a way that it produces the best sequence 
of pulses with maximum robustness against RF inhomogeneity Il57l[58l [591 . There are 
few techniques available to find suitable SMPs for a given target operation [|57l[58l[60ll . 
Most of the techniques rely on the numerical optimization of the overall transformation 
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by searching the available parameter space. Fortunato et al used the stochastic search 
methods to construct SMPs [l57ll . while Khaneja et al described the gradient ascent pulse 
engineering (GRAPE) technique ll60ll . Designing an SMP for a given target operator re- 
duces the numerical search problem to set the control parameters that maximizes the 
fidelity ||59l . In this thesis we have used SMPs in many cases in later chapters. Our 
SMPs are synthesized by the stochastic search methods powered by genetic search al- 
gorithm ItSSll . 

1.3.6 Measurement 

NMR being an ensemble technique, the measurement of a quantity (say D) is done by 
measuring its expectation value IfTSll . 

s{t) oc ^/(0<<a'(0I^I'A'(0) 

k 

= YjP^^^^Yj'^'r^^y]{t){r\D\s) 

k r,s 

= tr{Dp{t)}, (1.69) 

where p is the density matrix and D is the detection operator. The free evolution of 
spin system under Zeeman and coupling Hamiltonian is detected over a time scale. The 
detection period is normally decided by the relaxation rate and it is recorded in the time 
domain. By Fourier transform, the NMR spectra can be transformed into frequency 
domain. 

The quantum algorithms are designed such that the final output state is an eigenstate 
in the computational basis [|20ll . In NMR, the computational basis is generally same as 
Zeeman basis (product basis). The eigenstates of product basis (diagonal elements of 
density matrix) correspond to the population distribution of a pseudopure state. For a 
two spin system, the general population distribution can be written as, 

p = cJl + C2II + c^llUl, (1.70) 
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where ci,C2,C3 are real constants. However, the traceless density matrix p character- 
ized by {ci,C2,Ci,} does not give rise to any signal, since it corresponds to longitudinal 
magnetization. Applying a n/l pulse results, 

p' = cJl + C2II + C32/]/2. (1.71) 

From the above equation one can get the values of c\ and ci but not C3 since the last 
term does not contain any single quantum coherences. The signal received for the above 
applied njl pulse will be, 

s{t) = tr[e~'^'p'e'^'(ll+f^] (1.72) 
= rr|.-'-^'(ci/l+C2/g.'-^'(/|+/^)) (1.73) 

where the C3 term won't produce any signal since 'multiple quantum' under such a 
evolution remains 'multiple quantum'. Here [l\ + represents the detection operator 
and "K is the evolution Hamiltonian as given in Eq. 11.621 

The complete procedure of characterizing a quantum state is through density matrix 
tomography. By this method all the coherence orders can be measured by a series of 
experiments. A detailed analysis of density matrix tomography is given in Appendix A 
and B and its explicit application in practical situations is given in later part of Chapter 
2. 

1.3.7 Coherence order 

Decoherence time in NMR quantum computers is generally related to the spin-spin 
relaxation time (^2), although this is a simplification PTi . T2 generally represents the 
decoherence time of a single spin coherences and the decoherence times for a multiqubit 
systems can be quite different. However T2 gives a rough estimation of the decoherence 
time of the system and in most of the cases it is of the order of multiple seconds in a 
liquid state system. The decoherence time in NMR is one of the best among all the 
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available QC techniques. The coherence time can be prolonged by applying specific 
dynamical decoupling sequences. A thorough discussions given in Chapter 4. 



There are multiple issues which posses serious challenge to NMR-QIP tasks. 

1 . Scalability of liquid state NMR is a big challenge that we need to be overcome. 
Indeed, going beyond 10 qubits is a very difficult task. High resolution liquid state 
NMR-QIP relies on weakly coupled systems. For a large order qubit system (say 
10), the J-coupling constants between two farthest spins are very small. Lower 
coupling constant means lack of well resolved spectra and more time consuming 
'quantum gates'. This problem can partly be resolved by taking partially orient- 
ing spin systems in liquid crystal solvents and thus introducing dipolar couplings 
along with J-couplings. Solid state NMR qubits has the potential for becoming 
the scalable qubits. However, at the moment solid state NMR system produces 
complicated spectra and allows lesser controlling technique. 

2. Lack of creating a 'pure' state is another limitation in NMR-QIP. In liquid state 
NMR, spins are in a highly mixed initial state and preparation of a 'pure' state 
needs extreme experimental conditions. However, one can prepare a 'pseudopure' 
state (PPS) which mimics a 'pure' state and effectively able to perform as a 'pure' 
state. PPS can be written as. 



where e denoting the purity factor and at high temperature limit it can be written 
as, 



1.3.8 Limitations of NMR-QIP 




(1.74) 



Nhv 



(1.75) 



6 
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e is of the order of 10"^ under normal conditions. Hence it can be seen that 
the amount of magnetization (signal) decreases exponentially with the number 
of qubits. However for 15-20 qubits, the magnetization may not be enough to 
carry out QIP tasks. Many ideas come forward to tackle this issue. Carrying 
out quantum computation in mixed state is one of them ||6T]| . 'One clean qubit' 
protocol needs only one 'pure' qubit having rest of them in 'mixed' state Il62]| . 

3. One of the major requirement to perform some algorithm is to create an entangled 
state. Many believe that the power of quantum computers is largely depend on its 
entangling phenomena. But it is proved by Braunstein et. al. that, small scale 
liquid state system at normal conditions lacks any kind of 'entanglement'. But 
still NMR is the only technique which implements the Shor's factoring algorithm 
till date and this algorithm needs an entangling state. These two contradictory 
aspects, led many people to think that 'entanglement' might not the necessary 
condition to be fulfilled to perform QIP tasks. A new measure of calculating 
non-classical correlations known as 'discord' is introduced recently ll63l[64l[65l . 

4. Crowding of frequency space is another challenge as the number of qubits in- 
creased. A weakly coupled spin system of N spin-1/2 nuclei gives rise to Nl'^^^ 
resonance lines. Thus, a 10 qubit system would have 10 sets of 512 lines. Resolv- 
ing spectra is a daunting task for this kind of situation. This can partly be resolved 
by synthesizing special molecules and using sophisticated spectrometers. 

5. Decoherence can be a potential problem for large scale qubit and also for solid 
state qubits. However, the relevant parameter is not the decoherence time itself 
but the ratio of decoherence time to the gate-time. 
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1.4 NMR - An ideal platform for studying quantum me- 
chanical phenomena 

Although 'nature prefers quantum mechanics' and it is the most fundamental physics, 
capable of explaining almost everything starting from photosynthesis effect to blackhole 
formation, it is rather a challenging job to 'feel' it in our daily life Il66l [671 [681 [691 
[Toll . The lack of experimental proof of quantum mechanical phenomena in the early 
days of its introduction led many eminent scientists commenting on its 'utility' and 
'completeness' [[7n[72ll . Due to its very nature, even today it is highly challenging to 
perform quantum mechanical experiments in laboratory [[7211 . 

The nuclear spins in an NMR provides an excellent test bed for performing various 
kind of quantum mechanical experiments in a highly controlled way [[73ll . The principles 
of an NMR can only be fully understood by quantum theory [HI. In return, NMR can 
be used as a prototype quantum mechanical testbed. The outcome of the quantum me- 
chanical probabilistic calculations mostly produces counter intuitive results which are 
difficult to 'digest'. Nonetheless, NMR has proven to be one of the leading architecture 
in performing various quantum mechanical phenomena experimentally [[741 [7511 . 

In this thesis we have shown some important experimental implementations of quan- 
tum mechanical phenomena which earlier thought to be intractable. There are various 
examples where a quantum mechanical phenomenon does not have a classical analogue 
and in this kind of situation it is rather difficult to 'understand' it |[76ll . For example, 
quantum contextuality is a kind of quantum mechanical phenomenon which has been 
proved by various quantum platform including NMR recently [[771 . The experimental 
results clearly shows quantum mechanically expected values which are counter intuitive 
to our macrorealistic world [[771 . 

Quantum objects behave differently than a macrorealistic object and there are certain 
inequalities (e.g. Bell's inequality) which can only be violated by quantum objects [TTSl 
[791 . To prove this kind of violation one needs to have an excellent quantum platform. 
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In Chapter 5, we have shown the violation of Leggett-Garg inequality for nuclear spins 
as predicted by quantum theory iTTSll . 

Bohr's complementary principle is another famous description regarding quantum 
mechanical objects. A consequence of the complementary principle is that one can 
not observe both 'wave' and 'particle' nature of a quantum object simultaneously [[80ll . 
However, recently it was proposed that, by using certain special experimental setups, 
one can simultaneously observe wave and particle nature of a quantum system IfSTTl . 
This requires a reinterpretation of Bohr's complementary principle. In Chapter 6, we 
have shown a detailed experimental study of this new experimental proposal and our 
results clearly suggest that there is indeed a necessity of revisiting Bohr's principle [|74l . 

Many new fields related to experimental quantum mechanics are coming up due to 
the fact that now we have some excellent quantum platform which are capable of car- 
rying out experimental work in a highly controlled way. Quantum chemistry and quan- 
tum biology are such two emerging fields which are making lot of progresses Il69l[82ll . 
Understanding all these phenomena experimentally is vital in pursuit of understanding 
quantum mechanics and its practical applications at large. 
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Chapter 2 

Density Matrix Tomography of Long 
Lived Singlet States 

The lifetime of nuclear singlet states can be much longer than any other non-equilibrium 
states under suitable conditions. In section 2.2, we introduced long-lived singlet (LLS) 
states and it's preparation by standard methods. In section 2.3, we introduced a robust 
density matrix tomography scheme which is particularly suited to study homonuclear 
spin systems with small chemical shift differences. In section 2.4, we have applied 
the tomography scheme to characterize the singlet states under various experimental 
conditions, revealing interesting features of LLS. This chapter ends with a conclusion 
given in section 2.5. 

2.1 Introduction 

The long lifetimes of nuclear spin coherence enables NMR spectroscopists to carry 
out a variety of spin choreography |l2l [181. Nuclear spin coherences decay over time 
mainly due to spin-spin relaxation and magnetic field inhomogeneity. Often, coherences 
are converted into longitudinal nuclear spin orders to study slow dynamical processes. 
But even the longitudinal spin orders decay toward equilibrium state due to spin-lattice 
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relaxation. Hence for a typical NMR experiment consisting of preparing and measuring 
certain correlated spin states, the ultimate time barrier was assumed to be defined by the 
spin-lattice relaxation time constant T\ [[H. 

It has recently been demonstrated that there exist certain 'long-lived states' which 
decay slower than the Ti values of individual spins Il83l [841 [8511861 [87l [8811891 190H. The 
long lived singlet states (LLS) has found wide range of applications ever since it was 
discovered by Carravetta, Johannes sen, and Levitt in 2004 [83]. Overcoming the Ti 
barrier has led to several exciting applications in studying slow molecular dynamics and 
transport processes ([911 [921 . precise measurements of NMR interactions |[93ll . and the 
transport and storage of hyperpolarized nuclear spin orders [[941 [951 [961 [971 [981 [991 . 

Bodenhausen and co-workers have demonstrated that the singlet spin-lock can also 
be achieved by RF modulations which are used in heteronuclear spin-decoupling IIIOOII . 
Detailed theoretical analysis of zero-field singlet states as well as singlet spin-lock have 
already been provided by Levitt and co-workers |[85l IIOIII and by Karthik et al |[86[| . 
Recently, long-lived states in multiple- spin systems are also being explored |[94lll02ll . 

2.2 Long-lived singlet states 

Let us begin with a simplest model consisting of a pair of spin- 1/2 nuclei. These two 
spins are labeled as and I^. The free-precession Hamiltonian of this system at labo- 
ratory frame can be written as, 

•H = ojJl + 0^2/? + InJnl^ ■ /^ (2.1) 

where coi and 0)2 are denoting the resonant frequency of the two spins respectively, 
whereas J\2 denotes the spin- spin coupling (J-coupling) between the two spins. 

The quantum states of the system can always be expressed with the combination of 
superposition of Zeeman states, namely jOO), jOl), jlO), |1 1). Here |0) denotes the angu- 
lar momentum of hjl along the magnetic field direction ('up' direction) and |1) denotes 
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the angular momentum of -fi/ 2 along the exact opposite direction of the magnetic field 
('down' direction). The four Zeeman product states together lead to one singlet state 
and three triplet states, 

101) -110) 
|7^+i) = 100), 

,^ ^ _ 101) + 110) ^^-^^ 

Mo/ -p 

V2 

|r-i) = |ll). 

Singlet states have many diff'erent properties compared to its triplet counterparts. Two 
most important properties are: 



(a) Singlet state is anti-symmetric with respect to spin-exchange, whereas triplet 
states are symmetric. 



(b) Singlet state has a zero total nuclear spin angular momentum quantum number 
[/^|5o) = 0], whereas triplet states have non-zero total nuclear spin angular momentum 
quantum number [I^\Tm) = 2\Tm)]- 



In the case of magnetically equivalent nuclear pair, the singlet state and the triplet 
states form an orthonormal eigenbasis of the internal Hamiltonian "Kj = 2nJI^ ■ I^. 
Singlet states can be prepared between two assymetric spins by imposing equivalence 
condition (either by lifting the sample out of Zeeman field or by applying suitable RF 
field acting as 'spin-lock'). But, being a zero-quantum coherence, singlet state itself is 
inaccessible to macroscopic observable directly. The traditional methods by Caravetta et 
al (831, described the way to access the singlet states by breaking its magnetic symmetry 
to convert into observable single quantum coherences. In this context we may note that, 
protons in Hydrogen molecule or in water is already in magnetic equivalence, but there 
is still no way to break the symmetry. 
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2.2.1 Why singlet state is long lived ? 



Any quantum state, deviating from its thermal equilibrium conditions, will return to 
its stable thermal equilibrium state through a mechanism known as relaxation. Hence 
it is needless to say that in NMR any observable quantum state is in non-equilibrium 
condition and that is the reason each state has its own lifetime. There are two major 
factors behind relaxation, (i) spin-lattice relaxation {Ti) and (ii) spin-spin relaxation 
(72). In majority of the cases T2 relaxation is much faster than T\. So the upper limit 
of the nuclear spin memory is bounded by the T\ irrespective of any experimental safe 
guard. However, there are some specialized cases where exceptions can be found, such 
as in the case of 'parahydrogen', where the spin state isomers lived much longer than 
T\ I1103II . Though the major reasons behind Ti, T2 relaxation depend on individual 
molecular property, other controllable parameters such as magnetic field inhomogeneity, 
temperature fluctuations, sample concentration etc. also contributes to the relaxation 
mechanism. 

Levitt and co-workers have successfully demonstrated |[83l [841 that the singlet state 
lifetime can be made many orders of magnitude longer than Ti for 'ordinary' molecules 
in solution state of homonuclear system. Now we will discuss some physics behind this 
astonishing long-lifetime of singlet states [11041 . The Hamiltonian for a pair of spins in 
magnetically equivalent environment is written as bellow: 



where, coq = jBq denotes the Larmor frequency of both (equivalent) the spins and Bq 
is the applied static magnetic field. The matrix representation of the Hamiltonian in 




(2.3) 
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singlet-triplet basis can be written HI 0411 as follows: 
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From the earlier equation it is seen that at zero field (a»o = 0), the triplet states are 
degenerate with same energy eigenvalues (^7r7). The energy diff"erence between the 
singlet and the triplet sates is InJ which is independent of the field. Since the Hamil- 
tonian is diagonal, there will not be any mix-up of singlet state population with triplet 
states' populations II104L However, triplet states among themselves equilibrate quickly. 
Eventually there will be a singlet-triplet transition resulting in the re-establishment of 
thermal equilibrium much slower than T\ relaxation time scale. The time constant for 
singlet-triplet equilibration is loosely termed as 'singlet lifetime' (Js) [I104L 

We already know that singlet states are 'antisymmetric' with respect to spin ex- 
change, whereas triplet states are 'symmetric' with respect to the spin exchanges. The 
major relaxation processes, including intra-molecular dipolar relaxation mechanism, are 
'symmetry preserving' in nature. Hence these relaxation mechanisms will not aff"ect the 
singlet-triplet conversion which requires symmetry transformations. These conditions 
make singlet states as a 'special' state which is immune to intra-molecular dipolar re- 
laxation, though it is the major reason behind T\ relaxation I1104L 

Previous discussion shows how necessary it is to get a magnetically equivalent pair 
of nuclear spins to realize the LLS. We need to create such a magnetically equivalent 
condition to create and persist in singlet states, but to get signal out of singlet states we 
need to break the symmetry. In the next paragraphs we will discuss about the techniques 
for magnetically inequivalent pair of nuclear spins. The Hamiltonian for a pair of 
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chemically inequivalent nuclei in present of Zeeman field can be written as follows: 

= ojxl] + 0J2II + InJI^ ■ 

= aJo(l + 5i)ll + a;o(l + ^2)/^ + InJI^ ■ /^ (2.5) 

where and 62 are the two chemical shifts of the two spins. The matrix representation 
of this Hamiltonian in singlet-triplet basis can be expressed as II104II : 
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(2.6) 



where, 

Yj5 = 81+62, ^6 = 6^- 62. (2.7) 

In this case, the matrix is not a diagonal matrix, hence the singlet and triplet states 
are not the eigenstates of this Hamiltonian. The off-diagonal term in the matrix (^(x»oA5) 
represents the possible conversion of singlet-triplet transition. This transition is directly 
dependent on the chemical shift difference between the two spins. Hence, even if we are 
able to prepare singlet states in an inequivalent pair of nuclei, it will not be long lived 
till it has some dependency on the chemical shift differences. Still, it gives us a clue 
to experience long-lived singlet states if somehow the chemical shift difference (A(5) is 
suppressed [I104II . In the next subsection we will discuss this method of chemical shift 
suppression in detail. 

2.2.2 Singlet Preparation in NMR 

So far we have learn that singlet states can not be observed for magnetically equivalent 
pair of spins, as it does not give any observable NMR signal, and even for the magneti- 
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cally inequivalent spin pairs because of the chemical shift difference barrier. 

The key to LLS revelation is to switch the magnetic equivalence 'on' and 'off' by 
some experimental manipulations I1104II . There are at least two well established proce- 
dures to do so - (i) field cycling and (ii) radiofrequency spin-locking [[83l[84l . By field 
cycling method, we can switch-off the magnetic field manually so that magnetic equiv- 
alence is established and then once again switch-on the magnetic field to convert into 
single quantum coherences. The other method (radiofrequency spin-locking) is more 
practical with least manual work. We will discuss this method in detail. 

Getting pure singlet states may be seen as a three step procedure: 

(i) Building singlet population. 

(ii) Applying spin-lock. 

(iii) Singlet detection. 

Building singlet population 

With the application of suitable rf pulses and delays it is possible to create a density 
matrix operator which represents a part of singlet states in it. The density matrix for a 
singlet population can be represented by the Cartesian product operator formalism as 
follows: 

|5o)<Sol = ^(|01)-|10))«01|-<10|) 

= \ (101X011 - 110X011 - 101X101 + 110X101) 
t (2.8) 

= -l(/|/! + /!/^)-7^/? + ii. 

The earlier equation shows that singlet populations can be constructed from zero 
quantum coherences and longitudinal magnetization of both the spins. Hence a little 
trick with the excitation of zero quantum coherences with appropriate phase can leads 
us to the singlet populations II104L The following pulse sequence is found [83] suitable 
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to create singlet state populations starting from equilibrium condition. 

90o - Ti - I8O0 - (ti + T2) - 9O90 - T3, (2.9) 



where T\ = I /4J, T2 = 1 /2Av and T3 = 1 /4Av. J and v are denoting the spin-spin 
coupling constant and chemical shift difference in Hz respectively. The 'offset' should 
be placed at the middle of the two spins for simplification. The above written pulse 
sequence works as follows : 

Initial 90o pulse brings the longitudinal magnetization to transverse plane. 

90o 

/I + f > -/I - f 

^ ^z "-y 

followed by the spin-echo with only J evolution for the duration of 1/2 J : 

r,-180o-r, 
- f > Iff + Iff 

ly ly ' -^^X^Z ^ ^^Z^X- 

During the subsequent T2 interval, there will be evolution under the isotropic chemical 
shifts. This delay (t2 = l/2Ay) is relatively shorter and can be ignored for any signifi- 
cant 7-evolution during this time. The product operator formalism goes as follows: 

2/1/2 + 2/1/2 _JL^ 2ff - Iff. 
Now a 90 degree y pulse will bring the density operator into zero quantum coherences. 

90,0 

2/;/2 - 2/1/2 , 2llJl - IfJl = -i(f,f - fJl). 

A further chemical shift evolution required for a phase corrected zero quantum coher- 
ence. 

2/1/2 _ 2/1/2 _2/l/2 _ 2/1/2 ^ _(/! /2 + /l/2)_ 

This may be rewritten as follows: 



-fjl - fjl = -|01)<10| - |10)<01| 

= |So)(Sol-|ro)<ro|. 
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Hence from the above calculations it is seen that the resulting density operator is in 
fact combination of the singlet state and one of the triplet states' population. Now our 
aim is to filter out the singlet state from the singlet-triplet population distribution. This 
can be done by radio frequency spin locking as discussed below in detail. 

Radio frequency spin-lock 

A spin-lock is a low power on-resonant continuous radio frequency pulse along the spin 
magnetization in transverse plane. This low frequency rf pulse keeps the magnetization 
from precessing in transverse plane. Hence this pulse can be used as a possible way 
to suppress the chemical shift differences. It is popularly known as a 'spin-lock' as it 
arrests the spin precession. 

The duration of the spin-lock may last for several minutes, triggering the possibility 
of severe probe damage. Hence one must be careful to select the rf spin-lock power and 
duration. There are two basic kinds of spin-lock, (i) Unmodulated spin-lock, and (ii) 
modulated spin-lock. 

(i) Unmodulated rf field is commonly known as 'continuous wave' (CW) irradia- 
tion. CW irradiation has constant amplitude and has no phase modulation over time. 
CW has shorter bandwidth and hence not useful for large chemical shift differences. 
Theoretically it is possible to apply more power for higher chemical shift difference 
systems, but that may cause serious damage to rf probes. 

(ii) Modulated lock can be realized by using CPD (composite pulse decoupling) 
pulses. As the name suggests, it is a phase modulated composite pulse, routinely used 
as a decoupling pulse sequence. In many cases it can outperform CW pulses as a spin- 
lock sequence. The bandwidth of CPD pulses are much larger compared to CW pulses 
and hence useful for larger chemical shift difference singlets. Commonly used CPD 
pulses are WALTZ- 16, GARP etc. 

During spin-lock the three triplet states' populations equilibrate rapidly under nor- 
mal relaxation procedure, whereas singlet population being itself immune to rf spin- 
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lock, decays much slowly. After the fast decay of triplet states, singlet state achieve 
its maximum purity (singlet correlation may reach upto 0.997). Eventually singlet state 
also decays despite rf spin-lock shielding, but with much slower rate than any other 
states. 

Now here we can recap the fact that singlet state itself is a zero quantum coherences 
and can not be directly accessible. Hence we must transfer the zero-quantum coherences 
to a observable single-quantum coherences to detect it. The following section describes 
the method in detail. 

Singlet detection 

The simplest method to detect singlet is to evolve it for a 1 /(4Av) chemical shift evolu- 
tion and followed by a strong 90° pulse. The transformation of density matrix operator 
are as follows: 

-^{Ilf-I\ll)-Illl + \1. 

This can also be written in terms of Cartesian product operator formalism: 

Now a simple 90() pulse brings the magnetization into observable single quantum 
coherence. 

1 1 90o 1 I 

-{lllll + I]!]) - I'jl + -1 > -(211 - Illl) - I'yly + (2.1 1) 

These shows the antiphase transverse magnetization for the spin pair. The character- 
istic spectra for this kind of antiphase magnetization shows a typical "up-down-down- 
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up" pattern in the NMR peaks. 

However one might notice that this way of detecting singlet states has less qualitative 
information. A better quantitative study can be carried out through tomographic method. 
In this context we have developed a robust density matrix tomographic technique which 
is particularly suitable for this problem. In the following section we will discuss the 
'density matrix tomography' scheme in detail. Later we will apply this tomography 
sequence on singlets for its characterization in various experimental conditions. 

2.3 Density Matrix Tomography 

The delicate nature of 'quantum states' makes it vulnerable to macroscopic-world. The 
inevitable last step for most of the quantum information processing and quantum simula- 
tion is the measurement of derived quantum states. In the case of an ensemble quantum 
system, the states are presented by density matrix. In order to measure these density ma- 
trices, many sophisticated schemes have been envisaged. 'Density matrix tomography' 
technique has proven its utility for mapping any quantum states with high accuracy. It 
enables us to measure all the elements of a general density matrix at any time point. The 
knowledge of the full density matrix of any quantum state is important for many reason 
e.g. (i) one can find the error in the experiment, since we already have the knowledge 
about theoretical density matrix, (ii) Measuring density matrix at different time points 
of a dynamic quantum algorithm gives the pattern of population and coherence trans- 
fers. In the following section we have presented a robust tomographic scheme in the 
context of NMR [[T06]I . 

Earlier schemes of tomography were designed in the context of quantum information 
processing Ill07[|108l . They required spin- selective rotations and transition selective in- 
tegrations of spectra. In homonuclear spin systems, particularly in spin systems, 
it is hard to design high fidelity spin- selective rotations owing to the small differences 
in chemical shifts (on the other hand, the heteronuclear singlet state is predicted to be 
short-lived HlOlll ). These spin selective pulses generally tend to be long in duration, 
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still introduce significant errors. Integration of individual transitions is also problem- 
atic since the transitions, particularly those with mixed line shapes corresponding to a 
general density matrix, may severely overlap. Tomography based on two-dimensional 
NMR spectroscopy had also been proposed H 10911 . This is a general method in the sense 
only one 2D experiment is needed to be carried out irrespective of the size of the spin 
system. However, the 2D method is time consuming. Also since it relies on fitting the 
2D cross-sections (along the indirect dimension) to mixed Lorentzian, the accuracy is 
limited by the quality of the fit that is achieved. In the following we present a robust 
density matrix tomography for a homonuclear weakly coupled two spins- 1/2 system 
which needs only non-selective RF pulses and integrations over each spin instead of 
individual transitions II106L 

The general traceless deviation density matrix consists of 15 independent real num- 
bers: 

Po r3 -I- is3 ri + isi rj -I- iss 



pi rg -I- iS(, r2 + iS2 
P2 + is4, 



(2.12) 



Here real elements p^, are populations and the complex elements r^. -I- is^ correspond to 
single (fc = 1 to 4), double (fc = 5), and zero {k = 6) quantum coherences. The elements 
below the diagonal are determined by the Hermitian condition pjt = plj. Since only 
single quantum coherences are directly observable, four combinations Ri := (ri -I- r2), 
S 1 := (si + S2), R2 '■= (r^, + r4), and ^2 := (s^ + S4) can be obtained from the integration 
of complex line shapes of spins 1 and 2 respectively. Now consider an RF sequence 
with propagator U, that transforms the original density matrix p into p' = UpUK Single 
quantum coherences of p' will lead to different linear combinations of various elements 
in p. Thus, by applying different propagators on p, we can measure the values of differ- 
ent linear combinations of various elements of p. The real and imaginary values of the 
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integration of jth spin in fcth experiment will be labeled as i?^ and S ^ respectively. Fol- 
lowing six one-dimensional NMR experiments were found to be sufficient to tomograph 
a two-spin density matrix: 

1. 1 

2. 90, 

4. 45,^. 180,- 

5. 45,^. 180,.^ 

6- 2i;-45.v^-180,--5|7 

Here 1 is the identity i.e., direct observation without applying any extra pulses. Av and 
J are the chemical shift difference and the scalar coupling respectively (both in Hz). 
The offset is assumed to be at the center of the two doublets and the RF amplitudes are 
assumed to be much stronger than Av. By calculating the propagator for each of these 
experiments, 24 linear equations are achieved. Solving these equations, gives the all 
unknown parameters of the density matrix. A detailed analysis of the density matrix 
tomography scheme is given in Appendix A. 

Now we will use this tomographic scheme on long-lived singlet states [11061 . 

2.4 Singlet State Characterization 

2.4.1 Observing through antiphase magnetization 

The singlet state was prepared by the RF spin-lock method and converted into antiphase 
magnetizations as described by Carravetta and Levitt [[84l using the pulse sequence 
shown in Figure 12. 2h . The RF spin-lock was achieved by either CW irradiation or by 
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100 50 -50 -100 Hz 



Figure 2.1: Part of the H spectrum of 5-bromothiophene-2-carbaldehyde (inset), 
displaying the doublets corresponding to the two spins used to study the singlet 
state. The sample was dissolved in dimethyl sulphoxide-D6 and all the experiments 
are carried out at 300 K. The difference Av in chemical shifts is 192.04 Hz and the 
scalar coupling J is 4.02 Hz. Scalar coupling to aldehyde proton was too weak to be 
observed. The spin lattice relaxation time constants (TO for the two spins obtained 
from inversion recovery experiment are 5.2 s and 6.2 s respectively for the spins 1 and 
2. 

WALTZ- 16 modulations. The RF offset was set to the center of the two chemical shifts 
in these experiments. The total magnitude of the antiphase magnetizations decays at 
different rates depending on the spin-lock conditions (Figures [231 and [Z4l) . The decay 
constants with CW spin-lock are 16.6 s (Figure [23k ) and 13.4 s (Figure [23h ) respec- 
tively at RF amplitudes of 2 kHz and 500 Hz. Under WALTZ- 16 spin-lock, the decay 
constants are slightly smaller, 16.2 s (Figure [24k ) and 12.8 s (Figure [24h ) respectively 
at 2 kHz and 500 Hz. Nevertheless, these values are about 2 to 3 times the Ti values of 
the individual spins implying the preparation of long-lived singlet state. 

In this scheme, the integrated magnitude spectrum is usually monitored as a func- 
tion of spin-lock time. The contributions from the spurious coherences may not be 
eliminated in this process. Further, the double quantum coherences, if any, are not ob- 
served at all. Our interest is to quantify the singlet content in the instantaneous state p{t) 
during the spin-lock. One might guess that the singlet content is maximum in the begin- 
ning and exponentially decays with the spin-lock time. Further, one may also guess that 
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90 180 

(b) " ° 



90 



90 



^3* 



Tomography 



Figure 2.2: The pulse sequences for the preparation of singlet states and detection 
via (a) converting to antiphase single quantum magnetization and (b) tomography of 
singlet states. Here t\ = 1/(47), T2 = 1/(47) + l/(2Ay), and T3 = l/(4Av), with Ay 
and 7 being the chemical shift difference (in Hz) and the scalar coupling respectively. 
T4 is the duration of spin-lock. 



CW spin-lock is superior to WALTZ- 16 spin-lock at all timescales. But the following 
tomography results provide a different picture [11061 . 



2.4.2 Tomography under varying spin-lock duration 



The pulse sequence for the tomography of singlet states is shown in Figure 12.2b . The 
density matrix ofthe singlet state is |5o)(5'ol = |l+p.v, with the traceless part p., = -/'-/^ 
being the product of spin angular momentum operators of spins 1 and 2. The correlation 
of the theoretical singlet state operator in the instantaneous experimental density 
matrix p(t) (obtained from tomography), 

VA trace [p(t) ■ p,] 

(ps)(t) = (2.13) 
ytrace [p(0 J ■ trace [p;J 
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gives a measure of singlet content in p{t). The normalization used in the above ex- 
pression disregards the attenuation of p{t) itself. Similar definitions can be applied to 
calculate the correlations (/'), (|ro)(ro|), etc. We monitored the correlations as a func- 
tion of spin-lock time T4 from s to 30 s in steps of 0.5 s under different spin-lock 
conditions using the sequence shown in Figure [X2b . The results are shown in Figures 



l2.3l and [Z4l 3D bar plots of full density matrices at two particular spin-lock conditions 

are shown in Figure [231 

The Figures [23b . 12. 3i . 12.4b . and 12.41 indicate correlation (py) as a function of spin- 
lock time under various spin-lock conditions. In all the cases, the initial correlation is 
about 0.8. This is expected, since the initial state prepared by the pulse sequences in 
Figure [23[ just before the spin-lock is actually 



With CW spin-lock at a high RF amplitude of 2 kHz (Figure 12.3b - 12.3^ ). the singlet 
correlation (p,) quickly reaches to 0.95 in 0.5 s of spin-lock time (Figure [23b ). Most 
of the spurious coherences and the residual longitudinal magnetizations created during 
the preparation are destroyed by the RF inhomogeneity during spin-lock. Figures 12.3^ 
and l2.3h reveal that the initial correlation (|ro)(ro|)(0) is -0.7 ~ 1/ V2 which is just 
expected . Within 0.5 s, the \T()){T()\ content is rapidly reduced. But complete equili- 
bration of triplet levels takes about 5 s. Interestingly, there is a sudden build-up and 
gradual fall of double quantum coherence as seen in Figures [23f and Figures [23h i. As 
the singlet state gets purified, (pv) exceeds 0.99 in 6 seconds and reaches a maximum 
value of 0.994 at 9.5 s. After about 18 s, (p,) starts decaying below 0.99, probably due 
to the gradual conversion of singlet state to other magnetization modes via the triplet 
states by relaxation mechanisms. On the other hand, there is a gradual build up of y- 
and z-magnetizations (Figures 12.3d and 12.3b ) in a similar way as that of a steady state 
experiment [[Tl llllL Nevertheless, the singlet correlation remained above 0.95 till 30 
s. The x-magnetization and the double quantum coherence (Figures [23b and l2.3f ) re- 




p(0) = |5o)<5ol-|ro)<ro|. 



(2.14) 
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Figure 2.3: Data characterizing the singlet state under CW spin-lock at an RF ampli- 
tude of 2 kHz (a to g) and of 500 Hz (h to n). The spin-lock duration T4 was varied 
from s to 30 s in steps of 0.5 s in each case. Dots in (a) and (h) correspond to the 
total magnitude of antiphase magnetization obtained from the pulse sequence in Figure 
12.2b . Singlet decay constant Ts was obtained by using an exponential fit (smooth lines 
in (a) and (h)). During each fit, first two data points were omitted in view of strong 
spurious coherences created by the imperfections in the pulses. Remaining graphs are 
the results obtained from tomography using the pulse sequence shown in Figure 12.2b . 
They correspond to the correlations: (p.,) (b and i), (7^) (c and j), {ly) (d and k), (/f) 
(e and 1), (/|/+ + /l/^) (f and m), and Cq - {\Tq){Tq\) (g and n), with spin numbers 
p - {1,2} and triplet subscripts q = {"I'O^'^l). 
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Figure 2.4: Similar data as in Figure 1231 but under WALTZ-16 spin-lock at an RF 
amplitude of 2 kHz (a to g) and of 500 Hz (h to n). The graphs correspond to total 
magnitude of antiphase magnetization (a and h), (ps) (b and i), (/^ > (c and j), (ly) (d 
and k), </f > (e and 1), (/j/^ + (f and m), and c,^ - <|r,^>(r^|> (g and n), with' spin 
numbers p = {1,2} and triplet subscripts q = {-1,0, +1). 
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Figure 2.5: Bar plots showing (a) traceless part of the theoretical singlet state 
density matrix, (b) experimental state after 15 s of WALTZ- 16 spin-lock at an RF am- 
plitude of 2 kHz, and (c) experimental state after 14 s of WALTZ- 16 spin-lock at an 
RF amplitude of 500 Hz. The upper and lower traces correspond to the real and imag- 
inary parts respectively. The singlet conflations in (b) and (c) are respectively 0.997 
and 0.547. The density matrix in (b) shows significant decay, but still has high singlet 
content! The real part of the density matrix in (c) shows significant double quantum 
artifact. 

mained small during the period of high correlation. After the initial differences, the 
triplet states equilibrate in about 6s, and remain steady then onwards (Figures [23^ and 

With CW spin-lock at 500 Hz , the singlet correlation reaches only up to 0.94 again 
at about 9 s and then steadily drops to 0.71 at 30 s (Figure [231 ). The increased buildup 
of X-, y-, and z- magnetizations with the reduction of the spin-lock power can also be 
noticed (Figure [23] - 12.31 ). 

Under WALTZ- 16 spin-lock (Figure [Z4l) . all the graphs are characterized by oscil- 
lations that are either in-phase or anti-phase. The origin of oscillations probably lies in 
the cyclic nature of WALTZ- 16 modulation. 
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At an RF amplitude of 2 kHz, the maximum singlet correlation of 0.997 was reached 
at 15 s(Figure [24b ). The 3D bar plot of the density matrix corresponding to this case 
is shown in Figure 12.5b . More interestingly, (p.,) peaks seem to maintain above 0.99 
till T4 = 28.5 s, i.e., about 10 s longer than the CW case! Thus, for certain values of 
spin-lock durations, WALTZ- 16 provides purer singlet states than that of CW. 

The singlet correlation under WALTZ- 16 spin-lock at 500 Hz displays stronger os- 
cillations (Figure [241 ). Despite the oscillations, the singlet correlation reaches as high 
as 0.96 at 13.5 s. Again it can be noticed that good singlet content is held for longer 
periods by WALTZ- 16 than the CW of same amplitude. For example at 500 Hz RF 
amplitude, WALTZ- 16 gives a singlet correlation of 0.94 at T4 = 27 s, while that for 
CW it is only 0.79. 

2.4.3 Offset dependence 

Theoretical and numerical investigations on the offset dependence of singlet spin-lock 
has been have been carried out by Karthik and Bodenhausen If86]| and by Pileio and 
Levitt HIOIL Robustness of various modulation schemes with regard to offset of sin- 
glet spin-lock have been demonstrated by Bodenhausen and co-workers HIOOL Here we 
probe the offset dependence of singlet evolution using tomography [I106II . Figure IZ6l 
shows the experimental data obtained from a series of singlet state tomography exper- 
iments, each time varying the RF offset of the spin-lock. The RF offset was measured 
from the center of the two chemical shifts. Again the experiments were carried out 
under the following spin-lock conditions: 



(i) 



CW for 15s (Figure llSi-IIlf), 



(ii) 



CW for 28.5s (Figure [16^-1161), 



(iii) 



WALTZ-16 for 15s (Figure I23fai-l231-). and 



(iv) 



WALTZ-16 for 28.5s (Figure [IB-IIIK). 
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Figure 2.6: Correlations calculated using the density matrix tomography of singlet 
states prepared with different spin-lock conditions: (i) CW spin-lock at 2 kHz for 15s 
(a to f), (ii) CW spin-lock at 2 kHz for 28.5s (g to 1), (iii) WALTZ-16 spin-lock at 2 
kHz for 15s (m to r), and, (iv) WALTZ-16 spin-lock at 2 kHz for 28.5s (s to x). In each 
case, the horizontal axis indicates the RF offset Ay during the spin-lock. The offset is 
measured from the center of the two chemical shifts. The rows correspond to : (p.,) 
(b and i), (/^> (c and j), {I^.) (d and k), (/f > (e and 1), (/j/^ + l[f) (f and m), and 
= {\Tq){Tq\) (g and n), where spin numbers p - [\ (pluses), 2 (x's)} and triplet 
subscripts <7 = {-1 (dots), (circles), -i-l (pluses) 
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Figure 2.7: Pulse sequence for the creation of long lived singlet states in a 3 spin 
system {AMX). (a) Anti-phase singlet magnetization to be accessed via spin-3, (b) 
qualitative measure of singlet correlation is done by state tomography. t\ and T2 are 
optimized delays in a way that both /13, J23, get a 7r/2 /- evolution, T3 - 4^;^, and T4 
is the spin-lock duration. 



The graphs indicate that the WALTZ- 16 scheme is far superior compared to CW in 
preserving the singlet correlation at high RF offsets. The singlet correlations with 2 kHz 
CW drops below 0.5 for an offset of 50 Hz. However, WALTZ- 16 at 2 kHz amplitude 
maintains a high correlation of 0.97 at 28.5s, even with an offset of 2. 1 kHz. In the case 
of CW spin-lock, rapid build up of y-magnetizations can be noticed with the increase of 
the RF offset [Hoi. 



2.5 Long lived singlet states in multi-spin systems 
2.5.1 Long lived singlet states in a 3-spin system 

In this subsection we will describe the methods for preparing long lived singlet states 
in a 3-spin system {AMX). We have extended the procedure of the 2-spin system as 
described in the previous sections. The singlet population distribution between any two 
spins can be prepared in presence of a third spin. The pulse sequence relies on the 
refocusing of the unnecessary couplings. The NMR pulse sequence is shown in Fig. 

60 



2.5. Long lived singlet states in multi-spin systems 




40 -40 -240 Hz 




5 10 15 20 25 30 35 



Spin-lock duration in Sec. 

Figure 2.8: Experimental results of 3-spin LLS. (a) Acrylonitrile dissolved in CDCI3, 
where 3 protons acting as a three spin homonuclear system, (b) The reference 
spectra of Acrylonitrile in a 500 MHz spectrometer, (c) and (d) showing the antiphase 
spectra of spin-3 after a spin lock duration of 5 s and 40 s respectively, (e) The solid 
curve showing the antiphase magnetization decay and dotted curve showing the singlet 
correlations obtained from tomography over a duration of spin-lock (T4) time. 
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Figure 2.9: Pulse sequence for the creation of long lived singlet states in a 4 spin 
system {AMXY). (a) Singlet states accessed by transferring it into anti-phase magneti- 
zation, (b) qualitative measure of singlet correlation is done by state tomography. t\, 
Tj are optimized delays and T4 is the spin-lock duration. 9Qph denoting a optimized 
phase njl non-selective pulse. 



12.71 In this particular example we have prepared the singlet population between spin-1 
and spin-2. Singlet population is accessed by transferring the magnetization into spin- 
3 . The quantitative measure of singlet magnetization is done by the pulse sequence 
shown in Fig. I2.7h . The extensive tomographic method of accessing singlet correlation 
has also been performed. The experimental results are shown in Fig. 12.81 The decay 
of antiphase magnetization and tomographic correlation is shown in Fig. 12.8b . The Ti 
time for all the three spins are roughly 6 sec, while the LLS time (Tus) is found to be 
17.9 sec. Hence the ratio Tus /Ti « 3. The nature of this plot is similar to the spin-2 
system and the reason for this is given in previous section. The 3-spin density matrix 
tomography scheme is described in detail in Appendix B. 



2.5.2 Long lived singlet states in a 4-spin system 

We have prepared two pair of singlet states in a 4 spin AMXY system. The exact pulse 
sequence is shown in Fig. 12.91 We were able to prepare simultaneous singlet states 
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Figure 2.10: Experimental results of 4-spin LLS. (a) Aspirin dissolved in CD3OD, 
where 4 protons acting as a four spin homonuclear system, (b) The chemical shifts and 
J-coupling constants in Hz are shown in a table format, (c) The '// reference spectra 
of Aspirin is at the bottom trace. The antiphase spectra are shown in upper trace after 
a spin lock duration of 3 s and 10 s respectively, (d) The solid curve showing the 
antiphase magnetization decay and dotted^Surve showing the 1-2 singlet correlations 
obtained from tomography over a duration of spin-lock (T4) time. 
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in between spin-1 and spin-2 and also in between spin-3 and spin-4. The J-evolution 
delays (ri and tj) are calculated in a optimized way. The traditional method of accessing 
singlet states is by converting it into single quantum coherences. The antiphase spectra 
of aspirin are shown in Fig. 12.10b . We have done density matrix tomography to calculate 
the correlation at various spin-lock duration. The nature of antiphase decay and singlet 
correlation profile matches with the previous cases (2-spin and 3-spin) and hence got the 
similar explanation. We have found Ti « 3s and Tus ~ 65, hence the ratio Tlls /Ti « 2. 
This also proves the long lived nature of the prepared singlet states. 

2.6 Conclusions 

Study of singlet state is important not only because of the interesting Physics that makes 
it long-lived, but also because of its potential for a number of applications. We have 
studied the singlet state directly and quantitatively using density matrix tomography. 
A new set of tomography sequences have been introduced for this purpose. The den- 
sity matrix tomography provides a tool not only for characterizing various spin-lock 
schemes but also for understanding the spin dynamics during the spin-lock period. 

The singlet state is preserved with CW spin-lock as well as with WALTZ- 16 spin- 
lock at two different RF amplitudes: 2 kHz and 500 Hz. The results indicate that at high 
RF amplitudes, both CW and WALTZ- 16 achieve high singlet content. An important 
feature of singlet state is that it gets purified by itself during the spin-lock, simply be- 
cause of its longer life time compared to the spurious coherences. There exist optimum 
spin-lock values at which the singlet correlations are maximum. While WALTZ- 16 
shows significant oscillations in the singlet purity, for certain intervals of spin-lock it 
gives better performance than CW and holds the singlet content for longer intervals of 
time. The dependence of correlations with the RF offset during the spin-lock are also 
studied under both CW and WALTZ- 16 schemes. It is found that WALTZ- 16 is far 
superior in preserving the singlet state at large RF offsets. 
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Chapter 3 



Preparation of Pseudopure States 



Using Long Lived Singlet States 



An introduction to pseudopure states (PPS) is given in section [STl In section IX2l a gen- 
eral review of various methods proposed for preparation of pseudopure states in NMR 



have been described. In section [331 a new method for creating pseudopure states from 
long-lived singlet states is presented. In this context, we have developed a robust, scal- 
able quantum circuit for creating PPS for any number of qubits. Experimental results 
are shown and discussed for multi-qubit systems in the later part of the section [331 



A 'pure state' is a quantum state which can be represented by a single ket A quantum 
state which can not be represented by a single ket format in any of the basis is termed as 
'mixed state' . In a mixed state system, we are only able to say that the system has certain 
probabilities (say pi, p2... ) of being in different states (say |1), |2)... ). For example, 
let us consider a two level system of energy difference AE of about 2eV (e.g. electronic 
levels of sodium). At room temperature of 300 K, the ratio between ground and excited 



3.1 Introduction 



3.1.1 A pure state and a mixed state 
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states population is of the order of 1 : e ' ~ 1 : 10" . For all practical purposes, 
we can assume that such a state exists in a pure ground state at thermal equilibrium. In 
another case, if we consider a much shorter energy gap difference of the order of IfieV 
(such as the case of nuclear spin's energy gap for hydrogen at lOr external magnetic 
field), then the populations of the two states are roughly in the order of 1 : 0.9999. The 
later case described a perfect mixed state condition. 



While a pure classical state is represented by a single moving point in phase space 
having definite spatial (cji, ...qf) and momentum coordinates (/?,, .../?/), a classical mixed 
state is described by a non-negative density function p(^,, ...qf, pi, ...pf), such that the 
probability that a system is found in the interval dqi...dpf at time t is pdqi...dpf. The 
'quantum analog' of the classical pure state is represented by a single state vector, while 
the quantum analog of a classical mixed state is represented by the density matrix II112II . 

An ensemble is in an idealized pure state if all of its' members are in state \ip{t)). 
The density operator of a pure state \ilj{t)) = Yui Ci{t)\i) is given by, 

Ppure = = E Z Ci<^t)c){t)\i){j\. (3.1) 

' i 

A pure state density matrix, P\n) = \n){n\ is a projection operator since, 

P\n)Y^Ci\i) = Cn\n), and ^P|„) = l. (3.2) 

1 n 

Since all the subsystems behave identically, one has complete knowledge of all parts of 
the system. Pure states are also known as the 'states of maximum information'. A pure 
state has two basic properties : 



plure=P^ and /'r(pJ„J = l. 
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The density matrix of an ensemble in a mixed state is obtained by ensemble average : 

i i k 

where pk denotes the probability of the system of being in state \pk with YikPk = 1 and 
the bar denotes the ensemble average. p„„v is not a projection operator and is also known 
as state with 'less than maximum information'. It holds the property of : 

pL^P' and tripl^J < I. (3.5) 

Another way to represent a wave function and hence a density matrix is through the 
'Bloch sphere' . A pure state is a geometrical point on the Bloch surface. For a single 
qubit, can be written as : 

= cos(e/2)|0) + e'^ sin(0/2)| 1 ), (3.6) 

with < 6 < n and < < 2n. Here |0) and |1) are two orthonormal basis vectors of 
the two-level system (single qubit). Any arbitrary density matrix p can be expanded in 
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terms of Pauli operators cr, : 



p = -(1 + 



(3.7) 



where f is called as Bloch vector, the radius vector of the state from the origin. The 



For a pure state tr{p^) = 1, i.e., |r| = 1 and for a mixed state |r| < 1. Hence a pure state 
represents a point on the surface of a Bloch sphere, whereas any point other than on the 
surface represents a mixed state. The origin of the Bloch sphere (|r| =0) represents the 
maximally mixed state with tr{p^) = ^. 

3.1.2 Necessity of Pure states in QIP 

In order to carrying out information processing, a quantum register must satisfy a set of 
criteria laid out by DiVincenzo I1124L An important criterion is the ability to precisely 
initialize the register to a desired ket of the computational basis. It has been shown that 
highly mixed state may not be used to create an entanglement II113II in QIP. Again, a 
mixed input state leads to a mixed output state which are generally difficult for analysis 
||30]| . In one particular case it can be show that a mixed state input gives wrong out- 
put results. Grover search algorithm for a two-qubit system starting from a pure state 
gives the desired answer with full probability. While starting from a mixed state, in- 
stead of getting a definite state, one gets a result with probabilities of different states. 
On the otherhand, we can find some algorithms which have no requirement of starting 
from pure initial state, e.g. Deutsch-Josza algorithm II114L Deutsch-Josza algorithm 
decides whether the function is 'balanced' or 'constant' in one iteration. If the initial 
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eigenvalues of p are ^(1 + |r|). As density operator must be positive, we get |^ < 1. The 
purity of the density matrix can be measured by squaring it. 




(3.8) 



=^(l + l^'). 



3.1. Introduction 




Figure 3.2: Simulations of Graver's search algorithm for a two qubit system starting 
from pure-state and mixed-state condition. The state being selectively flipped is |01>. 
(a) Pure input state |00)(00|, (b) output state 1 10)( 10| with full probability, (c) a complete 
random mixed state population input, (d) the answer of the Grover's search algorithm 
with no definite result. 

state is a pure state, then the in-phase output signal denotes it as a balanced and the 
anti-phase output signal denotes it as a constant function UllSII . In other case, when the 
initial state is a mixed state, the output signal consists of all possible transitions for a 
constant function, whereas for a balanced function, at least one of the transition will be 
missing [|50l I116[ I117[lll4l . Since the Deutsch-Josza algorithm is used only to distin- 
guish between a balanced function and a constant function, mixed initial state should 
be sufficient for this purpose. Again there are some algorithms that take advantage of 
the mixed nature of an ensemble systems HI 18[ II 191 11201 . These algorithms utilize the 
power of parallel processing of the identical quantum registers. In a broad overview, 
however, initialization of quantum registers remain an important requisite for a large 
scale quantum computer. 

3.1.3 Pure states in NMR 

The density operator of an ensemble system characterized by the Hamiltonian "K at 
thermal equilibrium is given by : 

P = (3.9) 

where k is the Boltzmann constant and T is the temperature of the system. Z represents 
the partition function and is given by Z = tr ^e''^^''^^. For a N-spin homonuclear system 
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at room temperature, the eigenvalues of are much smaller than kT value, hence the 
equilibrium density operator can be approximated to. 



where y is the gyromagnetic ratio, Bq is the strength of the Zeeman field, I{ is the 
longitudinal component of spin operator for spin, and p^ev is the trace-less 'deviation 
density matrix'. The coefficient e = HyBo/kT is normally very small; e ~ 10"^ for 
hydrogen nuclei in a Zeeman field of 10 T, at room temperature (300 K). Thus a normal 
NMR system represents a highly mixed state [fTllIl 112511 . However, only the deviation 
density matrix contributes to the NMR signal, which is smaller by a factor of e compared 
to the system which is in a pure state, i.e., if all the spins were in one state. It seems 
from the Eq. 13.101 that the only way of preparing a pure state in NMR is by using 
extremely high magnetic fields at extremely low temperatures. While achieving higher 
magnetic fields than a few tens of Tesla, is still a technological challenge, carrying out 
experiments at low temperatures is associated with many technical problems including 
among others, liquid to solid phase transition. However, in Eq. I3.10[ we should notice 
that we have taken only the nuclear spin part. Eff'orts have been made to bring-in other 
interactions which will enable preparing pure states without the requirement of using 
unrealistically high magnetic fields or undesirable low temperatures. One technique 
is Optical pumping, wherein one tries to transfer electron polarization into the nuclear 
spin system. This includes for example, a polarization transfer from a laser polarized 
noble gas atoms like ^^^Xe to molecule of interest. While this technique is widely used 
in magnetic resonance imaging (MRI), it is expected to be relatively inefficient with 
regard to preparing a pure state [I121II . 

The spin temperature can nevertheless be reduced by using parahydrogens HI 2611 or 
by using Dynamic Nuclear Polarization (DNP) [|127II . In future, either or both of these 




(3.10) 



1 

= ^(1 + ^Pdev), 
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techniques may be available for preparing NMR quantum registers into almost pure 
states II122L The existing approach for initializing NMR registers is however based on 
specially prepared mixed states known as pseudopure states (PPS). 

3.1.4 Pseudopure states 

In 1997 it was suggested independently by Cory et al [|45ll and Chuang et al |f46] that 
a specially prepared mixed state can be prepared, known as pseudopure state (PPS), 
that can simulate a pure state in NMR. PPS are isomorphic to pure states for several 
computational problems [|45ll46ll . From Eq. I3.10[ it is seen that the equilibrium density 
operator (p) can be split into two parts, the identity part (1) and the deviation part (pdev)- 
The identity part comes as a uniform background and does not give any kind of NMR 
signals. The NMR signal solely depends on the deviational part of the density operator 
(Pdev)- Although preparation of pure states in NMR is a very difficult and technically 
challenging task, preparing a pure deviational density operator is rather a easy one. 
Hence a pseudopure state (PPS) is an ensemble with a pure deviational density operator. 
It is also known as 'effective pure state' since it mimics a 'pure state'. Fig. 13.31 shows 
representative population distributions of a two-spin system in equilibrium (Fig. 13. 3h ). 
the deviation part distribution (Fig. 13.3b ). and the pure deviational distribution (Fig. 
13.3b ). Fig. 13.3b represents the 'pseudopure state'. Equation 13.101 can be rewritten as 
follows : 

Peg =^(.'^ + ^Pdev) 

\ (3.11) 

= — (1 - 6)1 + 6|(A)<<AI = Ppp., 

where e is a measure of the magnetization retained in the pseudopure state and it usually 
gets halved with every additional qubit [I128II . The unit background is invariant under 
the Hamiltonian evolution, does not lead to NMR signal and is ignored [|45l : 

UpppsU'' = ^(1 - 6)1 + eU\ifr){iJ^\U\ (3.12) 
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Figure 3.3: Representative population distribution of a two-spin system at room tem- 
perature when the system is in equilibrium (a), the deviation distribution (b), and the 
deviation distribution corresponding to |00) pseudopure state. 



Thus the equilibrium density matrix of a single spin- 1/2 nucleus is always in a pseu- 
dopure state. Initializing a multi-spin system into a pseudopure state however is es- 
sentially a non unitary process ||59l . So far, PPS remains as the leading technique to 
simulate a pure state in NMR QIP. 



3.2 Methods for preparing pseudopure state 

Several methods have earlier been proposed for the preparation of pseudopure states. 
These methods involved in averaging the magnetization modes over the sample space 
(called 'spatial averaging' I1129II ). or over spin space (called 'logical labeling' Il46l[50l '). 
or over several transients (called 'temporal averaging' [I130II ). In some cases subsystem 
pseudopure states are easier to prepare either by transition selective pulses II135II or by 
coherence selection using pulsed field gradients II128L but these methods invariably 
result in loss of a qubit for further computation. 

Implementation of a bulk quantum computation can be described in a general way 
as the transformation of an initial density matrix po into a final density matrix p„„, ac- 
cording to lfT30H . 

Pou, = J]r,CPjPoP]C'^R], (3.13) 

'J 

where Pj is the preparation operator, C is the unitary transformation corresponding to 
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Figure 3.4: Preparation of |00) pseudopure state by using temporal averaging method. 
Representative deviation population distribution at thermal equilibrium (a). Modified 
deviation populations after inverting populations between |01) , |11) (b) and |10), |11) 
(c). The population distribution of |00) PPS is obtained by temporal averaging tech- 
nique i.e. by adding (a), (b), and (c). 



the computation to be performed, and i?,- is the post-processing unitary operator, usually 
the identity. 

It can be noticed that for a single spin- 1/2 system in thermal equilibrium, only two 
energy levels are possible. Hence, with an excess population in the lower energy level 
than higher energy level, the system is in a pseudopure state by default. However, for 
higher order spin systems (more than one) one has to suitably manipulate the spin order 
to achieve the desired pseudopure state. In the following, we will describe a few of the 
most successful methods for preparing pseudopure states in NMR. In section 13.31 we 
will introduce a new method for preparing pseudopure states by exploiting long-lived 
nature of singlet states. 



3.2.1 Temporal averaging 

Temporal averaging method is one of the earlier techniques that was proposed by Knill 
et. al. [I130L This method is based on the randomization of the probabilities of other 
states by adding multiple experiments. By averaging out all this experiments, we can re- 
tain the desired component whereas other undesirable terms cancel out. The basic idea 
of this method is quite similar to phase-cycling technique, routinely used in NMR. For a 
N-qubit system, the exhaustive averaging involves cyclically permuting the non-ground 
states in 2^ - 1 different ways. This permutation achieved by CNOT gates which can be 
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implemented by using spin-selective pulses and evolution of couplings or by transition 
selective pulses. 

The method for preparing |00) PPS for a 2 qubit system is described here (Fig. 13.41) . 
Figure l34b shows the deviational part of the equilibrium population distribution. After 
permuting the populations of |01) and |11) by applying a transition selective n pulse, 
we obtain the population distribution similar to Figure [X4b . Similarly, a permutation of 
populations of |10) and |11) can also be obtained as shown in Figure [34b . Now adding 
these three population distributions, we can obtain the |00) PPS as shown in Figure [34b . 

Though this method is useful for smaller number of qubits, it becomes quite labo- 
rious for higher number of spin systems. The number of experiment (cyclic permu- 
tations) increases exponentially with the number of qubits. For example, preparing a 
4-qubit PPS, would require 15 cyclic permutations. Because of this limitations var- 
ious improvements on this method have been proposed such as by using non-cyclic 
permutations and applying unequal weightings to different permutations. Utilizing this 
improvements, one can find out the 4-qubit PPS by only 5 experiments 111311 . Choosing 
random permutations may turn out effective for very large spin systems II130L 

3.2.2 Logical labeling 

The logical labeling technique [[461 II 0711 can be used to prepare conditional (or subsys- 
tem) pseudopure states. It exploits the fact that in an N-spin homonuclear system , there 
are ^Ca^/i levels having equal population distributions. The basic idea of this approach is 
to find the subset having similar kind of pattern of a pseudopure state and then 'relabel' 
these states. 

Let us consider a case for a three-spin homonuclear system with the equilibrium 
deviation populations as shown in Figure [33b ). By permuting the states |001) and |010) 
with 1 101) and |110) respectively, we obtain the population distribution similar to the 
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(a) (b) 

Figure 3.5: Preparation of pseudopure state by using logical labeling technique in a 
three qubit system, (a) Representative deviation population at thermal equilibrium (b) 
The three qubit system can be thought of consisting two subsystems corresponding to 
and 1 state of the first qubit. After inverting the transitions 001 101 and 010 <-> 110, 
subsystem is in |00> pseudopure state and the 1 subsystem is in -|11> pseudopure state. 

Figure [33b . Now these states correspond to the deviation density matrix, 

p = |0)<0| ® (2|00)<00| + 27) + |1)<1| ® (3/ - 2|11)<1 1|). (3.14) 

From the above equation it is clearly visible that depending on the state of the first 
qubit, we can have |00) and |11) pseudopure state of the remaining qubit. Again, the 
permutations can be carried out by spin-selective pulses and evolution of couplings or 
by transition selective pulses. Here the first qubit is acting as a ancilla qubit and this 
should not be disturbed during the computation. The two subsystems (two circles in 
Fig. 13.5b ) undergo independent and parallel evolution. 

For higher number of qubits, more labeling qubits may be required. However, it is 
not necessary that all the subsystems are in pseudopure states. One clear disadvantage 
of this technique is the requirement of one 'extra' qubit as ancilla. 
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Figure 3.6: Pulse sequence for preparing |00) pseudopure state by using spatial aver- 
aging method in a two qubit system (/i, I2). 6 representing particular pulse angles and 
subscript denoting phase of the pulses. J is the spin-spin coupling constant between 
two spins and Gz denoting pulse field gradient pulses. 



3.2.3 Spatial averaging 



Spatial averaging technique relies on the application of suitable selective pulses and 
magnetic field gradients [l45]|. This method was proposed by Cory et. al. and since 
then it is proven to be one of the popular techniques to prepare pseudopure states. 
Our aim is to prepare the state ll + 1} + 21} 1} from the state 71/J -I- yi^l- Depend- 
ing on the situation whether yi is greater, less, or equal to y2, we will have exact 
values for 9s. For a homonuclear system (71 = 72) or I71/272I < 1, we will have: 
(91,02,93) = in/3,n/4,n/4). On the other hand, suppose I71/272I > 1 e.g., -^^ C 
system, we will have the 9 values as (^i, 92, 9^) = (0, n/l2, 5n/\2). 



A particular case for a homonuclear spin system for the preparation of |00) pseudop- 
ure state is described below by product operator formalism. 
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iG2 

\ll + + '^llUl (3.15) 

The pulse sequence for three spin homonuclear weakly coupled system is also de- 
scribed ll5Tlll32iri33ll . Generalization of spatial averaging method for N-qubits is given 
by sakaguchi et. al. [I134II . 

In addition to these three main approaches , there are few other methods also avail- 
able to prepare pseudopure states such as spatially averaged logical labeling technique 
(SALLT) una, using 'cat states', etc lfT36]l . 

In the next section we propose a different approach that exploits long life-times of 
certain special states called 'singlet states' |[83l[84l . We don't need an extra qubit in this 
method as it was required in logical labeling technique [I137L The following section 
gives a detail description of the theory and pulse sequence required for the preparation 
of pseudopure state by this new approach. Later, experimental demonstrations on model 
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systems consisting of two, three and four-qubit NMR registers are also described. 

3.3 Preparation of pseudopure states using Long-Lived 
Singlet States 

Consider an ensemble of identical molecules each having n spin- 1/2 nuclei in a magnetic 
field. The Zeeman Hamiltonian 'H^ = h Vq/j, is characterized by the frequency of 
Larmor precession, and the z-component of the spin angular momentum operator li of 
spins j = 1 ■ ■ - n [j2|. The eigenstates | + 1/2) of are labeled as |0) and |1) states 
of a qubit, and the multi-spin eigenbasis {|00 • • ■ 00), |00 • • • 01), • • • } is treated as the 
computational basis. 

3.3.1 preparation of singlet states 

The Hamiltonian for an ensemble of spin- 1/2 nuclear pairs of same isotope, in the RF 
interaction frame, can be expressed as 



•K'^" =h 



Av , Ay , , 
—n - —II + JI' ■ f + vnl'/ 



(3.16) 



Here the RF frequency is assumed to be at the mean of the two Larmor frequencies, and 
Av, J and Vi2 correspond respectively to the difference in Larmor frequencies (chemical 
shift difference), the scalar coupling constant and the RF amplitude (all in Hz). 

The detail version of singlet state preparation is shown in Chapter 2.2. The standard 
pulse sequence for the preparation of singlet state is shown in Fig. 13.71 

The propagator form of this pulse sequence upto spin-lock can be written as follow 
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U U U I 

Figure 3.7: Pulse sequence for the preparation and detection of singlet states. All 
the pulses are non-selective pulses. Here ti = 1/(4/), T2 = 1/(4/) + l/(2Av), and 
T3 = l/(4Av), with Av and / being the chemical shift difference (in Hz) and the scalar 
coupling respectively. T4 is the duration of spin-lock. 

The singlet states by themselves are inaccessible to macroscopic observables, but 
can be indirectly detected by removing the equivalence and transforming to observable 
single quantum coherence using the following propagator (see Fig. 13.71) [[83l[84l . 

U]f = ^'^'''-^i^''-'"\ (3.18) 

A more detailed and quantitative analysis of singlet states was carried out using den- 
sity matrix tomography in Chapter 2.3. In the next section we described the preparation 
of pseudopure states starting from this pure singlet states. 

3.3.2 Initializing NMR Registers 
2-qubit register 

The pulse sequence for initializing a 2-qubit NMR register via singlet states is shown in 
Fig B.81 An initially imperfect singlet density matrix gets purified during the spin-lock 
period as a result of the long life time, while the artifact coherences are destroyed by 
relaxation process as well as the inhomogeneities in the spin-lock itself. There exist 
optimal spin-lock conditions at which one obtains singlet states with high fidelity II106[ 

[ml. 

Once the singlet state is prepared with high fidelity, the conversion \Sq^) |01) can 

79 



Chapter 3. Preparation of Pseudopure States Using Long Lived Singlet States 



90o I8O0 9O90 90o I8O0 90i 




Figure 3.8: Pulse sequence for the creation of |01> pseudopure state. Ay and / being 
the chemical shift difference (in Hz) and the scalar coupling respectively. A pulse field 
gradient (G^) used to destroy the unwanted coherences at the end of pulse sequence. 

1 

be easily achieved by the propagator U^' described by 
This propagator work as follows: 



|5o)<5ol = \1- 2l\ll - lllfy - IIUI 
i— 

2/1/2 _ 2/1/2 _ 2/1/2 

X y y X I z 

i90o 



2/1/2 _ 2/1/2 _ 2/1/2 

X z z X y y 



-/l+/2-2/j/2 



i90 



180 



n - If - llllf = 101 )(01 1 (3.20) 



Finally a pulsed field gradient can be used to destroy the residual single and 
multiple quantum coherences generated due to pulse imperfections. If necessary, other 
pseudopure states can be obtained simply by applying NOT gates. 
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Figure 3.9: Circuit diagram for the preparation of a 3 qubit pseudopure state. U^' 
(see Eq. 13.171) and J/j'^ (see Eq. 13.191 ) are unitary propagator as described previously, 
the C'NOT gates with open circles correspond to NOT operation if the control is and 
identity if the control is 1 denoting a pulse field gradient. 



3-qubit register 

For a 3-qubit system (see Fig B.l lb ), we first prepare a two-qubit singlet and apply 
CNOT(2,3), i.e., a NOT gate on qubit-2 controlled by qubit-3. Subsequent spin-lock 
and Ul'^ gate on qubits 1 and 2 initializes a three qubit system into |010) state II137L 

The circuit can be understood as follows. After preparing the singlet on qubits 1 and 
2 (at time point 2 in Fig. 13.91 ). the third qubit remains in a mixed state with a probability 
Po of being in state |0) and a probability pi of being in state |1). Since, 

I01)-_U0) CNOT(23), 101) -110) ^,^, . 

= |0) > = (gi |0), and 

V2 V2 

101) -110) ^,.. CNOT(23), 100) -111) ^,^^ 

(g)|l) > <8)|1), (3.21) 

V2 V2 



the CNOT(2,3) gate transforms the mixed state according to (at time point 3 in Fig. 13.9 



17 19 CNOT(2,3) 

\Sl''){S','-\ ® (polO)(0| + piUXH) i 



Po\Sl''){Sl''\ ® 10X01 + pi|0i'')<0i''| ® |1)<1|. (3.22) 
During the second spin-lock applied on qubits 1 and 2, the singlet part survives, 
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where as the second term consisting of the Bell state = (|00) - |11))/ V2 decays 
fast (at time point 4 in Fig. 13.91) . The singlet \S]^^) is ultimately transformed into |01) 

1 2 

by and thus we obtain |010) pseudopure state with a good approximation (at time 
point 5 in Fig. 13.91) . As mentioned in the introduction, at ordinary NMR conditions 
Po ~ P\, and therefore discarding the second term means a loss of magnetization by a 
factor of 2. Thus the magnetization is halved with every additional qubit in the register. 
This scaling behavior is similar to that of other traditional methods [I143II . 

4-qubit register 

Let us now analyze a 4-qubit register (see Fig. 13.101) . After the first spin-lock, we have 
the singlet state \Sq^) on 1 and 2. The third and fourth qubits are still in mixed states. As 
seen before, if the third qubit is |1), the CNOT(3,2) converts \Sq^) into l^i'^), which will 
be eventually dephased out during the second spin-lock. Therefore we shall consider the 
third spin to be in state |0). Now if the fourth qubit is in state |0), the pseudo-Hadamard 
h(3) gate followed by C'NOT(4,3) (which applies NOT on qubit-4 only if qubit-3 is in 
state |0)) will lead to a correlated state of singlet pairs: 

i^;^') ® 10) ® 10) ^ is^'^) ® ^ ® 10) 

i C'N0T(4, 3) 

\S];^)®\Sl'). (3.23) 

By similar analysis one obtains \S]f')®\4L:'^) if the 4th qubit is originally in state |1). This 
latter spin-order decays fast due to the shorter life-time of 10^^^). Finally the long-lived 
spin-order |5q ^)<8>|5q'^) can be converted into |0101) pseudopure state by the propagators 
and U^^\ 

The fact that only nearest-neighbor interactions are used is highly advantageous in 
practice II137II . Experimentally, the spin-lock of multiple singlet pairs can be achieved 
using sophisticated modulated RF sequences as described in the next section. 
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Figure 3.10: Circuit diagram for the preparation of a 4 qubit pseudopure state. U 



1,2 



(see Eq. 13.171 ). U^'^ (see Eq. 13.191 ). and U^''^ are unitary propagators. U^''* is a similar 

1 2 

propagator as U2 , acting on qubit 3 and 4. G, denoting a pulse field gradient to 
destroy unwanted coherences, the C'NOT gates with open circles correspond to NOT 
operation if the control is and identity if the control is 1. The h-gate corresponds to 

pseudo-Hadamard: |0> i (|0> - |1»/ a/2 and |1) A (|0> + |1»/ V2. 



t3,4 



Initialization for any number of qubits 



Initialization of NMR qubits through long-lived singlet states can be extended to any 
number of qubit in principle II137II . The circuits for initialization of registers with odd 
and even number of qubits are shown in Fig JB.l lb and 13. lib respectively. The basic 
idea is to divide the register into qubit-pairs and prepare a correlated state of singlets. 
Each pair must consist of two qubits of same nuclear species (homonuclear), but dif- 
ferent pairs may be made up of different species (heteronuclear). The correlated singlet 
states can be prepared by using CNOT gates and pseudo-Hadamard gates as shown in 
Fig J3.1 ir b-c). As each pair is converted into a singlet, an RF spin-lock is applied on all 
the singlet pairs. Under the RF spin-lock all the states except the singlet states decay 
rapidly. The circuits in Fig B.l lb and l3.1 lb differ only at the last qubit which is unpaired 
in the odd register. If two or more qubit-pairs are made up of same nuclear species, then 
it might be difficult to selectively spin-lock some of them leaving out others. However, 
as described by the optional spin-locks shown by the dashed boxes in Fig O.l 1[ apply- 
ing spin-locks on pairs which are not yet singlets has little effect on the overall scheme. 
Only exception is for the final spin-lock on an odd-qubit register. 
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Figure 3.11: The circuit diagrams for initializing (a) even-qubit register, and (b) odd- 
qubit register. Here SL and Gz denote spin-lock and the pulsed field gradient. U \ and 
U2 ai"e gates as described the text. In (a) and (b) the C'NOT gates with open circles 
correspond to NOT operation if the control is and identity if the control is 1 . The h- 

gate corresponds to pseudo-Hadamard: |0> (|0> - |1»/ V2 and |1> A (|0> + |1»/ V2. 
The dashed boxes indicate optional spin-locks. 
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which must be applied such that the last unpaired qubit is left undisturbed. It may 
be possible in some cases to overcome even this requirement (see 3-qubit experiment). 
Nevertheless it may be desirable to have the unpaired qubit in an odd-qubit register to 
be of different nuclear species than all others. 

3.4 Experiments 

Strongly modulated pulses are used for designing high fidelity local gates as well as 
CNOT gates [[571I581 . The spin-lock was achieved by WALTZ- 16 - a phase modulated 
RF sequence, which is routinely used in broadband spin decoupling I1106L Spectra 
corresponding to pseudopure states are obtained by linear detection scheme using small 
flip angle RF pulses. Since the diagonal pseudopure states have one energy level more 
populated than the equal distribution in all others, the spectrum should consist ideally 
of only one transition per qubit in each case. Quantitative analysis of the pseudopure 
states are carried out using extended versions of density matrix tomography described 
in chapter 1 (Appendix A) [I106L Finally, the success of the experimental state p in 
achieving a target pseudopure state ppps is measured by calculating the correlation IfSTl . 



Often only the diagonal elements of the density matrices are relevant and in such cases, 
the 'diagonal correlation' can be expressed by replacing all the operators in the above 
expression by their diagonal parts |[58ll . In the following we describe the individual 
cases of two-, three- and four-qubit registers II137L 

3.4.1 2-qubit register 

The two-qubit system, Hamiltonian parameters, and the corresponding pseudopure and 
the reference spectra are shown in Fig J3.12r a-d'). As shown in Fig J3.8[ the experiment 



(Ppps) = 



trace [p • Ppps] 



(3.24) 
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Figure 3.12: The molecular structure (a) and Hamiltonian parameters (b) of 5 - 
bromothiophene-2-carbaldehyde (5 mg dissolved in 0.5 ml of dimethyl sulphoxide-d6), 
forming a homonuclear two-qubit register. In (b) diagonal and off -diagonal elements 
correspond to the chemical shifts and the scalar coupling constant respectively (in Hz). 
The 'H spectra correspond to the pseudopure state (c) and the equilibrium mixed state 
(d). The baiplots (e-h) con^espond to the real (e,g) and imaginary (f,h) parts of theoret- 
ical (e,f) and experimental (g,h) pseudopure |01) state. 
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1 9 

involved preparing singlet using C/j' , followed by RF spin-lock with amplitude 2 kHz 
and duration 12.4 s, which are optimized for high singlet content [I106L The decay 
constant for singlet state was 16.2 s approximately three times the Ti values of the 
two spins. The singlet is then converted into |01) pseudopure state using U^2^. A final 
gradient pulse served to destroy the artifact coherences. The bar plots showing the real 
and imaginary parts of the theoretical and experimental density matrix are shown in 
Fig J3.12r e-h). A very high correlation of 0.995 is obtained with |01) pseudopure stat 

una. 

3.4.2 3-qubit register 

The three-qubit system, Hamiltonian parameters and the corresponding pseudopure and 
the reference spectra are shown in Fig JS.lSr a-d). The decay constant for singlet state 
of spins 1 and 2 was about 18 s, approximately three times of their T\ values. The 
pseudopure state was prepared using the circuit shown in Fig J3.9[ The off-set of the 
spin-lock was at the center of spins 1 and 2. As described in section II-C, the second 
spin-lock is to be applied ideally on qubits 1 and 2 only (without disturbing 3rd qubit), 
which was harder to achieve in this homonuclear system. However, we found that by 
carefully tuning the durations of the spin-locks, we can obtain identical result to that 
of a no spin-lock on the 3rd spin. This was possible due to (i) the cyclic nature of 
the WALTZ- 16 spin-lock, and (ii) the slower decay of |0^) state during the ofF-resonant 
spin-lock compared to that of artifact coherences. The two spin-locks consisted of 500 
Hz WALTZ-16 modulations whose durations were optimized to about 6.3 s. The CNOT 
gate was implemented using a 14 segment strongly modulated RF pulse of duration 
approximately 60 ms and of fidelity 0.96. The bar plots showing the real and imaginary 
parts of the theoretical and experimental density matrix are shown in Fig J3.13r e-h). The 
3-spin tomography is an extension of the technique described in reference II106II and is 
described in appendix B. The correlation of the experimental density matrix with the 
theoretical pseudopure state is |010) is 0.952. The correlation is smaller compared to 
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ii : 1 




Figure 3.13: The molecular structure (a) and Hamiltonian parameters (b) of acryloni- 
trile (5 mg dissolved in 0.5 ml of CDCI3), forming a homonuclear 3-qubit register. In 

(b) diagonal elements are chemical shifts (in Hz) and the off-diagonal elements are the 
scalar coupling constants (in Hz). The 'H spectra correspond to the pseudopure state 

(c) and the equilibrium mixed state (d). The bar plots are showing the real (e,g) and 
imaginary (f,h) parts of theoretical (e,f) and experimental (g,h) pseudopure |010)(010| 
state. 
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Figure 3.14: The molecular structure (a) and Hamiltonian parameters (b) of aspirin 
(5 mg dissolved in 0.5 ml of CD3OD), forming a homonuclear 4-qubit register. In (b) 
diagonal elements are chemical shifts (in Hz) and the off-diagonal elements are the 
scalar coupling constants (in Hz). The barplot in (c) displays the diagonal elements 
of the density matrix obtained by tomography of the pseudopure |1001) state. The 'H 
spectra correspond to the pseudopure state (d) and the equilibrium mixed state (e). 



the two-qubit case, mainly due to the errors in the CNOT gate and the non-selectivity 
of the second spin-lock. Nevertheless, the diagonal correlation is achieved as high as 
0.983 IfBTH . 



3.4.3 4-qubit register 

The four-qubit system, Hamiltonian parameters and the corresponding pseudopure and 
the reference spectra are shown in Fig J3.14[ The pseudopure state was prepared using 
the circuit shown in Fig JS.lOl The singlet decay constants were about 6 s, approximately 
twice the Ti values of the individual spins. We were able to carry out simultaneous spin- 
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da- 

t=0 t > T2 T2 « t « Tlls 



time 

Figure 3.15: Pictorial description of the evolution of PPS over time. At beginning 
(t=0), all the states are possible. After some time {t > T2), most of the states other than 
singlet-PPS decay considerably. After a time period when t :s> T2, only singlet- PPS 
remains while all other states decay close to zero. 

lock of two singlet pairs and initialize a four-qubit register. The two spin-locks were 
achieved by 2 kHz WALTZ- 16 modulations of durations 2 s and 4.5 s each. The two 
CNOT gates were made of 20 segments, approximately 61 ms duration and of fidelities 
about 0.94. The 10 segment h-gate was about 8.2 s long and of fidelity 0.98. Complete 
tomography of a 4-qubit density matrix is a laborious task. After the preparation of the 
pseudopure state, the non-zero quantum ofF-diagonal elements are efficiently destroyed 
by the final gradient pulse. Since only the diagonal elements are of main interest, we 
have carried out the four-qubit diagonal tomography ||58]| . The bar plot showing the 
diagonal part of the experimental density matrix is shown in Fig J3.14[ The diagonal 
correlation is estimated to be approximately 0.97 + 0.01 with |1001) pseudopure state. 
The ffist pair has collapsed to |10) state instead of |01), due to an additional 180 degree 
pulse that was applied on qubits 1 and 2 for refocusing purposes during U^'^ I1137L 

3.5 Conclusions 

An ensemble of nuclear spin-pairs under certain conditions is known to exhibit singlet 
state life-times much longer than other non-equilibrium states. This property of singlet 
state can be exploited in quantum information processing for efficient initialization of 
quantum registers. Here we have described a general method of initialization and ex- 
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perimentally demonstrate it with two-, three-, and four-qubit nuclear spin registers. The 
basic idea is to prepare correlated state of singlets which lives longer than other coher- 
ences and then convert this long-lived state into a desired diagonal pseudopure state (see 
Fig. 13.151) It is hard to initialize proton-based NMR registers using traditional methods. 
As a result many popular NMR registers were based on carbon spins or a combina- 
tion of protons and carbons which permitted initialization by traditional methods. The 
proposed method I1137II of using long-lived is useful since only the nearest neighbor 
interactions are used, the present method may facilitate initialization of larger registers 
with weak long range interactions. Molecules are of interest where in the inter pair 
dipolar couplings are sufficiently weak to keep the singlet states long-lived, while the 
covalent bond mediated scalar interactions among the nearest neighbor spins are suffi- 
ciently strong. The method may also be applicable to registers based on parahydrogens 
which naturally exist in singlet states. Similar techniques may be used for multi-qubit 
initialization in non-NMR systems exhibiting long-lived states. 
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Chapter 4 

Storing Entanglement Via Dynamical 
Decoupling 

In this chapter, we have described the experimental study of dynamical decouplings 
in preserving two-qubit entangled states using an ensemble of spin- 1/2 nuclear pairs in 
solution state. A brief introduction of decoherence and dynamical decouplings are given 
in section 4.1. In section 4.2, we have described Uhrig's dynamical decoupling and its 
usefulness in preserving coherence orders. In section 4.3, we have shown the creation of 
Bell states from long-lived singlet state. In section 4.4, experimental results are shown 
for storing coherence orders in Bell states including singlet state. We found that the 
performance of odd-order Uhrig sequences in preserving entanglement is superior to 
both even-order Uhrig sequences and periodic spin-flip sequences. We also found that 
there exists an optimal order of the Uhrig sequence using which the singlet state can be 
stored at high correlation for about 30 seconds. 

4.1 Introduction 

Harnessing the quantum properties of physical systems have several potential applica- 
tions, particularly in information processing, secure data communications, and quantum 
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simulators [|20ll . It is believed that such quantum devices may play an important role in 
future technology [I144II . But their physical realization is challenging mainly because of 
decoherence - the decay of the coherent states due to interaction with the surrounding 
environment II145[ 114611 . Therefore it is important to minimize the effects of decoher- 
ence using suitable perturbation on the quantum system II147L A technique, known 
as 'dynamical decoupling' involves protecting the quantum states from decoherence by 
driving the system in a systematic manner such that the effective interactions with the 
environment at different instants of time cancel one another. 

In the following subsections we will describe the effect of decoherence and saving 
coherence through dynamical decouplings. 

4.1.1 Decoherence 

As the name suggests, decoherence means loss of coherence order in a system which 
comprises with more than one entities. The idea of decoherence in quantum mechanics 
came much earlier in 1952 [I148II in the context of hidden variable model. It was first 
developed due to the possible explanation of the appearance of wave function collapse 
II148[ 1 14911 ■ Decoherence based explanation started getting acceptance in early 1980 
UlSOi 11511 . However, Decoherence became inevitable tool of understanding quantum 
mechanical interactions in the context of quantum information processing later. Soon it 
became clear that, decoherence is one of the biggest challenges to be overcome in order 
to realize a practical quantum computer. 

Let us consider an isolated 2-level quantum system (a single spin-1/2 system). The 
wave function can be represented by following : 



|iA) = Co|0) + Ci|l), 
with, |coP + |ciP = l. 
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The system can be best represented by the density matrix formalism : 



(01 
(11 



10) ID 



(4.2) 



The diagonal elements are representing the population distribution of the system in 
two states. The off-diagonal elements are the coherence terms. Here, we are interested 
in the evolution of the coherence terms once the system is no more an isolated quantum 
system and is interacting with the environment. The interaction of system-environment 
is a non-unitary process and hence irreversible. Below we see the effect of environment 
on the 2-level super-positioned state. 

mE) = (colO) + cm\E) ^ (co|0)|£o) + cilDl^i)). (4.3) 

Now it can be noticed that the output state is an entangled state and can not be written as 
system and environment separately (unless l^o) = ^"*|£^i))- In terms of density matrix, 
the situation can be represented as below : 

= mE)mE\ 

= Coc;|0)(0| ® \Eo){Eo\ + cic*|l)(l| ® \E,){E,\ 

+cocl\0){l\ ® \Eo){E,\ + cic;|l)(0| ® [EiKEol (4.4) 

Now tracing out the environment from the system gives the necessary information about 
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the system and can be written as : 



p, = tracer [p5£] 



= coc^lOXOI + cic;|l)<l| + {E,\Eo)c^c\m\\ + <£o|£i)cic*|l)<0| 



10) 



ID 



(01 



{Ei\E^)c^c\ 



(4.5) 



(11 (£o|£i)CiC* 



The above equation shows that the coherence terms obtains extra coefficients. Usu- 
ally, when the environment has a large degree of freedom, these coefficients decay ex- 
ponentially with time: 



Hence, after a certain time duration, the coherence terms decay to zero. 
4.1.2 Dynamical decoupling 

Dynamical decoupling is a technique by which it is possible to suppress, at least to some 
extent, the environmental effect on a open quantum system under study. The idea of 
dynamical decoupling has connections to the routinely used NMR decoupling sequences 
where unwanted couplings are averaged out with the applications of suitable modulated 
or unmodulated RF pulse sequences. The dynamical decoupling scheme relies on the 
application of n pulses at certain intervals. Preserving nuclear spin coherences by spin 
flips at regular intervals was long been known in NMR as the famous Carr-Purcell- 
Meiboom-GiU (CPMG) sequence I1152[|153L The CPMG sequence is widely used in 
NMR to measure the transverse relaxation time constants in the presence of spatial 
inhomogeneity of the static magnetic field and temporal fluctuations in the local fields 
arising due to the molecular motion [|2l. The sequence involves a set of N number 
of n pulses uniformly distributed in a duration [0, T] at time instants {t\,t2, ■ ■ ■ Jn]- 



|(^i(Ol^o(0)l = e 



.-Ml) 



(4.6) 
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Assuming instantaneous n pulses, j* time instant is linear in j. 

Of course, in practice the n pulses do have finite duration owing to the limited power of 
electromagnetic irradiation generated by a given hardware. Further, the constant time 
period between these spin flips should ideally be shorter than the correlation time of 
the spin-bath interaction. Even, this delay is limited by the maximum duty-cycle that is 
allowed for the hardware. Dynamical decoupling with such bounded controls have also 
been suggested HI 54111 55111 56111 571 1 . For instance Hao et al. have been able to calculate, 
using a particular type of atomic systems, the maximum delay between spin-flips in 
order to efficiently suppress decoherence due to a bath with a finite cut-off frequency 
[I158L By studying the efficiency of the decoupling as a function of the CPMG period 
often it is possible to extract valuable informations about molecular dynamics and such 
studies are broadly categorized under 'CPMG dispersion' experiments II159L 

Recently in 2007, Uhrig generalized the CPMG sequence by considering an optimal 
distribution {t\,t2, • ■ • ,tN} of N spin flips in a given duration [0, T] of time that provides 
most efficient dynamical decoupling [I160II . Using a simple dephasing model, Uhrig 
proved that the time instants should vary as a squared sine bell: 

tj = T sin^ ( V (4.8) 

UDD works well in systems having a high-frequency dominated bath with a sharp 
cutoff I1161[ I162[ 11631 . On the other hand, when the spectral density of the bath has a 
soft cutoff (such as a broad Gaussian or Lorentzian), the CPMG sequence was found 
to outperform the UDD sequence [UMl [HSl [M [EH [M [ffl [M. Suter and co- 
workers have studied these different regimes and arrived at optimal conditions for the 
dynamical decoupling [I171L 

Recently Agarwal has shown using theoretical and numerical calculations that even 
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entangled states of two-spin systems can be stored more efficiently using UDD [I172L 
Since entangled states play a central role in QIP, teleportation, data encryption, and 
so on, saving entanglement is crucial for the efficient physical realization of quantum 
devices [|20]| . More recently dynamical decoupling on an electron-nuclear spin-pair in a 
solid state system has been shown to prolong the pseudoentanglement lifetime by two 
orders of magnitude II173L 

While much of the experimental efforts have been on testing the loss of coherence 
due to T2 processes, here in this chapter, we presented the first experiments where we 
study not only the loss of coherences, but also the loss of entanglement due to both Ti 
and T2 processes. Though newer sequences have been suggested to decouple both of 
these processes, these are yet to be studied experimentally I1174[ 11751 . We have devel- 
oped experimental techniques where we can prepare Bell states with high fidelity and 
characterize these states with high precision [|106lll37ll . we explore the utility of differ- 
ent dynamical decoupling sequences on systems wherein both Ti and T2 relaxations are 
significant. 



4.2 Uhrig dynamical decoupling 

Uhrig dynamical decoupling (UDD) claims to be more efficient than the CPMG se- 
quence (which serves as the best known decoupling sequence for more than 50 years!) 
in preserving coherence orders. The efficiency of UDD over CPMG can be understood 
by various mathematical approaches. Uhrig HI 6011 explained the efficiency of UDD 
by considering the standard spin-boson model in ohmic bath. This model predicts the 
noise-spectrum with a sharp cut-off. Here we describe the 'filter function analysis' in 
brief for the study of UDD' s efficiency I1162L 
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4.2.1 Efficiency of UDD over CPMG 



We can write the Hamikonian of a system interacting with an environment as lll62[[T76ll . 

9i=^[Q+m]^z (4.9) 

where Q. is the unperturbed part representing the system and /3(t) is the time dependent 
fluctuating part due to environmental interaction. As in ref II176II . the time evolution of 
a superposition state initially oriented along Y under the aff"ect of this Hamiltonian can 
be written as 

where |0) and |1) representing the basis states and /3{f) adding the random phase errors. 
Accumulation of such phases lead towards decoherence. A fundamental technique for 
preserving coherence in NMR is 'Spin-echo' given by Hahn ni77L Spin-echo works 
as a refocusing technique by applying a n pulse in between two exact delays. Hahn 
echo became indispensable tool for coherence reorder and soon it was realized that the 
application of series of n pulses at regular interval would be most efi'ective in order 
to reduce dephasing 111521 1153II . Hahn echo acts as a high pass filter for an arbitrary 
noise spectrum S^((X>) and it neutralize the phase errors by slowly Fourier components 
of jS. Now this one n pulse logic can be extended to multiple n pulses technique as well. 
Application of multiple pulses on a qubit system, leads to coherence state as. 



W{T) 



(ctyXt') 



where, ;^ = - \ ^l-lF(cot)dio. (4.11) 
n Jo ^ 

Here, the filter function F{(jjt) contains all the necessary information regarding the ef- 
ficiency of pulse sequence for preserving coherence against the environment influence 
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Spir). Now, F(cot) can be calculated from 

F{cot) = Mcorf , (4.12) 

where \y„(a)T)\ is the Fourier transform of time domain filter function y„(0- Any mod- 
ification of filter function will give different efficiency power of that particular pulse 
sequence. CPMG sequence having n pulses at regular interval was modified by Uhrig 
by repositioning the n pulses at irregular intervals. Noise reduction is shown to be much 
more efficient for Uhrig sequence than CPMG I11601I162L 

Later, Agarwal has shown that this results of efficient UDD can be generalized for 
an entangled system as well II172II . Here our work mainly focuses on the experimental 
studies of UDD and CPMG on such an entangled states as well as on non-entangled 
states. 

4.3 Preparation of Entanglement 

We study storage of entanglement by dynamic decoupling on a pair of spin- 1/2 nuclei 
using liquid state NMR techniques. The sample consisted of 5 mg of 5-chlorothiophene- 
2-carbonitrile dissolved in 0.75 ml of dimethyl sulphoxide (see Figure |4TI) . The two 
protons of the solute molecule differ in the Larmor frequency by Av = 270.4 Hz and 
have an indirect spin-spin coupling constant of / = 4.1 Hz. The T2 relaxation time 
constants for the two protons are about 2.3 s and the Ti relaxation time constants are 
about 6.3 s. 

4.3.1 Preparation of singlet states 

High fidelity entangled states are prepared via long lived singlet states in a procedure 
described in chapter 2. 

The long-lived nature of singlet states under the equivalence Hamiltonian can be 
used to prepare high-fidelity Bell states. The experiment involves preparing an incoher- 

100 



4.3. Preparation of Entanglement 



100 







-Vi ^ 



' ' I ' ' 
-100 



(Hz) 



Figure 4.1: The NMR spectrum and the molecular structure of 5-chlorothiophene- 
2-carbonitrile. 
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Figure 4.2: NMR pulse sequence to study dynamical decoupling on Bell states. An in- 
coherent mixture of singlet and triplet states is prepared which under spin-lock purifies 
to singlet state. The resulting singlet state can be converted to other Bell states. Then 
dynamical decoupling sequence can be applied and the performance of the sequence 
can be studied by characterizing the residual state using density matrix tomography. 
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ent mixture of singlet and triplet states, 

p(0) = -I' ■ ^ |5o)<5ol - iroXrol (4.13) 

from the equilibrium state ll + by using the pulse sequence shown in Figure 1421 IfSTll . 
During the spin-lock \Tq) state rapidly equilibrates with the other triplet states. On the 
other hand, the decay constant of singlet state |5o) during the spin-lock is much longer 
than the spin-lattice relaxation time constant {T\) (and hence the singlet state is known 
as a long-lived state) [[831 [84)1 . Hence at the end of suitable spin-lock we left out with 
high fidelity singlet states. The goodness of the prepared singlet state is measured by 
the tomographic method as described in Chapter 2.3. The correlation of the singlet state 
is given by, 

<P„>W = , ■^[PW-P-l (4.14) 

ytrace [p(0 J • trace [p^J 

In the following we describe preparation of other Bell states from the singlet state in a 
two-qubit NMR system. 

4.3.2 Preparation of other Bell states from singlet states 

Other Bell states can be obtained easily from the singlet state: 

|So) ^ |«A+) = ^(|01) + |10)), 

V2 

\So) ^ I0-) = ^(|OO)-|11)), 

V2 

|5„) 10^) = ^(100) + 111)). (4.15) 

V2 



The z-rotation in the above propagators can be implemented by using chemical shift 
evolution for a period of 1 /(2Av), and qubit selective x-rotation can be implemented by 
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Figure 4.3: Density matrix tomography of Bell states : Real part of (a) singlet state 
\Sq) = -^QOl) - |10) with correlation 0.99, (b) |i/^+> = -^(|01> + |10» with correlation 

0.99, (c) |0^> = -^(|00> - |11» with coiTclation 0.98, and (d) |0+> - -^(|00> + |11» 
with correlation 0.97 
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using radio frequency pulses HI 3711 . Details of dynamical decoupling on the Bell states 
will be described in the next sections. In order to investigate the decoupling perfor- 
mance, it is necessary to quantify the decay of Bell states with decoupling duration. 
The Bell states by themselves are inaccessible to macroscopic observables, but can in- 
directly be detected transforming to observable single quantum coherences Il83l [84ll . 
Alternatively, a more detailed and quantitative analysis of Bell states may be carried out 
using density matrix tomography as described in Chapter2.3. [I106L We have utilized 
the density matrix formalism for the characterization of the Bell states. The goodness 
of prepared Bell states can be evaluated from the definition of correlation using expres- 
sions similar to (14.141) . The density matrices for all the four Bell-state have been shown 
in figure 1431 We achieved high fidelity Bell states with correlation around 0.99. In the 
following we have shown the experimental implementations of dynamical decoupling 
on such entangled states. 

4.4 Storage of entanglement by UDD 
4.4.1 Different orders of UDD 

As described earlier, the UDD scheme consists of a sequence of spin flips placed at time 
instants given by the expression (14.81) . Instead of applying the Uhrig's formula for the 
entire duration of decoupling, we have applied the formula for a short time interval {T) 
consisting of a small number (A'^) of pulses and then repeating the sequence. Figure l44l 
shows pulse sequences for various orders of Uhrig Dynamical Decoupling (we refer to 
an N-pulse UDD sequence as UDD-N). Note that UDD-1 (and UDD-2) are equivalent 
to CPMG sequences, in which repeating segment consists of [tcpmg --^-tcpmg]- In our 
experiments, tcpmg was set to 2 ms and the duration t„ of the n pulse was 27.2 yus. The 
total duration of UDD-N was set ioT = N{2TceMG +Tn), such that for an extended period 
of time, the total number of n pulses remain same irrespective of the order of UDD. Only 
the distribution of n pulses varies according to the order of UDD. For example, in one 
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Figure 4.4: Pulse sequences for various orders of Uhrig Dynamical Decoupling. Note 
that both UDD- 1 and UDD-2 are equivalent to CPMG. The time instants are calculated 
according to the expression (14. 8! ). with A'^ being the order of UDD and the total period 
T ^Nx 4.0272 ms. 
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second of decoupling, there will be about 250 n pulses in all UDD-N. Our investigation 
thus helps in studying the efficiency of decoupling over a fixed duration of time for a 
given number of n pulses dispersed according to different orders of UDD. 

4.4.2 Performance of UDD over CPMG sequence 

Now we describe the performances of UDD-N on the singlet state which was prepared 
as explained before (see Figure l4r2l) . After applying UDD-N for a fixed duration of time, 
we carried out density matrix tomography and evaluated the correlation of the preserved 
state with theoretical singlet density matrix. The correlations for various orders of UDD 
are displayed in Figure 14.51 As can be seen from the figure, the singlet state can be 
preserved for longer durations by UDD-1 (CPMG) than no-decoupling. It is also clear 
that all even-order UDD sequences result in significant fluctuations in the correlation 
of the singlet state. However, the odd order UDD preserve the singlet state for tens of 
seconds. For example, the correlation of the singlet state under UDD-7 at all the sampled 
time points till 20 seconds is above 0.96. This rather surprising even-odd behavior is 
likely due to the differences in the performances of the even and odd ordered sequences 
against the spatial inhmongeneity of the RF pulses. 

One way to quantify the efficiency of dynamical decoupling under various orders of 
UDD in figure 1431 is by counting the number of time instants in which the correlation 
of the preserved state exceeds a given threshold. The bar plot in Figure 14.61 compares 
the number of time instants during decoupling under various orders of UDD in which 
the correlation of the singlet state exceeded 0.9. It can be seen that there exists an 
optimal order of UDD (for a given Tcpmg and r^r), which performs the most efficient 
decoupling. The optimality may be because of the finite width of the n pulse. In a 
CPMG sequence the n pulses are uniformly dispersed, while in Uhrig sequence the n 
pulses are more crowded at the terminals (beginning and ending) of the sequence. For 
example, if there are too many n pulses, Uhrig 's formula will lead to an overlap of 
pulses. Experimentally, the overcrowding of n pulses may also lead to RF heating of 
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Figure 4.5: Experimental correlations (circles) of singlet state as a function of de- 
coupling duration of various orders of UDD. Also shown in the top-left figure is the 
con^elation decay under no dynamical decoupling (squares). 
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UDD order 

Figure 4.6: The number of time instants at which the correlation exceeded 0.9 for 
various orders of UDD-N. 

the sample and the probe. Thus the performance of the UDD sequence does not grow 
indefinitely with the order of the sequence, but instead will fall beyond a certain order. 
In our experimental setting, we find that UDD-7 is the optimal sequence for storing the 
singlet state. There are recent suggestions for decoupling using finite pulses, however 
these are yet to be studied experimentally Ill79[|180l l. 

4.4.3 Decay of magnetization during various dynamical decouplings 

It can be noticed that the attenuated correlation (expression (14.141) ) is insensitive to 
the decay of the overall magnetization {e in (13.101) ). but simply measures the overlap 
between and the theoretical density matrix \4'){4'\- An alternate method is to monitor 
the decay of magnetization (i.e., e) under dynamical decoupling. 

As already mentioned in Chapter 2, singlet state itself can not be measured directly, 
but can be converted to observable magnetization by using a chemical shift evolution 
for a duration l/(4Av) followed by a (|j ^ ^ pulse. Intensity of the resulting signal as 
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Figure 4.7: The decay of the singlet spin-order measured by converting it into observ- 
able single quantum magnetizations. The decay was studied under CPMG sequence 
(squares) as well as under Uhrig sequence (filled circles). The dashed and the soUd line 
correspond to the exponential fits for CPMG and UDD-7 data points respectively. 

a function of the duration of dynamical decoupling is shown in Figure 14.71 As can be 
seen, UDD-7 is no better than CPMG in preserving the overall spin-order. In fact the 
decay constant for CPMG and UDD-7 are 6.1 s and 5.9 s respectively. 

4.4.4 Efficiency of UDD over CPMG for a non-entangled state and 
various Bell states 

Now we compare the efficiency of the optimal sequence UDD-7 with UDD-1 (CPMG) 
for preserving product state (al + crl) and other Bell states. Figure |4~8] shows the vari- 
ation of correlation of product states and the Bell states as a function of the decoupling 
duration I1178II . Here, after preparing each of the initial state, the dynamical decoupling 
was applied for a fixed duration of time. To monitor the correlation, we have carried out 
the density matrix tomography as described earlier [11061 . In the case of no decoupling, 
we observe a rapid decay of the correlation. The UDD-1 (CPMG) sequence shows some 
improvement in the storage time. However, UDD-7 clearly exhibits much longer stor- 
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age times than the CPMG sequence. The superior peri'ormance of UDD-7 on the singlet 
state compared to other Bell states is presumably because of its antisymmetric property 
described in section II. 

4.5 Conclusions 

In this chapter, we have described the effect of decoherence on the quantum system 
and shows a method to tackle it in practical situations. Dynamical decoupling is a 
method by which we can reduce the effect of environment on the system and ultimately 
increase the coherence time scale of the system II178L CPMG is known to be the best 
known dynamical decoupling sequence both theoretically and experimentally for more 
than 50 years until 2007. In 2007, Uhrig introduced a new sequence where instead 
of applying the n pulses at regular intervals, one needs to apply n pulses at irregular 
intervals synchronizing with a sine-square bell. Theoretically it has been well proved 
that Uhrig dynamical sequence (UDD) performs better than CPMG sequence for saving 
coherence orders of a quantum system. The coherence order for an entangled state is 
also proved to be elongated by the application of UUD sequence compared to CPMG 
sequence. Stroboscopic spin flips have already been shown to prolong the coherence 
times of quantum systems under noisy environments. Uhrig's dynamical decoupling 
scheme provides an optimal sequence for a quantum system interacting with a dephasing 
bath. Several experimental demonstrations have already been verified the efficiency of 
such dynamical decoupling schemes in preserving single qubit coherences. 

Here we have shown the first experimental study of UDD sequence on an NMR sys- 
tem. We have studied the efficiencies of CPMG and UDD sequences on 2-qubit Bell 
states both in terms of magnetization as well as in terms of correlation decay II178II . 
While the Uhrig sequence is no better than CPMG sequence in terms of preserving the 
overall magnetization (or spin order), it clearly outperforms the CPMG sequence in pre- 
serving the correlation of the entangled as well as non-entangled states. We summarize 
three important features: 
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Figure 4.8: Experimental coiTelations of the product state and various Bell states as 
a function of duration under (i) no decoupling (open squares), (ii) CPMG sequence 
(filled circles), and (iii) UDD-7 (open circles). 
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(i) The even-order UDD sequences result is fluctuations in correlations. 

(ii) The odd-order UDD sequences out-perform the CPMG sequence. 

(iii) There exists an optimal length for the odd order UDD sequence which exhibits the 
most efficient decoupling. 

In our case, UDD-7 of 28.2 ms duration appeared to outperform all other sequences of 
both lower and higher orders. Further understanding on the subject can be achieved by 
carrying out investigations into the effects of other experimental issues like RF inho- 
mogeneity, resonance off- set, errors in calibration of pulse angle etc. These considera- 
tions may help in the theoretical and practical understanding of the optimal decoupling 
schemes. 
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Chapter 5 

Violation of Leggett-Garg Inequality 



In this chapter, we have performed experimental implementation of a protocol for testing 
the Leggett-Garg inequality (LGI) for nuclear spins in a NMR setup. The motivation and 
importance of this work is given in the introduction section [STl In section [521 we have 
laid out the mathematical formulation of Leggett-Garg inequality for a spin- 1/2 nucleus 
in external magnetic field. In section 15. 3[ we have presented the Moussa protocol for 
evaluating the expectation values of a target operator using an ancilla qubit. Section [541 
shows the experimental results for the 3-qubit and 4-qubit measurements respectively. 
The conclusion is given in section [531 

5.1 Introduction 

Distinguishing quantum from classical behavior has been an important issue since the 
development of quantum theory [r7n[78l[T8nil82[|204L This issue is also at the heart 
of physical realizations of quantum information processing (QIP) EOll . Experimental 
tests for confirming quantumness in physical systems are usually guided by the Bell- 
type inequalities (BI) [[781 and the Leggett-Garg inequality (LGI) [I204L BI places 
bounds on certain combinations of correlation coefficients corresponding to measure- 
ment outcomes for space-like separated systems which are assumed unable to influ- 
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ence one another {local realism). LGI, on the other hand, places bounds on combi- 
nations of temporal correlation coefficients between successive measurement outcomes 
for a system. Here the system at any instant of time is assumed to be in one or the 
other of many possibie states, and each measurement is assumed to be perfectiy non- 
invasive, in the sense that it has no effect on system's subsequent dynamics {macro- 
realism). In other words, violation of LGI indicates that the system's dynamics can- 
not be understood in classical terms. In recent years various protocols for implement- 
ing LGI and its refined versions have been proposed and experimentally demonstrated 
[fT83l[TM[T85l[TM[T87l[T88l[T89ll . 

Here we have implemented the LGI protocol for individual spin- 1/2 nuclei (from a 
liquid NMR sample) precessing in magnetic field and interacting with their local envi- 
ronments. A typical spin- 1/2 system is genuinely 'microscopic' and exhibits quantum 
behavior. However, it is well-known that, due to decoherence, microscopic quantum 
systems appear to behave classically and as a consequence QIP tasks relying on such 
candidate systems tend to fail [I146L Nuclear spins from an NMR sample are exam- 
ples of microscopic quantum systems that are in constant interaction with their local 
environment and are also candidate systems for QIP tasks. The interactions such as 
dipole-dipole and chemical- shift anisotropy are known to be leading to decoherence, 
dissipation and relaxation processes within the spin ensemble [I190L In experimental 
set-ups such as NMR, successful QIP implementation therefore demands confirmation 
of 'survival' of and determination of 'durability' of quantumness in candidate systems. 
While an LGI test was originally proposed for addressing the fundamental question 
about the ability of a macroscopic system to behave quantum mechanically, consider- 
ing its basic mathematical framework, we extend such a test to investigate survival and 
durability of quantumness within individual nuclear spins interacting with their environ- 
ments. The investigation also sheds light on the possible consistency of the assumptions 
of macrorealism with the 'decoherence perspective' I1191L 
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Although individual nuclear spins from an NMR sample are not directly addressable, 
the sample provides an easily accessible ensemble of nuclear spins from a large number 
of molecules. Therefore the experimental evaluation of a particular temporal correlation 
involves simultaneous implementations of the LGI protocol on a large number of nuclei 
(identical 'targets'). Further, an NMR read-out is an 'ensemble average' obtained in 
terms of magnetization signal. One thus needs to relate the required temporal correlation 
from an LG string with the NMR signal. A quantum network for encoding correlation 
between measurement outcomes of a target system in the phase of a probe system has 
recently been proposed by Moussa et al llTTll . With this network they were able to 
demonstrate quantum contextuality using nuclear spins from a solid state NMR sample. 
In this chapter, we exploit this network for testing LGI. 

Experimental results shown for values of LG-strings containing three and four tem- 
poral correlations as functions of delay between successive measurements llTSl . We have 
found good agreement between the quantum mechanically expected and experimentally 
observed values of the strings for short timescales over which the decay in correlations 
due to typical NMR relaxation processes are ineffective. Further, to demonstrate effect 
of decoherence on the state of individual target nuclei which leads to relaxation of the 
entire ensemble, we have also measured the values of LG strings over longer timescales 
and found that the LG strings gradually decay and ultimately fall within the classical 
bounds. 

5.2 Leggett-Garg inequality 

Consider a system (the 'target') whose state-evolution in time is governed by a particu- 
lar Hamiltonian. To perform an LGI test for the system, a particular system-observable 
(say Q) that can be taken as 'dichotomic', i.e. having two possible states with measure- 
ment outcomes Q = +\, requires to be identified. Next, from a set of 'n' measurement 
instants {t\,t2,h, ...,tn], pairs of instants tt and tj, such that j = i+l, and a pair containing 
the first (? = 1) and the last {j = n) instants are to be chosen. For each such pair, one is 
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then required to perform measurements of Q on the target system at the corresponding 
two instants and obtain outcomes Q{tj) and Q{tj). After repeating these two-time mea- 
surements over a large number of trials (say, A^), one can obtain the two-time correlation 
coefficient (TTCC) C,y for each pair given by the formula: 



where, r is the trial number. Finally, the values of these coefficients are to be substituted 
in the n-measurement LG string given by: 



Each coefficient from the r.h.s. of the above LG string would have a maximum value 
of -1-1 corresponding to perfect correlation, a minimum value of -1 corresponding to 
perfect anti-correlation, and for no correlation. Thus, the upper bound for Kn con- 
sistent with macrorealism comes out to be {n - 2), the lower bound is -n for odd 
n, and -(n - 2) for even n. With these considerations the LGI reads -n < K,, < 
(n - 2) for odd n, and - (n - 2) < < (n - 2) for even n. For example, -3 < ^3 < 1 
and -2<K4< 2. 

5.2.1 Spin- 1/2 precession 

The Zeeman Hamiltonian for the precession of a spin- 1/2 nucleus in a magnetic field 
about z-axis, is given by H = ^ojo-^, with co being the angular precession frequency 
and (T^ the Pauli-z operator. For the present work we choose the Pauli-x operator, i.e. 
(Tv, as the dichotomic observable. The quantum mechanical expression of Cjj for 
measurements on the nucleus is given by [I191II 




(5.1) 



Kn - C12 + C23 -I- C34 -I- .... -I- C( 

n-l)n 



(5.2) 




(5.3) 
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Figure 5.1: The protocols for evaluating ^3 = C12 + C23 - C13 (a) and K4 = Cu + 
C23 + C34 - Ci4 (b). In (a) three independent pairs of measurements are used to evaluate 
TTCCs C12, C23, and C13. Similarly (b) uses four pairs of independent measurements 
to evaluate Cu, C23, C34, and C14. 



In Heisenberg representation one can obtain this relation from: 



(5.4) 



Here, 1^)- e {|0), |1)}, is an eigenstate of the Pauli-z operator. If we divide the total 
duration from ti to t„ into {n - 1) equal intervals of duration At, we can express the LG 
string consistent with equation (3) as 



Kn = (n- l)cos{a»Af} - cos{(n - l)toAt]. 



(5.5) 



The protocols for evaluating ^3 and K4 are illustrated in Fig. 15.11 It can be seen that the 
quantum bounds for ^3 and K4 are [-3, +1.5] and [-2 V2, +2 V2] respectively. 



5.3 Evaluating TTCCs using network proposed by Moussa 
et al 



Suppose that we wish to evaluate correlations between the outcomes of repeated mea- 
surements of two commuting dichotomic unitary observables Si and 5 2 for a target 
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system (T). Consider an ancilla qubit (called 'probe' P) and a unitary transformation for 
the joint system 'T + P' , 

Us =lp® (PJt + ® (P-h. (5.6) 

Here P+ and P are the projectors onto the eigenspace of 5 e {5 1,52}, such that S = 
(P^h - (P_)t. 

Using equation |5.6[ it can be shown that the ensemble measurement of the 'probe' 
gives correlation between successively measured commuting observables of the 'target'. 
For evaluating TTCC's from an LG string, the observable-set for the target qubit is 
d-x{ti), o".v(?y)} and the corresponding unitaries to be applied to the joint (P + T) system 
at different time instants ti < tj are 

U^M,) = lp® P^t,) + (&,)p ® P.(t,). (5.7) 

Here d-x(tq) = P+i^q) - P~itq) and q = i, j for time instants tj and tj. The quantum 
network for implementing these unitaries is shown in Fig. I5.2r a). 

Let the target qubit 'T' be initially prepared according to p. If the probe qubit 'P' is 
initially in one of the eigenstates of the d-^ operator, say |+) = (|0) + 11))/ V2, the density 
matrix of the joint system is given by 

(P)p+T = (l+)<+l)p ® (P)r. (5.8) 
Due to the application of the unitaries (7) the joint density matrix evolves according to: 

(P)p+T U(tj, td(ph^TU\t:, tj) = (p')p+T, (5.9) 

where U(tj,ti) = Ua-,{ti)^a-Mi)- terms of the evolved joint density matrix, the 
probabilities of obtaining +1 outcomes for the Pauli-x measurements on the probe are 
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t=0 
Probe l+> 




100 Hz 

Figure 5.2: Quantum network for the evaluation of TTCCs (a) and the corresponding 
NMR pulse sequence (b). The ensemble was initially prepared according to (p)p®(p)t, 
where (p)p = (1 - 6p)I/2 + 6p|+>(+|, and (p)t = (1 - ej)!!! + erl+X+l. Here ep/r is a 
dimensionless quantity which represents the purity of the initial states. 



given by: 

p{+l) = rrp+T[(p')p+T{(l±)(±l)p<»lT}]. (5.10) 

By tracing over the probe states and using eqns. (15. 71 - 15.91) in eqn. 15.101 one obtains: 

p(±l) = trj[[p^(ti)P^(tj) + P_(tdP^(tj)](p)T]. (5.11) 

The ensemble average of the measurement outcome of joint (P + T) observable is given 
by: 

(((5-.v)p®It) = +M+1)-M-1)- (5.12) 
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Substitution of results I5TTT] in equation [5.121 gives : 

{(&,)p^h) = trT[d-M^x(tj)(p)T] 

= {&,iti)&Atj)) = Qj. (5.13) 

Comparing equations 15.11 and 15.131 it is clear that each TTCC in an LG string can be 
evaluated by applying unitaries (7) to the joint (probe + target) system followed by an 
ensemble measurement of Pauli-x operator on the probe. 

5.4 Experiment 

NMR sample consisted of 2 mg of '''C labeled chloroform ('^CHCls) dissolved in 0.7 
ml of deuterated dimethyl sulphoxide (DMSO). To implement the protocol described 
above, the spin- 1/2 nuclei of ^^C and 'H atoms are treated as the target spin and the 
probe spin respectively. All the experiments are carried out on a Bruker 500 MHz 
spectrometer at an ambient temperature of 300 K. The carbon RF offset was chosen 
such that the ^^C spin precesses at an angular frequency of (x> = 2;r x 100 rad/s under 
the effective longitudinal field in the rotating frame of the RF. The proton RF offset was 
chosen at the resonance frequency of 'H spin. The indirect spin-spin coupling constant 
(J) for these two spins is 217.6 Hz. The spin-lattice (TO and spin-spin (T2) relaxation 
time constants for 'H spin are respectively 4.1s and 4.0 s. The corresponding time 
constants for ^^C are 5.5 s and 0.8 s. 

The NMR pulse sequence for evaluating TTCCs is described in Fig. I5.2r b). Initial 
90 degree y-pulses on both probe and target prepares them in states. All the spin ma- 
nipulations including the C-NOT gates corresponding to Ua-^ operation are realized by 
specially designed strongly modulated pulses Il57l [58]| having Hilbert-Schmidt fidelity 
of over 0.995. These RF pulses are designed to be robust against the RF field inhomo- 
geneity in the range of 90% to 1 10% and static field inhomogeneity in the range of -5 
Hz to -1-5 Hz. The evolution of J-coupling during the intervals between the measure- 
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ments are refocused using n pulses on spin. Collective transverse magnetization of 
the probe spins induce an observable emf on a resonant Helmholtz-type coil which is 
amplified, digitized and stored as the probe signal. Quadrature detection of the probe 
signal enables us to measure the x-component of the probe magnetization as the real 
part of the complex signal. After Fourier transform, the probe signal is fitted to a mixed 
Lorentzian line shape to extract the absorptive content. A reference signal was obtained 
by an identical experiment with Af = 0. The correlation C,/AO was measured at each 
value of A? by normalizing the real part of the probe signal with the reference signal. 
Below, first we will prove the dichotomic nature of nuclear spin observable which is a 
requisite for the experimental verification of LGI violation. Later subsections shows the 
experimental results corresponding to LGI violations. 

5.4.1 Confirmation of dicliotomic nature of x-component of nuclear 
spin observable 

As the first step towards the implementation of any LGI protocol, one needs to identify 
a dichotomic observable for the target system - i.e. having only two possible outcomes 
scalable as +1 - for measurements of which temporal correlations are to be evaluated. 
Although Pauli-spin operators (relevant to systems such as spin- 1/2 nuclei) are routinely 
taken as dichotomic observables in NMR-QIP implementations, LGI test requires en- 
suring that this indeed is the case experimentally, despite the presence of dominant 
couplings of the target nucleus with its environment. 

The 'H and '^C spins in chloroform are coupled by indirect spin-spin interaction (J) 
with a strength of 217 Hz. The Hamiltonian for such a two-spin system in a doubly 
rotating interaction frame can be written as 

nn = hvH(r^/2 + hvccr'^Jl + hJa'^a^JA, 

where vh and vc are the precession frequencies of the two nuclei O. In the present 
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Figure 5.3: The energy level diagram of ^H-^^C system (a) displaying four levels 
coiTesponding to two coupled spin-1/2 particles, NMR spectra of (b) and '^C (c) 
nuclei showing splitting due to mutual interactions. The energy levels of '^C spin 
system after decoupUng spin (d), and the corresponding ^^C spectrum (e). 



experiment v// = Hz and vc = 100 Hz. The energy level diagram of such a system 
is shown in Figure 15.3b . The experimental spectrum of 'H spin consists of two lines 
corresponding to the two eigenstates of the ^^C spin (and vice- versa) (Figure IS^b -c). 
The effect of the probe spin (i.e., 'H) on ^^C can be removed by spin-decoupling. Under 
decoupling, the '^C spectrum displays just a single line (Figure [53b ) corresponding to 
a two-level system (Figure [53b ). 

We have also recorded the real part of the intensity of signal corresponding to x- 
magnetization of '^C spin (proportional to (cr^)), under 'H decoupling, as a function of 
precession duration (Figure [541) . The data clearly fits to a cosine oscillation of single 
frequency. 



Thus, given the above confirmations that '^C spin is indeed a two-level system and 
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Figure 5.4: Intensity of H decoupled C spectrum as a function of time. Offsets 
of the rotating frame are adjusted sucii that has zero precession frequency and ^^C 
has a precession frequency of 100 Hz. The continuous line is the cosine fit to the 
experimental data points (dots). 

the intensity of its signal corresponding x-magnetization has a cosine oscillation with a 
single frequency, we can say that ctx observable used for testing LGI in the present work 
is dichotomic. 



5.4.2 Violation of LGI for 3 measurement case 

The 3 -measurement LG string ^3 = C12 + C23 - C13 was evaluated for a»A? varying from 
to 60n, with A? incremented from to 300 ms in 360 equal steps. The results of the 
experiment are shown in Fig. 15.51 The maximum random errors in these experiments 
were found to be about 0.5%. It is clearly seen that the experimental ^3 data points 
violate the classical limit and hence macrorealism. Fig. 15.5b shows the ^"3 plot for an 
extended duration consisting of 30 periods. It can be observed that the experimental 
values of ^3 gradually decay at a time constant of about 288 ms predominantly due to 
Ti and T2 relaxations and due to inhomogeneities in the magnetic field, thus eventually 
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Figure 5.5: Correlations versus At: Cu (a), C23 (b), and C13 (c). ^3 is plotted for 
the range ojAt e [0, An] (d). Continuous lines are theoretically expecting plots with an 
exponential decay constant and crosses are experimentally achieved results at various 
time points. The horizontal Une in (d) demarcate the boundary between the classical 
and the quantum regimes. 

falling within the classical limit for coAt > 26n (« 42 ms). 



5.4.3 Violation of LGI for 4 measurement case 

Similarly, the 4-measurement LG string K4 was measured for ojAt varying from to 
16;r (i.e., for 8 periods), with At varying from to 80 ms. The results of the experiment 
are shown in Fig. 15.71 Unlike the 3-measurement case, where the classical and quantum 
mechanical lower limits for ^3 values match (i.e., -3), the 4-measurement case displays 
violation of the classical limit both in the positive as well as in the negative sides. Similar 
to the previous case, we observe an exponential decay of K4 with a time constant of 
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Figure 5.6: Decay of ^"3 w. r. t. time: ^3 is plotted for the range coAt e [0, SOn] 
. Continuous lines and crosses are used for theoretical (^^3 with decay) and experi- 
mental values respectively. The theoretical line was obtained by numerically fitting 
the ^^3 function given in (15.51 ) with an exponential decay to the experimental data. The 
horizontal line demarcate the boundary between the classical and the quantum regimes. 

about 324 ms. Decay of LG strings is faster than the measured T2 values of either spins 
mainly because Ti's have been measured using CPMG sequences which suppress the 
effects of static field inhomogeneity and local fluctuating fields. 



5.5 Conclusion 

The present investigation of LGI employs an ensemble of nuclear spins and alleviates 
the need for repeated experiments on single isolated systems iTTTll . Simultaneous imple- 
mentation of controlled operations on target-probe pairs enables evaluation of TTCCs 
and hence plotting of LG strings as functions of two-time measurement delays. The 
plots exhibit both violation and satisfaction of LGI respectively for delays shorter than 
and comparable to the relaxation timescales [iTSll . we qualitatively interpret them as fol- 
lows: For time scales, over which environmental effects on spin states are negligible, 
individual target spins can be taken as isolated quantum systems. The plots do reflect 
this fact in terms of violation of LGL However, the spin-environment interaction tends 
to destroy phase relationship characterizing superposition of quantum states of the target 
nuclear spin. As a result, each member from the ensemble, with its respective environ- 
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Figure 5.7: The individual correlations C12, C23, C34, and C14 are plotted in (a-d) and 
K4 = C12 + C23 + C34 - Ci4 is plotted in (e) for the range ojAt e [0, 16n]. Continu- 
ous lines and crosses are used for theoretical {K4 with decay) and experimental values 
respectively. The theoretical line was obtained by numerically fitting the K4 function 
given in (15.51 ) with an exponential decay to the experimental data. The horizontal line 
demarcate the boundary between the classical and the quantum regimes. 
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ment traced out, begins to appear as if pre-existing in either one of the two states (of 
a spin observable chosen for performing measurements, which is Pauli-x in the present 
work) but not in their superposition. Such a gradual transition from quantum to macro- 
realistic behavior of individual microscopic systems manifests itself in terms of decay of 
TTCCs. This ultimately leads to the satisfaction of LGI. Our experimental results thus 
not only demonstrate initial macrorealism-violating dynamics in genuine microscopic 
systems such as individual nuclear spins, but also bring forward their environment- 
induced emergent macrorealistic behavior, captured in terms of satisfaction of LGI and 
consistent with decoherence mechanism. 
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Chapter 6 

Quantum Delayed- Choice Experiment 



In this chapter, we have discussed Bohr's complementary principle and its implication 
on light quanta and subsequently on quantum systems. After giving a short introduction 
of wave-particle duality in section I6.1[ we discussed the various interferometer that is 
been used to study this strange property in section I6.2[ Then we described the theory 
of recently proposed quantum delayed choice experiment in section [63l In section [64l 
we have shown the experimental approach for the implementation of quantum delayed 
choice circuit in an NMR quantum information processor. The conclusion is given in 
section [631 

6.1 Introduction 

"Is light made up of waves or particles?" has been an intriguing question over past 
many centuries, and the answer remains a mystery even today. The first comprehensive 
wave theory of light was advanced by Huygens HI 9211 . He demonstrated how waves 
might interfere to form a wavefront propagating in a straight line, and he could also 
explain reflection and refraction of light. Soon Newton could explain these properties 
of light using corpuscular theory, in which light was made up of discrete particles II193II . 
The corpuscular theory held over a century till the much celebrated Young's double slit 
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experiment clearly established the wave theory of light [I194II . In the Young's exper- 
iment, a monochromatic beam of light passing through an obstacle with two closely 
separated narrow slits produced an interference pattern with troughs and crests just like 
one would expect if waves from two different sources would interfere. Other properties 
of light like diffraction and polarization could also be explained easily using the wave 
theory. The 20th century developments such as Plank's theory of black-body radiation 
and Einstein's theory of photoelectric effects required quantization of light into pho- 
tons II195[ 11961 . But the question remained whether individual photons are waves or 
particles. Subsequent development of quantum mechanics was based on the notion of 
wave-particle duality [I197II . which was essential to explain the behavior not only of the 
light quanta, but also of atomic and sub-atomic entities IfSOll . 

6.2 Studying wave-particle duality by interferometers 

6.2.1 Mach-Zhender Interferometer 

The wave-particle duality of quantum systems is nicely illustrated by a Mach-Zehnder 
interferometer (MZI) (see Fig. 16. 1|) [|1981I199]| . The intensity of the incident light is kept 
sufficiently weak so that photons enter the interferometer one by one. In the open-setup 
(Fig. 16. Ih ). it consists of a beam-splitter BSl, providing each incoming photon with two 
possible paths, named and 1. A phase-shifter in path-1 introduces a relative phase 
between the two paths. The two detectors DO and Dl help to identify the path traveled 
by the incident photon. Experimental results show that only one of the detectors clicks 
at a time [I200II . Each click can then be correlated with one of the two possible paths 
by attributing particle nature to the photons. Here the phase-shifter has no effect on the 
intensity of the photons measured by either detector, and therefore no interference is 
observed in this setup. 

In the closed-setup (Fig. 16.1b ). the interferometer consists of a second beam-splitter 
BS2, which allows the two paths to meet before the detection. Experimental results 
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again show that only one detector clicks at a time. But much to the astonishment of 
common intuition, the results after many clicks do show an interference pattern, i.e., 
the intensities recorded by each detector oscillates with (f> II200II . Since only one photon 
is present inside the interferometer at a time, each photon must have taken both paths 
in the interferometer and therefore this setup clearly establishes the wave property of 
photons. 

The naive question by the classical mindset is "whether the photon entering the 
interferometer decides to take one of the paths or both the paths depending on the ex- 
perimental setup?". Scientists who believed in a deterministic nature had proposed that, 
unknown to the current experimentalist, there exists some extra information about state 
of the quantum system, which in principle dictates whether the photon should take either 
path, or both the paths I1201L In other words, they assumed some hidden information 
availed by the photon coming out of BSl about the existence or non-existence of BS2. 

6.2.2 Wheeler's delayed-choice experiments 

In order to break this causal link between the two beam-splitters. Wheeler proposed 
a modification in the MZI setup (Fig. 16.1b '). in which the decision to introduce or 
not to introduce BS2 is to be made after the photon has already passed through BSl 
[|202[ 12031 1204II . This way, there is no causal connection between the selection of 
the paths by the photon and the presence of BS2. Although initially considered as a 
'thought-experiment', this proposal has recently been demonstrated by Jacques et al 
[I205L In their experimental setup, the second beam-splitter (RBS) was controlled by a 
random number generator (RNG), that choose to switch the beam-splitter ON or OFF 
after the photon has already passed through BSl. The results of this delayed-choice 
experiment was in agreement with Bohr's complementarity principle |[80l . That is, the 
behavior of the photon in the interferometer depends on the choice of the observable 
that is measured, even when that choice is made at a position and a time such that it is 
separated from the entrance of the photon into the interferometer by a space-like inter- 
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Figure 6.1: Different types of Mach-Zehnder interferometer setups (a-d) and equiva- 
lent quantum circuits (e-h). BSl and BS2 are beam splitters, cf) is phase shifter, DO and 
Dl are detectors. RBS is a beam-splitter switched ON or OFF by a random number 
generator (RNG) and QBS is a beam-splitter which is controlled by a quantum system 
in superposition. In the quantum circuits, H is the Hadamard gate and = e"'"'"'^ is 
used to prepare the state of ancilla qubit. 
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val. Breaking the causal link had no effect on the results of the wave-particle duality, 
thus ruling out the existence of hidden information 11205 L 

6.2.3 Quantum delayed-choice experiments 

Recently, lonicioiu and Temo have proposed a modified version (Fig. 16. Id ) of the 
Wheeler's experiment which not only demonstrates the intrinsic duality, but also shows 
that a photon can have a morphing behavior between particle and wave IISTl . In their 
setup, BS2 is replaced with a beam splitter which is switched OFF or ON depending on 
|0) or |1) state of a two-level quantum system. Using this modification, lonicioiu and 
Temo have been able to discard hidden variable theories which attempt to assign in- 
trinsic wave or particle nature to individual photons even before the final measurement. 
This proposed experiment is named as 'Quantum Delayed-Choice Experiment' IfSTll . 

Using nuclear magnetic resonance (NMR) techniques we study the behavior of a tar- 
get spin- 1/2 nucleus going through a similar situation as that of a photon going through 
an interferometer fi741 . Another spin- 1/2 nucleus acts as an ancilla controlling the sec- 
ond beam-splitter. In section 16.31 we briefly explain the theory and in section 16.41 we 
describe the experimental results. 

6.3 Theory 

In the following we shall use the terminology of quantum information. The two possi- 
ble paths of the interferometer are assigned with the orthogonal states |0) and |1) of a 
quantum bit. The equivalent quantum circuits for the different setups of MZI are shown 
in Figs. I6.ir e-h). Similar circuits have previously been used in 'duality computers' 
II206[ 12071 1208II . In these circuits the Hadamard operator H has the function of the beam 
splitter BSl. It transforms the initial state |0) to the superposition (|0) -i- |1))/ V2 such 
that both |0) and |1) states are now equally probable. The detection operators for the two 
detectors are Do = |0)<0| and Di = |1)<1|. 
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In the open setup (Fig. 16.1b ). the state after the phase shift becomes, \\]/p) = (|0) + 
V2. The intensities recorded by the two detectors are given by the expectation 

values, 

5'p,o = <iApl^o|(Ap) = ^ and 



5p,i = <iAp|£>i|iAp) = ^, (6.1) 



independent of the phase introduced. Therefore no interference can be observed and 
accordingly this setup demonstrates the particle nature of the quantum system. The 
visibility of the interference 

max(5)-min(S) 

V = 7^7— — r-— -, (6.2) 

max(5 ) + min(S ) 

is zero in this case. 

The equivalent quantum circuit for the closed interferometer is shown in Fig. 16. If . 
After the second Hadamard one obtains the state, |i/^w) = cos ||0) - ?'sin||l), up to a 
global phase. The intensities recorded by the two detectors are now. 



5'w,o = <iAwl^ol<Aw) = COS - and 

2 'P 

S^,\ = (lAwl^iliAw) = sm -. (6.3) 

Thus as a function of 0, each detector obtains an interference pattern with visibility 
V = 1 . This setup clearly demonstrates the wave nature of the target qubit. 



In the circuit corresponding to the Wheeler's experiment (Fig. 16.1^ ). the decision to 
insert or not to insert the second Hadamard gate is to be made after the first Hadamard 
gate has been applied. 

Here, we focus on the next modification, that is the quantum delayed-choice experi- 
ment [[8TII . In the equivalent quantum circuit (Fig. 16. Ih ). the second Hadamard gate is to 
be decided in a quantum way. This involves an ancilla spin prepared in a superposition 
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state cosalO) + sinall). This state can be prepared by rotating the initial |0) state of 
ancilla by an angle 2a about j-axis (using operator = )• The second Hadamard 
gate is set to be controlled by the ancilla qubit. If the ancilla is in state |0), no Hadamard 
gate is applied, else if the ancilla is in state 1 1 ), Hadamard gate is applied. The combined 
state of the two-qubit system after the control-Hadamard gate is 

I'Awp.a) = cos Qr|i/rp)|0) + sinar|i/^w)|l), (6.4) 

wherein the second ket denotes the state of ancilla. After tracing out the ancilla, the 
reduced density operator for the system becomes, 

Pwp = cos^ ar|(/^p><iAp| + sin^ a\i^^){i//^\. (6.5) 

Again, the intensity recorded by each detector can be obtained by calculating the expec- 
tation values. For example, the intensity at the detector DO is, 

S^pfiia,^) = tr[DoPwp] 

= tr[Do|iAp)(iApl]cos2a + 

tr[Z)o|iAw)(<Awl] sin^a 
= Sp.oCOS^Or + 5w,osin^Qr 

1 2 2 2 

= - COS a + cos — sin a. (6.6) 

2 2 

It can be immediately seen that the visibility v for the above interference varies as sin^ a. 
When a = 0, the quantum system has a particle nature and when a = n/2, it has a wave 
nature. In the intermediate values of a, the quantum system is morphed in between 
the particle and the wave nature. In the following section we describe the experimental 
demonstration of morphing of a quantum system between wave and particle behaviors. 
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6.4 Experiment 



The sample consisted of '^CHCl3 (Fig. 16. 2h ) dissolved in CDCI3. Here and '^C spins 
are used as the target and the ancilla qubits respectively. The two spins are coupled by 
indirect spin-spin interaction with a coupling constant of 7 = 209 Hz. All the experi- 
ments were carried out at an ambient temperature of 300 K in a 500 MHz Bruker NMR 
spectrometer. 

6.4.1 Open and closed interferometers 

The pulse- sequences corresponding to open and closed setups of MZI are shown in Fig. 
I6.2r b-c). In these cases, the circuits (Fig. I6.ir e-f)) need only a single target qubit and no 
ancilla qubit. Here 'H spin is used as the target qubit, and its interaction with ^^C spin 
is refocused during the MZI experiments. Ideally both of these setups need initializing 
the target qubit to |0) state. In thermal equilibrium at temperature T and magnetic field 
Bq, an ensemble of isolated spin- 1/2 nuclei exists in a Boltzmann mixture, 



e = yfiBo/kT is a dimensionless constant which depends on the magnetogyric ratio 7 
of the spin. At ordinary NMR conditions e ~ 10^^ and therefore Peq is a highly mixed 
state. Since preparing a pure |0) state requires extreme conditions, one can alleviate this 
problem by rewriting the equilibrium state as the pseudopure state 



The identity part does neither evolve under the Hamiltonians, nor does it give raise to 
NMR signals, and is therefore ignored. Thus the single qubit equilibrium state effec- 
tively mimics the state |0). 

In all the cases (Fig. I6.2r b-d)). the first Hadamard gate on the target qubit is followed 



Peq = ^e^^'lOXOI + ^-e 



\m 



(6.7) 




(6.8) 
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Figure 6.2: Molecular structure of chloroform (a) and pulse-sequences (b-d) for differ- 
ent setups of MZI. Figs, (b) and (c) correspond to the open and closed setups respec- 
tively, and (d) corresponds to the quantum delayed-choice experiment. The unfilled 
rectangles ai^e n pulses. Shaped pulses are strongly modulated pulses corresponding to 
Hadamard gate (H), Yq, gate, and control-Hadamard (cH) gate. 7r/2 detection pulses 
are shown in dotted rectangles. J is the coupling constant and r is the phase-shifting 
delay. Gi and G2 are two pulsed-field-gradients for destroying coherences. In (d) two 
separate experiments for and '^C are recorded after applying respective 7r/2 detec- 
tion pulses, peq, Pp = |iAp><<Apl> Pw = li^wXiAwl, and pwp = I'/'wpXiAwpl represent the 
states at different time instants. 
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by the phase shift. A 100 Hz resonance off-set of spin was used to introduce the 
desired phase shift 0(t) = 200:7rT, with the net free-precession delay t. Experiments 
were carried out at 21 linearly spaced values of (p in the range [0, In]. The '^C spin was 
set on-resonance and the 7-evolution during r was refocused with a n pulse on ^''C. 

Unlike the open interferometer (Fig. 16.2b ). the closed interferometer (Fig. 16.2b ') 
has a second Hadamard gate. In both of these cases, the intensity recorded by Dl 
detector corresponds to the expectation value of Dq = |0)(0| operator, which is a di- 
agonal element of the density operator. To measure this element, we destroy all the 
off-diagonal elements (coherences) using a pulsed field gradient (PFG) Gi, followed by 
a {n/2)y detection pulse. The most general diagonal density operator for a single qubit 
is p = ^1-1- ccr,, where c is the unknown constant to be determined. After applying the 
(n/2)y detection pulse, we obtain ^1 -I- ccFx- The corresponding NMR signal is propor- 
tional to c. The experimental NMR spectra for the open and closed setups are shown 
in Fig. 16. 3[ These spectra are normalized w.r.t. equilibrium detection. Since both the 
pathways created by BSl are equally probable in the open MZI, c = and therefore 
spectrum vanishes. On the other hand, because of the second beam-splitter (BS2) in 
closed MZI, c becomes (p dependent, and hence the interference pattern. 

The corresponding intensities 5'p(w),o = c-l- 1/2 are shown in Fig. 16.41 The theoretical 
values from expressions (16.11) and (16.31) are also shown in solid lines. The experimental 
visibility of interference in the particle case is 0.02 and that in the wave case is 0.97. As 
explained in the previous section, the open setup demonstrates the particle nature and 
the closed setup demonstrates the wave nature. 

6.4.2 Quantum delayed-choice experiment 

The circuit for quantum delayed-choice experiment is shown in Fig. 16. Ih and the cor- 
responding NMR pulse- sequence is shown in Fig. 16. 2d . This circuit requires one target 
qubit (^H) and one ancilla qubit ('^C). The equilibrium state of the two-qubit system 
does not correspond to a pseudopure state and therefore it is necessary to redistribute 
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particle 




Figure 6.3: The experimental spectra obtained after the open (top trace) and closed 
(bottom trace) setups of MZI. Each spectrum (pair of lines) corresponds to one of the 
21 linearly spaced values of (p in the range [0, 2n]. 

the populations to achieve the desired pseudopure state. We used spatial averaging tech- 
nique to prepare the pseudopure state HI 2911 

Ppps = ^l + 6'|00)<00|, (6.9) 

where e' is the residual purity. 

All the gates on the target and the ancilla were realized using strongly modulated 
pulses (SMPs) [[57l[58]| . The SMPs were constructed to be robust against RF amplitude 
inhomogeneities, which normally have a distribution of about 10 % about the mean. 
Robust pulses were achieved by calculating the Hilbert-Schmidt fidelity between the 
desired operator and the experimental operator for different possible RF amplitude dis- 
tributions, and then maximizing the average fidelity Il59ll . An average fidelity of over 
0.995 was achieved for each gate. After the control-Hadamard gate, the state of the 
two-qubit system is expressed by the density operator p^p (eqn. 16.51) up to the unit 
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(|) (rad) 

Figure 6.4: The experimental intensities Sp_o (particle) and Sw,o (wave) at various 
values of (p. 

background. 

The interference 5wp,o (in eqn. 16.61) due to the detection operator Dq = |00)(00| can 
be obtained by measuring the first diagonal element of the density matrix, and hence 
complete density matrix tomography is not necessary [11061 . As in the single qubit case, 
we apply a PFG G2 which averages out all the coherences and retains only the diagonal 
part of the density matrix. The most general diagonal density matrix of a two-qubit 
system is of the form 

p = il (g) 1 + CiCTj 1 + C2I <8> (T- + CsCTj ® (T,, (6.10) 

with the unknown constants ci, C2, and C3. 

Recording the target spectrum after a in/2)y pulse on the above state gives two sig- 
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Figure 6.5: The experimental spectra obtained after the quantum delayed choice ex- 
periment with (7r/2)y detection pulse on target (^H) qubit. These spectra are recorded 
with 21 equally spaced values oi(p e [0, 2n] and at different a values (as indicated). In 
each spectrum, only one line is expected due to the preparation of pseduopure state. 
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nals proportional to c\ + C3 and c\ - C3. The spectra of the target qubit at various values 
of (p and a are shown in Fig. I6.5[ The signals obtained after applying a {nl2)y pulse 
on either qubit after preparing the |00) pseudopure state are used to normalize these 
intensities. In each spectrum, the left transition (corresponding to the |0) state of an- 
cilla), vanishes because of the particle nature (similar to the top trace of Fig. 16.31) and 
the right transition (corresponding to the |1) state of ancilla) displays the interference 
pattern because of the wave nature (similar to the bottom trace of Fig. 16.31) . 

Similarly, recording the ancilla spectrum after a {n/2)y pulse gives two signals pro- 
portional to C2 + C3 and C2 - C3. From these four transitions one can precisely determine 
all the three unknowns ci , C2, and C3, and obtain the population S wp,o = 1 /4 + ci + C2 + C3 . 
Calculated experimental intensities 5wp,o are shown in Fig. 16. 6h . The intensities were 
measured for five values of a in the range [0, n/2], and for 21 values of (p in the range 
[0, In]. The theoretical values from expression (16.61) are also shown in solid lines. The 
experimental values were found to have small random errors with a standard deviation 
less than 0.01 . The significant systematic errors are due to experimental limitations such 
as radio-frequency inhomogeneity and spectrometer non-linearities. 

The visibility v calculated at different values of a are plotted in Fig. 16.6b . The 
theoretical visibility varies as sin^ a as explained in the section n. There appears a gen- 
eral agreement between the quantum mechanical predication (solid-line) and the exper- 
iments (symbols). 

6.5 Conclusions 

In this chapter, we have studied the open and closed setups of Mach-Zehnder inter- 
ferometer using nuclear spin qubits, and demonstrated the particle-like and wave-like 
behaviors of the target qubit. Previously NMR interferometer has been used to study 
dipolar oscillations in solid state NMR I1209II and to measure geometric phases in multi- 
level systems [|210[ 121 1112121 . We have reported the first experimental demonstration of 
the quantum delayed-choice experiment using NMR interferometry. 
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71/2 71/4 371/8 71/2 

a (rad) 

Figure 6.6: The intensities S wp,o(a, (p) versus phase (p for different values of a (a) and 
the visibility v versus a (b). The theoretical values are shown in solid lines and the 
experimental results are shown by symbols. 
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Chapter 6. Quantum Delayed-Choice Experiment 



Bohr's complementarity principle is based on mutually exclusive experimental ar- 
rangements. However, the quantum delayed-choice experiment proposed by lonicioiu 
and Temo llSTIl . suggests that we can study the complementary properties like particle 
and wave behavior of a quantum system in a single experimental setup if the ancilla is 
prepared in a quantum superposition. This experiment is the quantum version of the 
Wheeler's delayed-choice experiment. The quantum delayed-choice experiment sug- 
gests a reinterpretation of complementarity principle: instead of complementary exper- 
imental setups, the new proposal suggests complementarity in the experimental data. 

The NMR systems provide perfect platforms for studying such phenomena fJM . In 
our experiments we found a general agreement between the intensities and the visibil- 
ities of the interference with the theoretically expected values. These experiments not 
only confirm the intrinsic wave-particle duality of quantum systems, but also demon- 
strates continuous morphing of quantum systems between wave and particle behavior 
of the target qubit depending on the quantum state of the ancilla qubit. 
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Appendix A 



Density Matrix tomography for a pair 
of spin- 1/2 homonuclear system 



A density matrix describes the statistical state of a quantum system and is considered to 
be the most profound way of representing a quantum state. A pair of spin- 1/2 system can 
be expressed fully by a density martix of 4 x 4 order. In general the size of the density 
matrix depends on the number of spins, such as a n-spin system can be represented fully 
by a 2" x 2" density matrix. The density matrix of any quantum state is a Hermitian 
matrix with trace 1 . In a NMR spin system, the diagonal elements of the density matrix 
represent populations, whereas the off-diagonal elements represent coherence orders. 
The density matrix for a typical two spin- 1/2 system is given below: 



100) 101) 110) 111) 



(001 



Poo SQ SQ DQ 



(Oil 



Poi ZQ SQ 



(A.l) 



(101 



(HI 



V 



Pn 



J 



Here, Poo,Poi,Pio, and Pn are representing the populations in jOO), |01), jlO), |11) 
states respectively. Single, double, and zero-quantum coherences are denoted by S Q, 
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DQ, and ZQ respectively. In NMR only single quantum coherences are directly observ- 
able. Hence to determine all these elements, we need to perform some unitary transfor- 
mations by which we can effectivle transfer other non-SQ terms into a SQ coherence 
and hence can be readout from the output signal. 

The general traceless deviation density matrix consists of 15 independent real num- 
bers: 



P = 



Po 



r-i + IS3 ri -I- isi rs + IS5 
pi re -I- ise r2 + is2 



P2 



r^ + is^ 
-Zio Pi , 



(A.2) 



Each coherence element has a real (r) and a imaginary (s) part in it. The elements below 
the diagonal are determined by Hermitian condition : pji, = p^j. Since, spin systems 
have smaller coupling constants, often it is difficult to readout each of the line seperately 
incase of a dispersive spectra. Here, we devoloped the tomographic technique where we 
need to do total integratation for one spin. This way integration errors for a dispersive 
spectra can be reduced significantly. Four combinations of different readout elements 
can be obtained from the integration of complex line shapes of spin 1 and 2 respectively. 

• Real part of spin 1 : Ri = (ri +7:2) 

• Imaginary part of spin 1 : Si = (S1-1-S2) 

• Real part of spin 2 : R2 = (r3-l-r4) 

• Imaginary part of spin 2 : S2 = (s3-i-S4) 

Now consider an RF sequence with propagator representing by U, that transforms 
the original density matrix p into p' = UpUK Single quantum coherences of p' will lead 
to different linear combinations of various elements in p. Thus, by applying different 
propagators on p, we can measure the values of different linear combinations of various 
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elements of p. The real and imaginary values of the integration of 7th spin in A;th exper- 
iment will be labeled as i?^ and S J respectively. Following six one-dimensional NMR 
experiments were found to be sufficient to tomograph a two-spin density matrix: 

1. 1 

2. 90, 

3. i7-180,.^ 

4. 45,^. 180,- 

5. 45,^.180,-i 

6- 45,^- 180,-^ 

Here 1 is the identity i.e., direct observation without applying any extra pulses. Av and 
7 are the chemical shift difference and the scalar coupling respectively (both in Hz). 
The offset is assumed to be at the center of the two doublets and the RF amplitudes are 
assumed to be much stronger than Av. 

Now, we will operate each of this propagator one by one and read out the four 
detectable element orders. We will follow the spin-operator formalism for mathemical 
treatments, /j defines as spin-operator, j can be x, y, or z depending on the phase of a 
particular pulse and k can be either 1 or 2 depending on the spin, which is going through 
the evolution. E. g. I\ denotes evolution of spin-1 w. r. t. x axis. The single spin 
operator can be written as: 



7, = cr,/2 = - 



'0 1^ 



( 



; Iy = cryl2 = - 



-i 

1 



( 



I, = ajl = - 



1 
-1 



(A.3) 



Where cr„ o-y, cr^ are Pauli matrices. Now, for a two spin system matrix operator is of 
the order of 4 x 4. Hence /_J = 7v ® ^d, where Id denotes identity matrix of 2 x 2. Again 
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l\ = /(/ ® Ix- Similarly l], I^, I}, can be written. Obviously for a two-spin system 

h = ll+lj,Iy=^+I^,L=ll+I^ 

The spin-operator formalism for the 6 tomographic experiments described in detail. 



I. 

The first experiment is the identity operator and hence the density matrix will be same 
and the readout elements are as written above: 

R\ = (ri + r2), R'2 = (ra + r^), S\=isi + S2), S\ = (*i + ^2). (A.4) 

The left hand side of each of these equations are known values (achieved by integrating 
the spectra) and righ hand side are unknowns to be determined. 

11. 

The second unitary operator is a 90x pulse. The matrix form of this operator can be 
written as: 

110 
10 1 
10 1 
110 

Now operating U2 on p can be written a&: p2 = U2 ■ p ■ U\. The readout elements 
from p2 are : 

R\ = {n+r2), Rl = {r^ + r^\ S] = {p^+p,), Sl = ipo+P2). (A.5) 



U2 



1 

2 



148 



III. 



The third experiment is a 1/2J evolution of J-coupling, whereas chemical shift is refo- 
cussed since the offset is set in between the spins. Mathematically the operator form for 
this unitary evolution can be written as: 














-1 +i 


1 








-l-i 





v! 





-1 -i 










-1 +i 












Similarly as above, p-i = Ui, ■ p ■ U\. The readout elements from ps are : 



R\ = {s,-S2\ Rl = {s^-sA S] = {r,-r2\ Sl = {r^-r^). (A.6) 

IV. 

The fourth experiment comprises of a 45,^ pulse and a 1/2 J evolution similar to the 
previous one. Mathematically, the operator can be writen in a similar way as shown 
previously. = e^-iW'l) .e^-i-''-i-^\e^-ii'^-^^^^^^^ Then, p4 = f/4-p-f/|. The 

readout elements from p4 are: 



^{Po -P\- Pi) + ]^ir5 - rg + 5i - 52 - *3 + *4), 
^(A) -P\- Pi) + ]^{r5 - rg - 5i + 52 + *3 - *4), 

1 

— (ri - r2 - 55 + 56), 
V2 



1 



(r3 - r4 - 55 - 56). 



(A.7) 
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V. 

The operator form of the fifth experiment can be written as, t/j = g(~'f ^-^c^c) g(-'-^A)_g(-it-2./i/?) g(-i.;r/4./,,)_ 
Then, p=, = Us ■ p ■ U\. The readout elements from p=, are: 



R\ = —pisi - S2- S5- se), 
. 1 

R2 = — p(*3 - S4- S5 + 56), 

Sl = ^(Po- Pi - P2) + i^^n - r2- r3 + r4-r5- r(,), 

^ 2 = ^(Po -Pi- P2) + ^{-ri + r2 + r3 - r4 - rs - r^). (A. 8) 

VL 

The operator form of the sixth and last experiment can be written as, f/g = e("'t-2-^j^?).e(-'-^-'A)_g(-if ^./.'/f) g(-i.;r/4./t 
Then, p^ = ■ p ■ Ul- The readout elements from are: 



6 1 1 

^1 = -^(Po - Pi - P2) + :^(ri - r2 + r^ - r4 + rs + r(,), 

6 1 1 

^2 = ■^iPo-Pi-P2) + :^i-ri+r2-r3 + r4 + r5 + re), 

S\ = -^(-si + S2- S5 + se), 
1 



Sl = —(s3-S4-S5 + se). (A.9) 
V2 



So, at the end of these 6 experiments, one should have a total 24 linear equations 
(Eqs. A.3-A.8) with 15 unknowns. All these equations can be written together in a 
matrix format, as written below: (see Eq. lA.lOl) (AX = Y, where A is a constant matrix, 
X is an unknown matrix to be calculated and Y is a known matrix). 
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Pq 
Pi 

P2 

n 

ri 
r3 
r4 
rs 
r6 

Sl 
Sl 
S3 
S4 
S5 
S6 



-1 -1 

V2 V2 

-1 1 

V2 V2 



Certainly, this is a over-determined problem since 24 linear equations are to be used 
for 15 unknowns. This redundancy however works to increase the precision of the solu- 
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tion by reducing the condition number of the constraint matrix. For the 24 x 15 matrix 
in equation lA.lOi condition number is about 3.7, meaning the solutions are precise to 
5 significant digits. The equation lA. 101 can be solved either by singular value decom- 
position (SVD) or by Gaussian elimination method (both of which are implemented in 
MATLAB). 
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Appendix B 

Density Matrix tomography for a three 
spin- 1/2 homonuclear system 

The method for three-spin tomography is similar to one that we have described for a two- 
spin system. Since the number of unknowns here for three-spin system is much higher 
(63 unknowns) than two-spin (15 unknowns) system, we need to have more number of 
experiments in order to find out all the unknowns faithfully. The 8x8 general density 
matrix (p) for a three-spin system can be written as follows: 





1000) 


1001) 1010) 


1011) 


1100) 


1101) 


1110) 


1111) 


(0001 


Po 


rg -1- isg rs + iss 


ri3 -1- Z5i3 


ri -1- is I 


ri4 -1- Z5i4 


^5 + iSl5 


r25 + iS25 


(0011 




Pi ri6 + isx(, 


rg -1- Z56 


rn + isn 


ri + is2 


rie + is26 


fl8 + isig 


(0101 




Pi 


rio -1- is 10 


r\9 + isi9 


rn + isn 


r3 + is3 


r20 + is20 


(0111 






Ps 


ri?, + is2?, 


rii + isii 


rii + is22 


r4 + is4 


(1001 








P4 


rn +isu 


r-i + isj 


r23 + IS23 


(1011 










P5 


r24 + iS24 


rs + isg 


(1101 












Pe 


ru + isn 


(1111 














(B.l) 
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The lower triangle of the density matrix can be filled by applying the Hermitian 
property of it {pj^ = p*p. The diagonal elements {Pj, j = ^ 6) are representing pop- 
ulation distributions of the density matrix. Applying the traceless property (or identity 
trace property), one of the unknowns can be reduced. All other off-diagonal elements 
are representing the various coherence orders. Each coherence elements has a real (r) 
and imaginary {s) part in it. Elements rj and sj, with j = 1 ^ 12 representing the real 
and imaginary part of single quantum coherences. Whereas rj and sj, with j = 13 ^ 28 
representing the real and imaginary part of double, triple, or zero quantum coherences. 
Only single quantum coherences are directly accessible in NMR. As described in the 
2-spin tomography method, we have to find suitable unitary transformations which can 
transfer the double, triple, zero, and population orders into single quantum coherences. 
Consider a propagator U, that transforms the original density matrix p into p' = UpW . 
Following 13 unitary transformations were found to be sufficient to tomograph a three- 
spin homonuclear system. 

(1) . 1 

(2) - i 

(3) - Jk 

(4) - i 

(5) . 2^ ' ^^90 

(6) . i; ■ 9045 

(7) - 2k ■ 

(8) . 2k ■ 45o 

(9) . i • 6O45 
(10). ^-45.35 
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(11) . 

(12) . i-^Oo-ifc-Wo 

(13) . ^ . 6090 • • 90i35 

Here 1 represents the identity operator i.e., direct observation without applying any 
extra pulses. Jn, J23, and are the scalar couplings between spin 1 & 2, spin 2 & 
3, and spin 1 & 3 respectively (in Hz). The offset is assumed to be at the center of 
the spin-1 and spin-2 and the RF amplitudes are assumed to be much stronger than Av. 
Hence, all the pulses used are non-selective RF pulses. For this particular case Jn is not 
used since it has very small coupling constant (a small coupling constant leads to larger 
duration of evolution which intum makes the result more error prone). 

Interms of unitary operator the delays (j^ or ■^) can be written as bellow. For 
example, let us take the 2'"' experiment (j^), 

U2 = e''^"'''""^-^ , (B.2) 

where, H^^ and Hj denoting Hamiltonian due to chemical shifts and J-couplings. 

3 33 
Hcs = J] Vili; Hj = J] Y/nJuliH. (B.3) 

A combination of pulses and delays can also be seen as a required tomographic experi- 
ments. A pulse can be easily be written as a unitary transformation as shown in detail in 
previous appendix. For example, we can take the experiment named 6O90 (experiment 
no. 5). The unitary operator for this pulse can be written as. 

Us = (ri¥il^il*^, (B.4) 

We need to apply this unitary operator one by one on the primitive density matrix 
(B. 1). By doing individual integration on each of the transition (12 transition for each 
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experiment) and taking the imaginary values as well, we can get a total 312 linearly 
dependent equations. These equations can be solved by singular value decomposition 
(SVD) method and all the 63 unknowns can be find out. 
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